Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,599,240 T→A 5.2% K174I (AAA→ATA)  dapA ← dihydrodipicolinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,599,2400TA 5.2% 322.7 / 2.9 114K174I (AAA→ATA) dapAdihydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  ref base T (52/52);  new base A (5/2);  total (60/54)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

TCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTCAGTAT  >  NC_000913/2599093‑2599387
                                                                                                                                                   |                                                                                                                                                   
tCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTg                                                                                                                                                   <  2:1135226/150‑1 (MQ=255)
tCTGCTGCAAGTTTGCACATATGGGCCATATAACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCAAATTGCATGAAGTCCAGCCAGCTCTCAACATAGAAGCTCAGCATAACAAAATCATCTCAAACCAGCTCTATg                                                                                                                                                   >  2:1124399/1‑150 (MQ=255)
 ctgctgCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGa                                                                                                                                                  <  2:944402/150‑1 (MQ=255)
  tgctgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGAt                                                                                                                                                 >  1:1459499/1‑150 (MQ=255)
     tgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTATCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGCAGTCCAGCGCGCTCGCAGCATCGCCGCTCAGCAGACCAACATCCTCTGAAACAAGCTCTGTGATCTg                                                                                                                                              >  1:944650/1‑150 (MQ=255)
     ggcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCTCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTg                                                                                                                                              <  2:1465128/149‑1 (MQ=255)
      gcCAGTTTGCACATCTGGGCCTTATCACGCGCTGCGACGTTAGCCGTAACGGAAATATCCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTgg                                                                                                                                             <  1:17504/150‑1 (MQ=255)
      gcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTgg                                                                                                                                             >  2:688231/1‑150 (MQ=255)
      gcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTgg                                                                                                                                             <  2:1735675/150‑1 (MQ=255)
      gcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCGCTTTGATCTgg                                                                                                                                             <  2:1530647/150‑1 (MQ=255)
      gcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACACCATGCCCGCCCAATTGCTTGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATATgg                                                                                                                                             >  2:589872/1‑150 (MQ=255)
      gcCAGTTTGCACATCTGCGCCATATCACGCGCTGCGACGTAAGCCCTAACGGAAATAACCCCATGACCGCCCAATTGCAAGAAGTCCAGCGCGCTCTCATCTTAGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTgg                                                                                                                                             <  1:243067/150‑1 (MQ=255)
       ccAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGt                                                                                                                                            <  2:445102/150‑1 (MQ=255)
         aGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCTCCCAATTGCATGAAGGCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGttt                                                                                                                                          <  2:924990/150‑1 (MQ=255)
            ttGCACATCTCGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTaca                                                                                                                                       <  1:353432/150‑1 (MQ=255)
             tGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCAAAGAGCTGGtttcaga                                                                                                                                      <  1:1046184/150‑5 (MQ=255)
              gCACATCTTGGCAATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcg                                                                                                                                     <  2:121587/150‑1 (MQ=255)
               cacaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgc                                                                                                                                    >  1:1082411/1‑150 (MQ=255)
                acaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAAATGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgcg                                                                                                                                   <  2:1691489/150‑1 (MQ=255)
                   tCTGGGCCATATCCCGCGCCGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCGCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAAGCAGCTCTTTGATCTGGTTTACACGCGTTa                                                                                                                                <  2:356310/150‑1 (MQ=255)
                   tCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTa                                                                                                                                >  1:1837385/1‑150 (MQ=255)
                     tGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATAGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGACCAAAATCAACTGAAACCAGCTCTTTGATCTGGTTTAGACGCGTTAAg                                                                                                                              >  1:71876/1‑150 (MQ=255)
                      gggCCATATCACGAGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCACGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAATCCAGCTCTTGGATCTGGTTTACACGCGTTATGt                                                                                                                             >  1:1706461/1‑150 (MQ=255)
                       ggCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGtt                                                                                                                            >  1:1504386/1‑150 (MQ=255)
                       ggCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGtt                                                                                                                            >  1:46680/1‑150 (MQ=255)
                        gCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCAAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTAGATCTGGTTTTCACGCGTTAACTTc                                                                                                                           >  1:273024/1‑150 (MQ=255)
                          cATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTccc                                                                                                                         >  2:145110/1‑150 (MQ=255)
                          cATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTccc                                                                                                                         <  1:1779840/150‑1 (MQ=255)
                          cATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTccc                                                                                                                         >  2:421417/1‑150 (MQ=255)
                                      tGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCttt                                                                                                             <  1:1260252/150‑1 (MQ=255)
                                       gCGACGTTAGCCGTAACGGAAATTACCCCATGCCCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTg                                                                                                            <  2:960030/150‑1 (MQ=255)
                                       gCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTg                                                                                                            <  1:1495125/150‑1 (MQ=255)
                                       gCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTg                                                                                                            >  1:1273547/1‑150 (MQ=255)
                                          aCGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGAtt                                                                                                         <  1:1764148/150‑1 (MQ=255)
                                            gTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTcc                                                                                                       <  1:1167850/150‑1 (MQ=255)
                                             ttAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCTCTCAGCAGTACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCg                                                                                                      <  2:336316/150‑1 (MQ=255)
                                               aGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGat                                                                                                    <  2:1701893/150‑1 (MQ=255)
                                               aGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGat                                                                                                    <  1:812215/150‑1 (MQ=255)
                                                gCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCAGTTGCCTCTTTGATTCCGata                                                                                                   <  2:481437/150‑1 (MQ=255)
                                                   gTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGataata                                                                                                <  2:199971/150‑1 (MQ=255)
                                                      aCGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAttt                                                                                             <  2:766538/150‑1 (MQ=255)
                                                       cGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAtttt                                                                                            >  2:821112/1‑150 (MQ=255)
                                                        ggtaaTAACCCCATGAAAGTCCAATTGCATGAAGTCCAGCGCGCTCGCATCTTCGAATATCATCAGAACAAAATCTTCTGATACCATCTCTTTGATCTGGTTTACTAGCGTTAAGTGCCCTGTTGCCTCTTTGATTCCGATAATAttttt                                                                                           <  1:1124399/147‑1 (MQ=255)
                                                         gAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAtttttt                                                                                          <  1:384253/150‑1 (MQ=255)
                                                         gAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTGTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAtttttt                                                                                          >  2:329700/1‑150 (MQ=255)
                                                          aaaTAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTa                                                                                         <  2:921192/150‑1 (MQ=255)
                                                          aaaTAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTa                                                                                         >  2:244782/1‑150 (MQ=255)
                                                            aTAACCCCATGACCGCCCAATTGCATGATGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACt                                                                                       >  1:259560/1‑150 (MQ=255)
                                                             tAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACtt                                                                                      >  2:1723949/1‑150 (MQ=255)
                                                              aaCCCCATGACCGCCCAATTGCATGAAGTCCAGCTCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACttt                                                                                     >  1:996307/1‑150 (MQ=255)
                                                                ccccATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCg                                                                                   <  2:560489/150‑1 (MQ=255)
                                                                  ccATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCGTCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGcc                                                                                 >  2:741022/1‑150 (MQ=255)
                                                                  ccATGACCGCCCAATTGCATGAAGTCCAGCGCGCTAGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTAAACGCGTTAAGGTCCCTGGTTCCTCTTTGATTCCGATAATATTTTTTACTTTCGcc                                                                                 <  1:163848/150‑1 (MQ=255)
                                                                    aTGACCGCCCAAATGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAg                                                                               <  2:720218/150‑1 (MQ=255)
                                                                        ccGCCCAATTGCATGAAGTCCAGCGCGCTCGCCTCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTGCGATAATATTTTTTACTTTCGCCAGTCgg                                                                           >  1:1514451/1‑150 (MQ=255)
                                                                        ccGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCTCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACgg                                                                           >  1:1543975/1‑150 (MQ=255)
                                                                        ccGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACgg                                                                           >  1:879372/1‑150 (MQ=255)
                                                                           cccAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCGCAGCAGAACAAAATCATCGGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGccc                                                                        <  2:125969/150‑1 (MQ=255)
                                                                             cAATTGCATGAAGTCCAGCTCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGATACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCAc                                                                      <  2:692067/150‑1 (MQ=255)
                                                                                 tGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGtt                                                                  <  1:688231/150‑1 (MQ=255)
                                                                                  gCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGCTCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGttt                                                                 >  1:806993/1‑150 (MQ=255)
                                                                                     ttacggccAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCGGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCg                                                              <  2:335797/144‑1 (MQ=255)
                                                                                      gAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTGTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCGCTTTGATTCCGATAAGAGTTTTTACTTTCGCCAGACGGCCCAACGTTTCCgg                                                             >  1:153862/1‑150 (MQ=255)
                                                                                        aGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGa                                                           >  2:1218240/1‑150 (MQ=255)
                                                                                         gTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAg                                                          >  1:849316/1‑150 (MQ=255)
                                                                                         gTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAg                                                          <  1:851978/150‑1 (MQ=255)
                                                                                           ccAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCa                                                        >  1:1644258/1‑150 (MQ=255)
                                                                                            cAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAg                                                       <  1:404252/150‑1 (MQ=255)
                                                                                            cAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAg                                                       >  1:775688/1‑150 (MQ=255)
                                                                                              gcgagcTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGAt                                                     <  1:45821/150‑1 (MQ=255)
                                                                                                gcgcTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCg                                                   <  1:579988/150‑1 (MQ=255)
                                                                                                gcgcTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCg                                                   <  2:439823/150‑1 (MQ=255)
                                                                                                   cTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCATATCGCAg                                                >  2:1773978/1‑150 (MQ=255)
                                                                                                   cTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAg                                                <  1:294056/150‑1 (MQ=255)
                                                                                                      gCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTTCTTTCGCCAGACGGCCCACCGTATCCGGGAGCAGATCGCAGCCa                                             >  2:107965/1‑150 (MQ=255)
                                                                                                          catcGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTAc                                         >  2:1648715/1‑150 (MQ=255)
                                                                                                          catcGCCGCTAAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCATCTTTGATTCCGATTATATTTTGTACTTTCGCCAGCCGGCCCACCGTGTCCGGGCGCAGATCGCGGCCAGtcc                                         >  2:171196/1‑148 (MQ=255)
                                                                                                           atcGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACg                                        >  2:666986/1‑150 (MQ=255)
                                                                                                             cgccgcTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACggg                                      >  1:1838552/1‑150 (MQ=255)
                                                                                                              gccgcTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGa                                     >  1:1864579/1‑150 (MQ=255)
                                                                                                                 gcTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACgg                                  >  1:1247209/1‑150 (MQ=255)
                                                                                                                    cagcagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGcac                               <  1:1630043/150‑1 (MQ=255)
                                                                                                                      gcagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACAt                             >  2:1314848/1‑150 (MQ=255)
                                                                                                                      gcagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGGTAAGTTCCCTGTTGCCTCTGTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACAt                             >  2:939204/1‑150 (MQ=255)
                                                                                                                          aaCAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATtata                         >  2:754683/1‑150 (MQ=255)
                                                                                                                           aCAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATAc                        <  2:1801590/150‑1 (MQ=255)
                                                                                                                           aCAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGGTAAGTTCCCTTTTGCCTCTTTGATTCCGATAATATTTTTTACTATCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATAc                        <  1:1253066/150‑1 (MQ=255)
                                                                                                                            cAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACa                       <  1:1370806/150‑1 (MQ=255)
                                                                                                                            aaaaaTCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACa                       <  2:1407621/149‑1 (MQ=255)
                                                                                                                                  gcgtgtaaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAAttt                 >  1:191116/7‑150 (MQ=255)
                                                                                                                                  gcgtgtaaaCCAGATCAAAGAGCTGGTTTACACCAGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTAACTTACGCCAGAGCGCCCACCGTTTCCGCGAGCAGAACGCAGACAGTACGGGACGGCACATTATACAGAAatt                 >  1:206322/7‑150 (MQ=255)
                                                                                                                                  catcTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAAttt                 <  2:1082411/150‑1 (MQ=255)
                                                                                                                                    gtgtaaaCCAGATCAAAGAGCTGGTTTAAACGCGTTAAGTTCCCTGTTGCCTCTATGATTCCGATAATATTTTTTACTTTCGACAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGCGACGGCACATTATACAGAATTTGc               >  2:863326/5‑150 (MQ=255)
                                                                                                                                      tGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCgg             >  1:704473/1‑150 (MQ=255)
                                                                                                                                         aaCCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCCCCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGCCGGCACATTATACAGAATTTGCGGCAg          >  2:1875261/1‑150 (MQ=255)
                                                                                                                                         aaCCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAg          >  1:1849583/1‑150 (MQ=255)
                                                                                                                                          aCCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAgg         >  2:254336/1‑150 (MQ=255)
                                                                                                                                              gCTCTTTGATCTGGTTTACACGCGTTACGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGCCGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGCCGGCACATTATACAGAATTTGCggcagggcag     <  1:621201/150‑1 (MQ=255)
                                                                                                                                              gCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTCat     >  1:172087/1‑149 (MQ=255)
                                                                                                                                                 ctTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGACAGTAt  >  2:1111463/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTCAGTAT  >  NC_000913/2599093‑2599387

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: