Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,501 T→A 13.7% intergenic (+56/+31) ilvC → / ← ppiC ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,959,5010TA13.7% 519.8 / 16.7 192intergenic (+56/+31)ilvC/ppiCketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base T (53/101);  new base A (27/2);  total (88/104)
Fisher's exact test for biased strand distribution p-value = 1.89e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

AACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCT  >  NC_000913/3959353‑3959650
                                                                                                                                                    |                                                                                                                                                     
aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                                                                                                                      >  2:3723909/1‑150 (MQ=255)
   gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCGACCCACGCCCACACGGAGTGGGAGAAAACACACAAg                                                                                                                                                   >  1:681763/1‑150 (MQ=255)
   gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCACTACCCCACGGTGAGGGAGAAAACACACAAg                                                                                                                                                   >  1:2295048/1‑150 (MQ=255)
   gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCAACCCGACCCACACCCACCGGGGGAGGGAGGAAACAAACAAg                                                                                                                                                   >  1:473386/1‑150 (MQ=255)
   gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCCACCCACTCCCCCCGGGCGAGGGGGGAAACACACAAg                                                                                                                                                   >  1:3207299/1‑150 (MQ=255)
    aaGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                                                                                                                  <  1:864880/150‑1 (MQ=255)
         aTTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCGCACCCCCACCCACTCCCACACGGGGAGGGGGAAAAAACACAATGCCAtc                                                                                                                                             >  1:1484181/1‑150 (MQ=255)
            cGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                                                                                                                          <  1:2026905/150‑1 (MQ=255)
             gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                                                                         <  2:1327947/150‑1 (MQ=255)
              cAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                                                                                        <  2:2527191/150‑1 (MQ=255)
               aGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                                                                                                                       <  1:550704/150‑1 (MQ=255)
                 ccATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGGGCGCTGATGGCCTCACCCCGACCCTCTACCACAGGGGTGGGGGGAAAACACTAAATTGCTTGTGCTGgag                                                                                                                                     >  1:2188875/1‑150 (MQ=255)
                 ccATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCGACCCACGCCCCCCCGGCGCGGGAGAAAAACAACAAGGCCGTCTCCgggcg                                                                                                                                     >  1:3889264/1‑148 (MQ=255)
                   aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGCGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                   >  2:4011015/1‑150 (MQ=255)
                        aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                                              <  1:1715366/150‑1 (MQ=255)
                            cgCAGGTAGGTAAGAAACTGCGCGGCTAGATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCGACCCTCTCCCCCAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCAGGAGAAGGCCTTgag                                                                                                                          <  2:3782524/149‑3 (MQ=255)
                             gCAGGTAGGTAAGAACCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                                                                                                                         <  2:3098837/150‑1 (MQ=255)
                                ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                                                                                      <  1:3967637/150‑1 (MQ=255)
                                ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                                                                                      <  1:2413931/150‑1 (MQ=255)
                                    ggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCCACCCGCTCCCCCAGGGAGAGGGAGAAAACACGCAAGGCCGTCTCCGGGACCAGGCCTGGCTATTAGTTg                                                                                                                  >  1:527488/1‑150 (MQ=255)
                                      taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCAACCCCACCCTCACCCCCAGGGGGAGGGAGAAAACACTCCAGGCCTTATCCGGGAGAAGGCCGTGCGATTAGTTGCg                                                                                                                >  2:216380/1‑150 (MQ=255)
                                      taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCGACCCGCTCCCACAGGGGGTGGGGGAAAACACTCAAGGCCTTCTCCGGGGGGAGGGCTTGCTATTATTTGCg                                                                                                                >  2:3575307/1‑150 (MQ=255)
                                      taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCACAACCCCACGGGGTGGGGGAAAAAACCCAACGCCTTCACCAGGAGAAGGCCGTGCGATTAttttag                                                                                                                >  1:1479283/1‑147 (MQ=255)
                                       aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCTCACCCCGACCCTCTCCCACAGGGTGTTGGTGTAAACAATCAATGCCTTCTACTGGAGAAGGCCTTGGTATTTGtttggg                                                                                                               >  1:903921/1‑146 (MQ=255)
                                       aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGGCCGAAACCCCACCCACGCCCCCCCGGCGAGGGGGGAAAAACACAAGGCCAGCGCCGGGGGAAGGGCGTGGGATTATTTGAgg                                                                                                               >  2:1768606/1‑150 (MQ=255)
                                        aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTCAGTGCGCGCTGATGGCCACACCCCGACCCGCACCCACAGGGAGTTGGTGGAAACAAAGAAGGCCTTCTCCGGGTGAAGGCCGTGCTATTAGTTTCggg                                                                                                              >  2:1205977/1‑149 (MQ=255)
                                            aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                                                                                                          <  2:718496/150‑1 (MQ=255)
                                               gcgcGGCTATATGACAGATATGAAACGGATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCGCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTGCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAc                                                                                                       <  2:925905/150‑1 (MQ=255)
                                                 gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGCGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGGACAGCCCCt                                                                                                     >  2:2692756/1‑150 (MQ=255)
                                                 gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGTGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTTCTATTAGTTGCGGTTCATCACCt                                                                                                     >  2:1015061/1‑150 (MQ=255)
                                                   ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTa                                                                                                   <  2:1938720/150‑1 (MQ=255)
                                                    gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGGCCACACCCCCAACCACGCCCACAAGGAGGGGGTGAAAAAACACAACGCCGTCGGCTGGAGAAGGCCGGGGTGTTTGTTTGGGGACAGCAccaaaa                                                                                                  >  1:733119/1‑146 (MQ=255)
                                                        tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAGGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgat                                                                                              <  2:1838827/150‑1 (MQ=255)
                                                        tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgat                                                                                              <  1:141924/150‑1 (MQ=255)
                                                         atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCACCACCCCACCCACACCCACAGCGGGAGGGGGAAAAAAAAAAACGCCATCAAAAGGAGAAGGCCGTGCTAATATTTTGGGTTCCCCACCGTAatgaat                                                                                             >  1:3498780/1‑148 (MQ=255)
                                                          tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGAACCACACCCACACGGAGAGGGAGAAAACAAACAAAGCCTTCTCCCGGAGAATGCCTTGCGATTAGTTTCGGTACAGCACCAAAAAGAtgt                                                                                            >  2:841950/1‑150 (MQ=37)
                                                            aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTCAGGGAGAGGGAGAAAACACTCACGGCCGTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTACTGATGTGa                                                                                          <  2:3690866/150‑1 (MQ=255)
                                                              aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGata                                                                                        >  2:3397337/1‑150 (MQ=255)
                                                               gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGatat                                                                                       <  1:469065/150‑1 (MQ=255)
                                                                atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCCCCCCCACCCACACCCCCCGGGGGAGGGGGAAAACCAACAACGCCGGAAACAGGCGATGGCCGTGCGATTTGTTTGGGGACAGCACCATAAAGATGTGATATc                                                                                      >  2:1860763/1‑150 (MQ=255)
                                                                 tatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATAccc                                                                                     <  2:2065075/150‑1 (MQ=255)
                                                                 tatGAAACGTATTGCTGTAGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACACAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc                                                                                     <  2:1442992/150‑1 (MQ=255)
                                                                   agAAACTTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGCGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGa                                                                                   <  1:3491392/149‑1 (MQ=255)
                                                                     aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGTGGGAGAAAACACGCAATGCCTTCTCCTGGTGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAAGATGTGATATCCGAAc                                                                                 >  2:2350258/1‑150 (MQ=38)
                                                                     aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTATCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                                                                                 >  2:3784359/1‑150 (MQ=255)
                                                                     aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCACACCCACCCGGAGAGGGAGAAAAAACACAAGGCCATCTCCGCGAGAAGGCCATGCGATTAGTTGCGGGAAAGCACCAAAATGATGTGATATACGaaa                                                                                 >  2:1096861/1‑149 (MQ=255)
                                                                     aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACAACCCCACCCACACCCACAGGGGGAGGGGGAAAACACACAAGGCCGTCACCGGGGGAAGGCCGTGCGATTAGTTGCGGTACAGCCCCTTAAAGATGTGATATACGaaa                                                                                 >  2:845799/1‑149 (MQ=255)
                                                                     aaaCGTATTGCTGTTGCGGGTTAAGTGAGCGCTGACGCCCTCACCCCGACACACTCCCACACGGAGAGTGAGAAAACACACAAGGCCTTAACCTGGACAAGGCCTTGATATAAGAAGAGGTACAGCACCTTAATGATGTGATATCAGAAc                                                                                 >  2:3301298/1‑150 (MQ=255)
                                                                       aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCATCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCATAATGAAGTGATATCCGAACTg                                                                               >  2:3064371/1‑150 (MQ=255)
                                                                        cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                              <  1:2686623/150‑1 (MQ=255)
                                                                        cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                              <  2:3526195/150‑1 (MQ=255)
                                                                        cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                              <  2:3526608/150‑1 (MQ=255)
                                                                        cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                              <  1:1516388/150‑1 (MQ=255)
                                                                        cGTATTGCTGTTGCGGGTGAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                              <  2:2620491/150‑1 (MQ=255)
                                                                          tagtGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                                                                            <  1:3753676/147‑1 (MQ=255)
                                                                           aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtg                                                                           <  2:280902/150‑1 (MQ=255)
                                                                             tGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGtgtg                                                                         <  2:4030460/150‑1 (MQ=255)
                                                                               cagTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGTGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCAATTAGTTGCGGTACAGCACCTTAATGTTGTGATATCCGAACTGGGTGTGCa                                                                       <  2:395020/148‑1 (MQ=255)
                                                                                tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTACCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTATTGATGTGATATCCGAACTGGGTGTGCAg                                                                      <  1:1040267/150‑1 (MQ=255)
                                                                                  ttGCGGGTTAAGTGCGCGCTGATTCCCTCACCCCGACCCACACCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCATAATGATGTTATAAACGAAAAGGGTGTGCAGCg                                                                    >  2:3741015/1‑150 (MQ=37)
                                                                                  ttGCGGGTTAAGTGCGCGCTGATGCCCTCACCGCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTATTGATGTGATATGCGAACTGGGTGTGCAGCg                                                                    <  2:4118143/150‑1 (MQ=255)
                                                                                      gggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGcc                                                                <  1:1750284/150‑1 (MQ=255)
                                                                                         ttAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGt                                                             <  2:1397322/150‑1 (MQ=255)
                                                                                         ttAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAGTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGt                                                             <  1:1060322/150‑1 (MQ=255)
                                                                                          tAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTc                                                            <  1:2217460/150‑1 (MQ=255)
                                                                                               gcgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                                                       <  2:4054352/150‑1 (MQ=255)
                                                                                                cgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc                                                      <  1:3279180/150‑1 (MQ=255)
                                                                                                 gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                                     <  1:3097148/150‑1 (MQ=255)
                                                                                                 gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                                     <  1:3003653/150‑1 (MQ=255)
                                                                                                 gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                                     <  1:2159105/150‑1 (MQ=255)
                                                                                                     tGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAc                                                 <  1:4102150/150‑1 (MQ=255)
                                                                                                      gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAcc                                                <  1:3397337/150‑1 (MQ=255)
                                                                                                      gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGGACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGGCGGCTCCAGTAcc                                                >  2:501377/1‑150 (MQ=255)
                                                                                                       aTGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCg                                               >  2:2473319/1‑150 (MQ=255)
                                                                                                           ccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACa                                           >  2:3000617/1‑150 (MQ=255)
                                                                                                           ccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACa                                           >  2:3466917/1‑150 (MQ=255)
                                                                                                             tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAag                                         >  2:28474/1‑150 (MQ=255)
                                                                                                                   cGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaaa                                   <  1:1724002/150‑1 (MQ=255)
                                                                                                                       ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc                               <  2:3889264/150‑1 (MQ=255)
                                                                                                                           tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTa                           >  2:3642125/1‑150 (MQ=255)
                                                                                                                             ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCGTTGCTATTAGTTGCGGTCCAGCACCTTAATGATGTGATATCCGAACTGGGGGTGCAGCGGGCCGGTCGCCTCCAGTACCGGACAAGAGAAAACCACtgtctc                         >  1:2715496/1‑145 (MQ=255)
                                                                                                                             ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                         >  2:1116057/1‑150 (MQ=255)
                                                                                                                             ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                         >  1:3730531/1‑150 (MQ=255)
                                                                                                                             ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTAGCTATTAGTAGCGCAACAGCACCATAATGATGTGATATCCGAACTGGGTGTGAAGCGGGCCGGTCGGCTCCAGTAACGGACAAGAGAAAACCACTTTATc                         >  1:1212856/1‑150 (MQ=255)
                                                                                                                             ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCGTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCCGTACCGGACAAGAGAAAACCACTTTATc                         >  2:2387053/1‑150 (MQ=255)
                                                                                                                              cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCg                        >  1:3139824/1‑150 (MQ=255)
                                                                                                                                ctgggagtggtataAAACACTCAAGGCCTGCTCCGGTATATGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATTATGTGATATCCTAAGTGGGTGTGCAGGGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGaa                      <  1:1963992/148‑1 (MQ=255)
                                                                                                                                caGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGATAACCACTTTATCTaa                      <  2:1858936/150‑1 (MQ=255)
                                                                                                                                   ggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGc                   >  1:3524996/1‑150 (MQ=255)
                                                                                                                                   ggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGc                   >  1:2407728/1‑150 (MQ=255)
                                                                                                                                         ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAc             >  1:3866196/1‑150 (MQ=255)
                                                                                                                                         ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAc             >  1:2984991/1‑150 (MQ=255)
                                                                                                                                          gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAcc            >  1:3051217/1‑150 (MQ=255)
                                                                                                                                          gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAcc            >  1:3373890/1‑150 (MQ=255)
                                                                                                                                           agaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCa           <  1:3294175/150‑1 (MQ=255)
                                                                                                                                            gaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAccct          >  1:2455331/1‑148 (MQ=255)
                                                                                                                                            gaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCAt          >  1:2455866/1‑150 (MQ=255)
                                                                                                                                            gaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCAt          >  2:163739/1‑150 (MQ=255)
                                                                                                                                             aaaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATc         >  1:1648785/1‑150 (MQ=255)
                                                                                                                                                    tCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGGGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCt  <  1:1128913/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCT  >  NC_000913/3959353‑3959650

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑