Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,460,618 C→T 6.8% intergenic (‑149/+32) fadI ← / ← yfcZ beta‑ketoacyl‑CoA thiolase, anaerobic, subunit/UPF0381 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,460,6180CT 6.8% 447.0 / 14.4 133intergenic (‑149/+32)fadI/yfcZbeta‑ketoacyl‑CoA thiolase, anaerobic, subunit/UPF0381 family protein
Reads supporting (aligned to +/- strand):  ref base C (60/64);  new base T (6/3);  total (66/67)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 133 total aligned reads displayed.

AACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTC  >  NC_000913/2460470‑2460766
                                                                                                                                                    |                                                                                                                                                    
          tGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGTCAGTTGCTGGCAggg                                                                                                                                           <  2:2212961/150‑1 (MQ=255)
           gTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAgggg                                                                                                                                          >  2:1792474/1‑150 (MQ=255)
           gTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAgggg                                                                                                                                          >  1:1719071/1‑150 (MQ=255)
           gTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAgggg                                                                                                                                          >  1:1547840/1‑150 (MQ=255)
            tAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGc                                                                                                                                         >  1:1690669/1‑150 (MQ=255)
            tAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGc                                                                                                                                         >  2:1249904/1‑150 (MQ=255)
             aaaaaaaTAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACACCCCTGCCAGCAACTGACAGGGGcg                                                                                                                                        >  1:665214/1‑150 (MQ=255)
               aaaaaTAACTACCGGGGTCTGACGTGTTCATAGTATGAACCTTTTTTTTACTTTGCGCCATGTGGCGAAACGGGACTGGGGGAGCTATGCCCCCGAGAGAACAGAAGACGCGAACCGAACACGCCCCTGCCAGCAACTGACAGGGGcgta                                                                                                                                      <  2:519541/150‑1 (MQ=255)
               aaaaaTAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgta                                                                                                                                      >  1:1717552/1‑150 (MQ=255)
                   ataaatacCGTGGTCTGACCGTTTCATCGTCTTAACCATTTTTTTACATTTCGCAAAGTGGTTAAACGGGATTGTGGGTGCTATGACAAACATATAAAAGAAGAAGCGAAAACACCACGCCCCTGTCAGTTGCTGGCAGGGGcgttttct                                                                                                                                  <  2:1445862/150‑6 (MQ=255)
                   ataaataaAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCGAAGTGGAGAAATGGGATTGTGGGAGCTATGACACTGAGAGTAAAGAAGATGAGAATAGTAATCGCCCGTGCCAGCAAGTGACAGGGGcgtacgta                                                                                                                                  >  2:565977/1‑150 (MQ=255)
                   ataaataaAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgtacgta                                                                                                                                  >  2:2368399/1‑150 (MQ=255)
                   ataaaaaaAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgtacgta                                                                                                                                  >  1:273472/1‑150 (MQ=255)
                           aGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGTAGAAGAGAAAAGAAATCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGagag                                                                                                                          >  1:336555/1‑150 (MQ=255)
                             tGGTCTGACCTGATCTTATTCTTACCCATTTTTTTACTTTTAGCCAAGTGGAGAAATGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGa                                                                                                                        <  2:1781870/150‑1 (MQ=255)
                             gggTCGGACCGGATCATAGTCTAAACCATTTTTTTACATTTCGCCATGTGGAGAAAAGGGAATGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAACACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGa                                                                                                                        <  2:321881/149‑1 (MQ=255)
                                tCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAAtt                                                                                                                     >  1:1668949/1‑150 (MQ=255)
                                     ccTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGc                                                                                                                >  1:727651/1‑150 (MQ=255)
                                        gATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAACAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGa                                                                                                             >  1:501611/1‑150 (MQ=255)
                                        gATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGa                                                                                                             >  1:501601/1‑150 (MQ=255)
                                               gTCTTAACATTTTTTTTAATTTTAGCCACGTGGAGAAAAGGTAATGTGGTAGCTATGACGCTGAGTGAAAAGAATAATAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCt                                                                                                      <  2:870793/150‑1 (MQ=255)
                                                tCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTg                                                                                                     >  2:2317113/1‑150 (MQ=255)
                                                 cTTAACATTTTTTTTAAATTTAGCCAAGTGTAGAAAAGGGAAAGTGGGCGCTATTACACAGAGATAAAAGAAGAAGCGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTgg                                                                                                    <  2:1700095/150‑1 (MQ=255)
                                                     acaaTTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGaaaa                                                                                                <  2:1775383/147‑1 (MQ=255)
                                                       cATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATc                                                                                              >  1:1374974/1‑150 (MQ=255)
                                                       cATTTTTTTAAATTTAGCAAAGTTGAGAAAAGGTAATGTGGGAGCTATGACACAGCGACAAAAGAAGAAGACAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATc                                                                                              <  2:1887578/150‑1 (MQ=255)
                                                       aaTTTTTTAAAATTTAGAAAATTTTATAAAAGGCAAATTGGTAGCTATTACACACACACAAAAGAACAATACAAAACAACACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATc                                                                                              <  2:1887294/149‑1 (MQ=255)
                                                        aTTTTTTTACATTTAGCCAAGGGGAGAACCGGTATTGGGGGCGCTATGACACCGAGAGAACAGAAGAAGACAAACGAACACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCa                                                                                             <  2:880078/150‑1 (MQ=255)
                                                              ttACATTTCGCCAAGTGGAGAAAAGGGACAGTGGGAGCTATGCCACAGAGAGAAAAGAAGAAGAGAAAAGAACACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCAttt                                                                                       <  2:359820/150‑1 (MQ=255)
                                                              ttAAATTTAGCCAAGTGGATAAATGGGAATGTGGGAGCTATTACACAGATATAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCAttt                                                                                       <  2:1873912/150‑1 (MQ=255)
                                                                aCATTTAGCCATGTGGATAAAAGGGAATGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCa                                                                                     <  2:2119884/150‑1 (MQ=255)
                                                                aCATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCa                                                                                     <  2:1952834/150‑1 (MQ=255)
                                                                    ttAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCtt                                                                                 >  2:1524389/1‑150 (MQ=255)
                                                                      aGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAATCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCg                                                                               >  2:70039/1‑150 (MQ=255)
                                                                        ccATGTGGATAAAAGGTAAAGTGGGAGCTATTACACAGAGATAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCa                                                                             <  1:402557/150‑1 (MQ=255)
                                                                        ccACGTGTATAAAAGGTAATGTGGGATCTATCACACAGAGATAAAACAAGAACACAAAACAACACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCa                                                                             <  1:461876/150‑1 (MQ=255)
                                                                        ccAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCa                                                                             <  2:1558682/150‑1 (MQ=255)
                                                                               gagaAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCaaaa                                                                      >  1:1456648/1‑150 (MQ=255)
                                                                                 gaAACGGGAAAGTGGGAGCTATGACACAGCGAGAAAAGAAGAAGACAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGt                                                                    <  2:1153500/150‑1 (MQ=255)
                                                                                  aaaaGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTg                                                                   >  2:1977998/1‑150 (MQ=255)
                                                                                   aaggggAAAGTGGGATCTATGACACAGAGATAAAAGAAGACTAGAAACGAACACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTTACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGa                                                                  <  1:1569416/147‑1 (MQ=255)
                                                                                   aaaGGTAATGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGACAAACGAACACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGa                                                                  <  1:1774553/150‑1 (MQ=255)
                                                                                       ggAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATc                                                              >  2:1734639/1‑150 (MQ=255)
                                                                                           cgTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGata                                                          <  1:2402827/149‑1 (MQ=255)
                                                                                                gAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAg                                                     >  1:953170/1‑150 (MQ=255)
                                                                                                     aTGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGca                                                >  1:949991/1‑150 (MQ=255)
                                                                                                     aTGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGca                                                >  1:2312880/1‑150 (MQ=255)
                                                                                                     aTGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGca                                                >  1:1881523/1‑150 (MQ=255)
                                                                                                        acacAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCacac                                             >  2:1755047/1‑150 (MQ=255)
                                                                                                        acacAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCacac                                             >  2:1828047/1‑150 (MQ=255)
                                                                                                          acaGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCg                                           >  1:1204645/1‑150 (MQ=255)
                                                                                                          acaGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCg                                           >  1:793415/1‑150 (MQ=255)
                                                                                                            agagagAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTc                                         >  2:1036099/1‑150 (MQ=255)
                                                                                                             gagagaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCa                                        <  1:475308/150‑1 (MQ=255)
                                                                                                             gagagaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCa                                        <  2:605144/150‑1 (MQ=255)
                                                                                                             gagagaAAAGAAGAAGACAAAAGAACACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCa                                        <  1:1923099/150‑1 (MQ=255)
                                                                                                               gagaAAAGAAGAAGAGAAAAGAACCCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGa                                      <  1:191731/150‑1 (MQ=255)
                                                                                                                agaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGAc                                     >  2:312235/1‑150 (MQ=255)
                                                                                                                agaAAAGAAGAAGAGAAAAGAAAACGCCCCAGCCAGCAACAGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGAc                                     >  1:980357/1‑150 (MQ=255)
                                                                                                                    aagaagaagAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCtt                                 >  2:2281096/1‑150 (MQ=255)
                                                                                                                       aagaagaGAAAAGAAAACGCCCCTGCCAGCACCTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATATCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAg                              >  2:2241190/1‑150 (MQ=255)
                                                                                                                       aagaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAg                              >  1:182808/1‑150 (MQ=255)
                                                                                                                         taagaGAAAACAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTg                            <  1:2006547/149‑1 (MQ=255)
                                                                                                                         gaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTg                            <  1:2138032/150‑1 (MQ=255)
                                                                                                                            gagaAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGaaa                         <  2:316485/150‑1 (MQ=255)
                                                                                                                            gagaAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGaaa                         <  1:2397842/150‑1 (MQ=255)
                                                                                                                                 accaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCgg                    <  1:1653759/147‑1 (MQ=255)
                                                                                                                                 accaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATGAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCgg                    <  1:441766/147‑1 (MQ=255)
                                                                                                                                 aagaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCgg                    <  2:1570708/150‑1 (MQ=255)
                                                                                                                                 aagaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCgg                    <  1:1214042/150‑1 (MQ=255)
                                                                                                                                  agaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGaa                          >  2:1034852/1‑143 (MQ=255)
                                                                                                                                  agaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGaa                          <  1:1034852/143‑1 (MQ=255)
                                                                                                                                   caaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCgggg                  <  1:1086687/149‑1 (MQ=255)
                                                                                                                                    aaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGt                 >  2:1037978/1‑150 (MQ=255)
                                                                                                                                     aaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTg                >  2:1840088/1‑150 (MQ=255)
                                                                                                                                      taCGCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGa               >  1:1851474/2‑150 (MQ=255)
                                                                                                                                      taCGCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGa               >  1:645287/2‑150 (MQ=255)
                                                                                                                                         gCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATCTCCAGGCGGACACCGTCAGGCTCTTCAGGGCGAGTGGGGGGGAtgt            >  1:1708693/1‑148 (MQ=255)
                                                                                                                                         gCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAttt            >  1:1623690/1‑150 (MQ=255)
                                                                                                                                         gCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAttt            >  2:293971/1‑150 (MQ=255)
                                                                                                                                         gCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAttt            <  1:738851/150‑1 (MQ=255)
                                                                                                                                         gCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAttt            <  1:358105/150‑1 (MQ=255)
                                                                                                                                          ccccTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCGGACTCTTCAGTGAAAGTCGGGGTGAtttt           >  1:1820329/1‑150 (MQ=255)
                                                                                                                                          ccccTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtttt           <  2:2014416/150‑1 (MQ=255)
                                                                                                                                          ccccTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtttt           >  1:1369489/1‑150 (MQ=255)
                                                                                                                                           cccTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTg          >  1:499101/1‑150 (MQ=255)
                                                                                                                                           cccTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTg          >  2:872321/1‑150 (MQ=255)
                                                                                                                                            ncTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGACTTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGc         >  1:677463/2‑150 (MQ=255)
                                                                                                                                            ccTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGc         >  1:1662825/1‑150 (MQ=255)
                                                                                                                                            ccTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGc         >  1:1388720/1‑150 (MQ=255)
                                                                                                                                              tGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTACCGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCaa       <  2:442350/150‑1 (MQ=255)
                                                                                                                                              tGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCaa       <  2:1092374/150‑1 (MQ=255)
                                                                                                                                              tGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCaa       <  2:1092359/150‑1 (MQ=255)
                                                                                                                                              tGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCGTCAGTGAAAGTCGGGGTGATTTTGCaa       <  2:1417284/150‑1 (MQ=255)
                                                                                                                                                   gCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTc  <  2:671917/150‑1 (MQ=255)
                                                                                                                                                   gCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTc  <  2:672173/150‑1 (MQ=255)
                                                                                                                                                   gCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTc  <  1:22595/150‑1 (MQ=255)
                                                                                                                                                   gCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTc  <  1:1984615/150‑1 (MQ=255)
                                                                                                                                                   gCAACTGACAGGGGCGTACGTACCCGAGAGGAAGTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTc  <  2:1837250/150‑1 (MQ=255)
                                                                                                                                                   gCAACTGACAGGGGCGTACGGACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTc  <  1:170286/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTC  >  NC_000913/2460470‑2460766

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