Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,950,690 C→A 6.1% K518N (AAG→AAT)  ptrA ← protease III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,950,6900CA 6.1% 872.7 / 26.7 230K518N (AAG→AAT) ptrAprotease III
Reads supporting (aligned to +/- strand):  ref base C (85/131);  new base A (8/6);  total (93/137)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 247 total aligned reads displayed.

ATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACAAAGTAAGC  >  CP009273/2950541‑2950839
                                                                                                                                                     |                                                                                                                                                     
aTCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATttc                                                                                                                                                       <  2:3156540/150‑1 (MQ=255)
aTCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATttc                                                                                                                                                       <  2:2232886/150‑1 (MQ=255)
  cACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATttctt                                                                                                                                                     >  1:3907263/1‑150 (MQ=255)
   aCCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATttcttc                                                                                                                                                    >  1:803574/1‑150 (MQ=255)
      tGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTctct                                                                                                                                                 <  1:2185022/150‑1 (MQ=255)
        aTTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGa                                                                                                                                               <  1:3099346/150‑1 (MQ=255)
             gggCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTaa                                                                                                                                          >  2:3407432/1‑150 (MQ=255)
             gggCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTaa                                                                                                                                          >  2:2889295/1‑150 (MQ=255)
               gcgcTTTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATc                                                                                                                                        <  2:4444480/150‑1 (MQ=255)
                cgcTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCa                                                                                                                                       <  1:333716/150‑1 (MQ=255)
                 gcTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAg                                                                                                                                      <  2:2234174/150‑1 (MQ=255)
                 gcTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAg                                                                                                                                      <  1:2106666/150‑1 (MQ=255)
                   tGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCg                                                                                                                                    >  1:44959/1‑150 (MQ=255)
                   tGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCg                                                                                                                                    >  2:2465731/1‑150 (MQ=255)
                     tCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCgag                                                                                                                                  >  2:765215/1‑150 (MQ=255)
                     tCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCgag                                                                                                                                  <  2:2968765/150‑1 (MQ=255)
                         tGGCTTTCGGATTACGCAAAATCAGGCTTACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAat                                                                                                                              <  2:3350291/150‑1 (MQ=255)
                         tGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAat                                                                                                                              <  2:2808125/150‑1 (MQ=255)
                          ggCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAatc                                                                                                                             >  1:891084/1‑150 (MQ=255)
                           gCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAatca                                                                                                                            <  2:3193402/150‑1 (MQ=255)
                                  gATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGaa                                                                                                                     <  2:3092437/150‑1 (MQ=255)
                                      aCGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAAtata                                                                                                                 <  2:3080920/150‑1 (MQ=255)
                                        gCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAg                                                                                                               <  2:1808212/150‑1 (MQ=255)
                                          aaaaTCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAggg                                                                                                             >  2:4663090/1‑150 (MQ=255)
                                           aaaTCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGt                                                                                                            <  1:4537832/150‑1 (MQ=255)
                                           aaaTCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGt                                                                                                            >  1:4620367/1‑150 (MQ=255)
                                            aaTCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGtt                                                                                                           <  2:4335100/150‑1 (MQ=255)
                                             aTCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTa                                                                                                          >  2:4448799/1‑150 (MQ=255)
                                             aTCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTa                                                                                                          <  2:1329805/150‑1 (MQ=255)
                                              tCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTaa                                                                                                         <  2:2507787/150‑1 (MQ=255)
                                              tCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTaa                                                                                                         <  2:196169/150‑1 (MQ=255)
                                               cAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAg                                                                                                        <  2:1050827/150‑1 (MQ=255)
                                               cAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAg                                                                                                        <  2:410620/150‑1 (MQ=255)
                                                aGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGc                                                                                                       >  1:2746396/1‑150 (MQ=255)
                                                aGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGc                                                                                                       >  2:2189045/1‑150 (MQ=255)
                                                 ggCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGct                                                                                                      <  1:2235850/150‑1 (MQ=255)
                                                   cTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGctct                                                                                                    >  2:1951885/1‑150 (MQ=255)
                                                    tGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTg                                                                                                   >  2:1480400/1‑150 (MQ=255)
                                                        aTCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCaa                                                                                               >  1:4067335/1‑150 (MQ=255)
                                                          cAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAag                                                                                             >  2:3592138/1‑150 (MQ=255)
                                                            gCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAagag                                                                                           <  2:1154132/150‑1 (MQ=255)
                                                              tttGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAagagag                                                                                         <  1:773701/150‑1 (MQ=255)
                                                              tttGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAagagag                                                                                         <  2:351840/150‑1 (MQ=255)
                                                                  ggCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCa                                                                                     <  2:1001625/150‑1 (MQ=255)
                                                                   gCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCaa                                                                                    <  2:2251451/150‑1 (MQ=255)
                                                                    cTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAAt                                                                                   <  1:3999501/150‑1 (MQ=255)
                                                                    cTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAAt                                                                                   <  2:2392837/150‑1 (MQ=255)
                                                                      cGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGt                                                                                 >  1:425844/1‑150 (MQ=255)
                                                                        cTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTcg                                                                               <  1:3853548/150‑1 (MQ=255)
                                                                        cTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTcg                                                                               >  1:1472923/1‑150 (MQ=255)
                                                                          ggCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTcggc                                                                             >  1:3368244/1‑150 (MQ=255)
                                                                             aaaaTAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCggcggc                                                                          >  2:3835037/1‑150 (MQ=255)
                                                                             aaaaTAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCggcggc                                                                          >  1:3927386/1‑150 (MQ=255)
                                                                                aTAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCttt                                                                       <  2:3774275/150‑1 (MQ=255)
                                                                                 tAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCtttt                                                                      <  2:1417020/150‑1 (MQ=255)
                                                                                  aaCGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCttttt                                                                     >  1:3748437/1‑150 (MQ=255)
                                                                                   aCGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCtttttt                                                                    <  2:3642941/150‑1 (MQ=255)
                                                                                     ggCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCt                                                                  <  1:1376842/150‑1 (MQ=255)
                                                                                     ggCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCt                                                                  <  1:3051880/150‑1 (MQ=255)
                                                                                     ggCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCt                                                                  <  1:3400442/150‑1 (MQ=255)
                                                                                      gCTTGGCGCATACACCACGCGCAGATTCGACTCATAAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTg                                                                 <  1:4364601/150‑1 (MQ=255)
                                                                                        ttGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGcc                                                               <  2:3825372/150‑1 (MQ=255)
                                                                                               aTACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGc                                                        >  1:1652392/1‑150 (MQ=255)
                                                                                               aTACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGc                                                        >  2:3882262/1‑150 (MQ=255)
                                                                                                     cacGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGt                                                  <  1:3896522/150‑1 (MQ=255)
                                                                                                       cgcgCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGttt                                                >  1:1172101/1‑150 (MQ=255)
                                                                                                        gcgcAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTtg                                               >  2:2947819/1‑150 (MQ=255)
                                                                                                        gcgcAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTtg                                               >  2:576408/1‑150 (MQ=255)
                                                                                                         cgcAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTtgt                                              >  2:648403/1‑150 (MQ=255)
                                                                                                         cgcAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTtgt                                              <  2:2081969/150‑1 (MQ=255)
                                                                                                         cgcAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTtgt                                              <  1:1952472/150‑1 (MQ=255)
                                                                                                          gcAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTtgtg                                             >  1:1309462/1‑150 (MQ=255)
                                                                                                           cAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTgc                                            >  2:2584609/1‑150 (MQ=255)
                                                                                                                tCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGa                                       >  1:3162657/1‑150 (MQ=255)
                                                                                                                    cTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGAtttt                                   <  2:672854/150‑1 (MQ=255)
                                                                                                                     tcatcaACAATCGGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTa                                  <  1:1637087/150‑1 (MQ=255)
                                                                                                                       atcaACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATc                                >  1:703794/1‑150 (MQ=255)
                                                                                                                       atcaACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATc                                <  1:1113532/150‑1 (MQ=255)
                                                                                                                        tcaACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCg                               <  1:2036748/150‑1 (MQ=255)
                                                                                                                        tcaACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCg                               <  1:551733/150‑1 (MQ=255)
                                                                                                                        tcaACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCg                               <  1:551864/150‑1 (MQ=255)
                                                                                                                           acaaTCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCACTGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGAcc                            >  2:401625/1‑150 (MQ=255)
                                                                                                                               tCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGAt                        >  2:2511608/1‑150 (MQ=255)
                                                                                                                                 aGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATAc                      <  2:3427700/150‑1 (MQ=255)
                                                                                                                                 aGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATAc                      >  2:2994248/1‑150 (MQ=255)
                                                                                                                                  gCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACg                     >  2:876463/1‑150 (MQ=255)
                                                                                                                                    tctGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCCGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGgc                   >  2:2020058/1‑150 (MQ=255)
                                                                                                                                        gATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCAt               >  1:4689315/1‑150 (MQ=255)
                                                                                                                                           ggTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGa            <  2:2530543/150‑1 (MQ=255)
                                                                                                                                            gTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGAc           >  2:1399308/1‑150 (MQ=255)
                                                                                                                                            gTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGAc           <  1:4291715/150‑1 (MQ=255)
                                                                                                                                             tCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACa          <  1:2318954/150‑1 (MQ=255)
                                                                                                                                             tCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACa          <  1:1724185/150‑1 (MQ=255)
                                                                                                                                                tATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACAAAg       >  1:3473152/1‑150 (MQ=255)
                                                                                                                                                 aTTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACAAAGt      >  2:4711873/1‑150 (MQ=255)
                                                                                                                                                 aTTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACAAAGt      <  1:3566953/150‑1 (MQ=255)
                                                                                                                                                 aTTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACAAAGt      <  1:2545444/150‑1 (MQ=255)
                                                                                                                                                    tcttctCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACAAAGTAAg   >  1:718648/1‑150 (MQ=255)
                                                                                                                                                    tcttctCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTATGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACAAAGTAAg   >  1:3971597/1‑150 (MQ=255)
                                                                                                                                                     cttctCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACAAAGTAAGc  >  2:859464/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACAAAGTAAGC  >  CP009273/2950541‑2950839

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: ‑