Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,291,194 C→A 40.4% C182F (TGC→TTC)  napC ← quinol dehydrogenase, electron source for NapAB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,291,1940CA40.4% 292.4 / 271.6 240C182F (TGC→TTC) napCquinol dehydrogenase, electron source for NapAB
Reads supporting (aligned to +/- strand):  ref base C (97/46);  new base A (55/41);  total (152/88)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 260 total aligned reads displayed.

AGCGTAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTCATTCT  >  CP009273/2291047‑2291341
                                                                                                                                                   |                                                                                                                                                   
aGCGTAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGcaa                                                                                                                                                   <  2:1393897/150‑1 (MQ=255)
aGCGTAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGaaa                                                                                                                                                   >  2:617558/1‑150 (MQ=255)
 gCGTAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGaaat                                                                                                                                                  <  1:2560587/150‑1 (MQ=255)
  cGTAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGaaatc                                                                                                                                                 <  2:687077/150‑1 (MQ=255)
   gTAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCAGATAAAAACCTGGCTCTACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGaaatca                                                                                                                                                >  1:2048202/1‑150 (MQ=255)
    tAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGcaatcaa                                                                                                                                               >  2:2185007/1‑150 (MQ=255)
    tAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGaaatcaa                                                                                                                                               >  2:183175/1‑150 (MQ=255)
    tAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGaaatcaa                                                                                                                                               <  2:461066/150‑1 (MQ=255)
     aaaTGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGaaatcaat                                                                                                                                              <  1:2484560/150‑1 (MQ=255)
      aaTGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATa                                                                                                                                             >  1:2885013/1‑150 (MQ=255)
      aaTGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATa                                                                                                                                             >  1:2884922/1‑150 (MQ=255)
          aCAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAgg                                                                                                                                         >  1:1342673/1‑150 (MQ=255)
               ccACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGc                                                                                                                                    >  1:809911/1‑150 (MQ=255)
               ccACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGc                                                                                                                                    >  1:20924/1‑150 (MQ=255)
               ccACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGc                                                                                                                                    >  1:1985762/1‑150 (MQ=255)
                 aCTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGccc                                                                                                                                  >  1:3698851/1‑150 (MQ=255)
                  cTAAATAAGGTTCGTTCATCCCGCTCACATAGTAACTCTCTGGCTTCATGCATACCCACGCAATAACCCTGTTAAAAACCTGGATCGACTTCACGCATATCGGGCAGCTTGTGCGCGATCCCTTTCTGGTAATCAATACAGGTTTgaaca                                                                                                                                 <  2:2128268/150‑5 (MQ=255)
                  cTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCa                                                                                                                                 >  2:358489/1‑150 (MQ=255)
                  cTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCa                                                                                                                                 <  2:2235425/150‑1 (MQ=255)
                  cTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCa                                                                                                                                 <  1:426502/150‑1 (MQ=255)
                    aaaTAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATc                                                                                                                               >  1:541734/1‑150 (MQ=255)
                    aaaTAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATc                                                                                                                               >  2:3176644/1‑150 (MQ=255)
                     aaTAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCt                                                                                                                              >  1:2273214/1‑150 (MQ=255)
                        aaGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTc                                                                                                                           >  1:756642/1‑150 (MQ=255)
                         aGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCa                                                                                                                          >  1:3901225/1‑150 (MQ=255)
                               gttcATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCt                                                                                                                    <  2:1523388/150‑1 (MQ=255)
                                ttcATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTg                                                                                                                   <  2:4052596/150‑1 (MQ=255)
                                   aTCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTc                                                                                                                <  1:3238553/150‑1 (MQ=255)
                                    tCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCa                                                                                                               <  1:2680945/150‑1 (MQ=255)
                                      ccGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATg                                                                                                             <  1:3199538/150‑1 (MQ=255)
                                       cGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGc                                                                                                            <  2:3243540/150‑1 (MQ=255)
                                        gCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCa                                                                                                           >  1:2138606/1‑150 (MQ=255)
                                           cacaAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCt                                                                                                        >  1:1780527/1‑150 (MQ=255)
                                           cacaAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCTCCTGGTCATGCATCt                                                                                                        <  2:1186994/150‑1 (MQ=255)
                                           cacaAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCt                                                                                                        <  2:2771218/150‑1 (MQ=255)
                                           cacaAAGTAACTCTCTGGCTTCAAGCAGACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCt                                                                                                        <  2:3332346/150‑1 (MQ=255)
                                            acaAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCtt                                                                                                       <  2:177863/150‑1 (MQ=255)
                                              aaaGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTcg                                                                                                     >  2:2476402/1‑150 (MQ=255)
                                                 gTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCgg                                                                                                  >  1:4502330/1‑150 (MQ=255)
                                                  tAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGc                                                                                                 >  2:167922/1‑150 (MQ=255)
                                                  tAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGc                                                                                                 >  2:2968933/1‑150 (MQ=255)
                                                      tctctGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAAcc                                                                                             <  1:4107602/150‑1 (MQ=255)
                                                      tctctGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAAcc                                                                                             >  1:1295653/1‑150 (MQ=255)
                                                       ctctGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCg                                                                                            >  2:4277156/1‑150 (MQ=255)
                                                       ctctGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCg                                                                                            >  2:4277181/1‑150 (MQ=255)
                                                       ctctGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGAGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCTCCTGGTCATGCATCTTCGCGGCACCCg                                                                                            >  2:222144/1‑150 (MQ=255)
                                                        tctGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGa                                                                                           >  2:1985122/1‑150 (MQ=255)
                                                              ttCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTg                                                                                     >  2:4094132/1‑150 (MQ=255)
                                                              ttCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTg                                                                                     >  2:3483021/1‑150 (MQ=255)
                                                               tCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTgg                                                                                    >  2:3244124/1‑150 (MQ=255)
                                                                cAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTggg                                                                                   <  2:554384/150‑1 (MQ=255)
                                                                 aaGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGc                                                                                  >  1:4576292/1‑150 (MQ=255)
                                                                       aCCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGt                                                                            >  1:3609988/1‑150 (MQ=255)
                                                                             gCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCat                                                                      <  1:1206608/150‑1 (MQ=255)
                                                                              cAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCata                                                                     >  2:3409098/1‑150 (MQ=255)
                                                                                aTAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCAtata                                                                   >  2:3066948/1‑150 (MQ=255)
                                                                                  aaCCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACt                                                                 >  2:2569713/1‑150 (MQ=255)
                                                                                  aaCCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACt                                                                 <  1:1993673/150‑1 (MQ=255)
                                                                                   aCCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTc                                                                >  1:3472886/1‑150 (MQ=255)
                                                                                    cccTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCg                                                               >  2:4235255/1‑150 (MQ=255)
                                                                                      cTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGaa                                                             >  1:2962347/1‑150 (MQ=255)
                                                                                         ttAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGtt                                                          <  2:4203369/150‑1 (MQ=255)
                                                                                          tAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATAATCGAAGTtg                                                         >  1:3882957/1‑150 (MQ=255)
                                                                                             aaaCCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGa                                                      <  1:3851684/150‑1 (MQ=255)
                                                                                              aaCCTGGCTCGACTTCACGCATATCGGGCAGCTAGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGAc                                                     <  2:3840778/150‑1 (MQ=255)
                                                                                               aCCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACa                                                    >  1:2456089/1‑150 (MQ=255)
                                                                                                ccTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAg                                                   >  2:823191/1‑150 (MQ=255)
                                                                                                ccTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAg                                                   <  1:2693728/150‑1 (MQ=255)
                                                                                                  tGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCAGGTCATGCATCTTAGCGGCAACCGATTTCTGGGCGGTTGTATCCATAAACTCGAAGTTGTGACAGtt                                                 >  1:1317448/1‑150 (MQ=255)
                                                                                                   ggCTCGACTTCACGCATATCGTGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTa                                                <  2:442565/150‑1 (MQ=255)
                                                                                                   ggCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTa                                                <  2:3745849/150‑1 (MQ=255)
                                                                                                   ggCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTa                                                <  2:577899/150‑1 (MQ=255)
                                                                                                     cTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACg                                              >  1:2738478/1‑150 (MQ=255)
                                                                                                       cGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGc                                            >  1:2819980/1‑150 (MQ=255)
                                                                                                       cGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGc                                            >  2:350730/1‑150 (MQ=255)
                                                                                                        gACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCa                                           >  1:167589/1‑150 (MQ=255)
                                                                                                        gACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCa                                           >  1:4048657/1‑150 (MQ=255)
                                                                                                           ttCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTc                                        <  1:1415437/150‑1 (MQ=255)
                                                                                                             cACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCt                                      <  1:3121065/150‑1 (MQ=255)
                                                                                                             cACGCATATCGGGCAGCTTGTGCGCTATCCCTTNATGGAAATCNATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCt                                      <  1:1898178/150‑1 (MQ=255)
                                                                                                               cGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGc                                    <  1:3395624/150‑1 (MQ=255)
                                                                                                                gCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCg                                   <  1:332972/150‑1 (MQ=255)
                                                                                                                 cATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGa                                  <  1:2311216/150‑1 (MQ=255)
                                                                                                                   tatCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGt                                <  2:2834083/150‑1 (MQ=255)
                                                                                                                     tCGGGCAGCTTGTGCGCTATCACTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTa                              <  2:2119121/150‑1 (MQ=255)
                                                                                                                      cGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTAt                             >  1:1648835/1‑150 (MQ=255)
                                                                                                                         gCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGt                          <  1:4336781/150‑1 (MQ=255)
                                                                                                                          cAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTc                         >  1:825027/1‑150 (MQ=255)
                                                                                                                            gCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCt                       <  2:1278166/150‑1 (MQ=255)
                                                                                                                                gtgCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCAt                   >  1:3507543/1‑150 (MQ=255)
                                                                                                                                 tgCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATg                  <  2:3017253/150‑1 (MQ=255)
                                                                                                                                  gcgcTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATgc                 >  1:4563422/1‑150 (MQ=255)
                                                                                                                                  gcgcTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATgc                 >  1:983637/1‑150 (MQ=255)
                                                                                                                                   cgcTATCCCTTTCTGGCACTCCATACAGGTTGGCCCATCTTTCACCGCCTGGTCATGCATCGTCGCGGCAACCGATTTCTGGGCGGTTGGATCCATATACTCGAAGTTGTGACTGTTACGGCACGCCTGCGTGTTATTGTCCTTCATgcg                <  2:3345359/150‑1 (MQ=255)
                                                                                                                                     cTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGcgcc              >  2:2694495/1‑150 (MQ=255)
                                                                                                                                       aTCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGcgccgc            >  1:1016810/1‑150 (MQ=255)
                                                                                                                                          ccTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCAc         >  2:2696493/1‑150 (MQ=255)
                                                                                                                                             ttATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTCa      >  1:809620/1‑150 (MQ=255)
                                                                                                                                               aTGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTCAtt    >  1:694310/1‑150 (MQ=255)
                                                                                                                                                tGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTCATTc   >  1:2369788/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCGTAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTCATTCT  >  CP009273/2291047‑2291341

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