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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,950,6900CA24.7% 565.3 / 168.6 238K518N (AAG→AATptrAprotease III
Reads supporting (aligned to +/- strand):  ref base C (70/109);  new base A (30/29);  total (100/138)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 251 total aligned reads displayed.

ATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACAAAGTAAG  >  CP009273/2950541‑2950838
                                                                                                                                                     |                                                                                                                                                    
  cACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGtatttatt                                                                                                                                                    >  1:1280329/1‑150 (MQ=255)
   aCCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATttcttc                                                                                                                                                   >  2:4366664/1‑150 (MQ=255)
   aCCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATttcttc                                                                                                                                                   >  1:2540524/1‑150 (MQ=255)
    ccTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTtcttct                                                                                                                                                  >  2:2517730/1‑150 (MQ=255)
     cTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTtctc                                                                                                                                                 >  2:557477/1‑150 (MQ=255)
      tGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTctct                                                                                                                                                >  1:1412155/1‑150 (MQ=255)
      tGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTctct                                                                                                                                                <  1:1307798/150‑1 (MQ=255)
       gATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTg                                                                                                                                               <  1:1349633/150‑1 (MQ=255)
       gATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTg                                                                                                                                               <  1:3683600/150‑1 (MQ=255)
       gATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTg                                                                                                                                               >  1:2046776/1‑150 (MQ=255)
        aTTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGa                                                                                                                                              >  1:3130582/1‑150 (MQ=255)
         ttGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGAc                                                                                                                                             <  1:4308698/150‑1 (MQ=255)
            cGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTa                                                                                                                                          <  2:3474405/150‑1 (MQ=255)
             gggCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTaa                                                                                                                                         >  1:592314/1‑150 (MQ=255)
             gggCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTaa                                                                                                                                         >  1:2619330/1‑150 (MQ=255)
              ggCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAAt                                                                                                                                        <  2:3529234/150‑1 (MQ=255)
              ggCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAAt                                                                                                                                        <  1:558207/150‑1 (MQ=255)
              ggCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAAt                                                                                                                                        <  2:3438332/150‑1 (MQ=255)
               gcgcTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATc                                                                                                                                       <  1:4204265/150‑1 (MQ=255)
                 gcTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAg                                                                                                                                     <  1:526761/150‑1 (MQ=255)
                  cTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGc                                                                                                                                    <  1:4656088/150‑1 (MQ=255)
                    gTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCga                                                                                                                                  <  2:155219/150‑1 (MQ=255)
                      ccATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCgaga                                                                                                                                >  1:4347385/1‑150 (MQ=255)
                      ccATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCgaga                                                                                                                                >  1:1347673/1‑150 (MQ=255)
                       cATGGCTTTCGGATTACGCAAAATCAGGCTGACCTCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGaa                                                                                                                               <  2:1428672/150‑1 (MQ=255)
                       cATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGaa                                                                                                                               <  1:3323624/150‑1 (MQ=255)
                         tGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAat                                                                                                                             <  2:2910710/150‑1 (MQ=255)
                         tGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAat                                                                                                                             <  2:477917/150‑1 (MQ=255)
                              ttCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAAtcatca                                                                                                                        >  1:1897065/1‑150 (MQ=255)
                                cGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAgg                                                                                                                      <  2:1438263/150‑1 (MQ=255)
                                cGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAgg                                                                                                                      <  2:4638976/150‑1 (MQ=255)
                                cGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAgg                                                                                                                      >  2:3786173/1‑150 (MQ=255)
                                 ggATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGa                                                                                                                     >  2:3848209/1‑150 (MQ=255)
                                   aTTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAat                                                                                                                   <  1:2875195/150‑1 (MQ=255)
                                     tACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAatat                                                                                                                 <  1:3561916/150‑1 (MQ=255)
                                       cGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATaa                                                                                                               <  1:3900700/150‑1 (MQ=255)
                                         cAAAATCAGGCTGACATCAGCTTTTGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAgg                                                                                                             <  2:3538316/150‑1 (MQ=255)
                                           aaaTCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGACTATAAGGGt                                                                                                           >  2:4215107/1‑150 (MQ=255)
                                           aaaTCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGt                                                                                                           <  1:773204/150‑1 (MQ=255)
                                            aaTCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGtt                                                                                                          <  1:2914215/150‑1 (MQ=255)
                                             aTCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTa                                                                                                         <  2:2276523/150‑1 (MQ=255)
                                             aTCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTa                                                                                                         >  1:4457823/1‑150 (MQ=255)
                                              tCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTaa                                                                                                        >  2:1950355/1‑150 (MQ=255)
                                                   cTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGctct                                                                                                   <  2:1223293/150‑1 (MQ=255)
                                                   cTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGctct                                                                                                   >  1:3877963/1‑150 (MQ=255)
                                                      aCATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGc                                                                                                >  2:1831685/1‑150 (MQ=255)
                                                      aCATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGc                                                                                                >  2:4672824/1‑150 (MQ=255)
                                                      aCATCAGCTTTGGGCTCGCTGGAAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGc                                                                                                <  2:3002173/150‑1 (MQ=255)
                                                        aTCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCaa                                                                                              >  1:998174/1‑150 (MQ=255)
                                                             cTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAagaga                                                                                         >  2:4048359/1‑150 (MQ=255)
                                                                       gCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTc                                                                               <  1:3008590/150‑1 (MQ=255)
                                                                        cTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTcg                                                                              <  2:1783016/150‑1 (MQ=255)
                                                                           gCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTcggcg                                                                           <  2:2739055/150‑1 (MQ=255)
                                                                             aaaaTAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCggcggc                                                                         >  1:1754909/1‑150 (MQ=255)
                                                                               aaTAACGGCTTGGCGCATACCCCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCtt                                                                       <  1:2431975/150‑1 (MQ=255)
                                                                               aaTAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCtt                                                                       <  2:1162549/150‑1 (MQ=255)
                                                                                 tAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCtttt                                                                     <  1:1392132/150‑1 (MQ=255)
                                                                                 tAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCtttt                                                                     <  1:3232026/150‑1 (MQ=255)
                                                                                 tAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCtttt                                                                     <  2:1144331/150‑1 (MQ=255)
                                                                                  aaCGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCttttt                                                                    <  1:704858/150‑1 (MQ=255)
                                                                                   aCGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCtttttt                                                                   <  2:457234/150‑1 (MQ=255)
                                                                                   aCGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCtttttt                                                                   <  1:199964/150‑1 (MQ=255)
                                                                                   aCGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCtttttt                                                                   <  1:4679136/150‑1 (MQ=255)
                                                                                   aCGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCtttttt                                                                   <  1:2225947/150‑1 (MQ=255)
                                                                                     ggCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCt                                                                 >  1:4099611/1‑150 (MQ=255)
                                                                                       cTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGc                                                               >  2:1534349/1‑150 (MQ=255)
                                                                                               aTACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGc                                                       >  2:2142179/1‑150 (MQ=255)
                                                                                               aTACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGc                                                       >  1:1881027/1‑150 (MQ=255)
                                                                                                tACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCCGTCGGCg                                                      >  2:3459745/1‑150 (MQ=255)
                                                                                                 acacCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGa                                                     >  2:175312/1‑150 (MQ=255)
                                                                                                  caccacGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGaa                                                    >  2:2334815/1‑150 (MQ=255)
                                                                                                     cacGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGt                                                 <  2:1582715/150‑1 (MQ=255)
                                                                                                     cacGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGt                                                 >  1:2099160/1‑150 (MQ=255)
                                                                                                       cgcgCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGttt                                               >  2:192246/1‑150 (MQ=255)
                                                                                                       cgcgCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGttt                                               <  2:4099611/150‑1 (MQ=255)
                                                                                                        gcgcAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTtg                                              <  1:1804429/150‑1 (MQ=255)
                                                                                                           cAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTATGCCAGTCGGCGAAAGTTTGTgc                                           >  1:148790/1‑150 (MQ=255)
                                                                                                            aGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTgcg                                          >  2:897202/1‑150 (MQ=255)
                                                                                                              aTTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCt                                        <  2:722290/150‑1 (MQ=255)
                                                                                                              aTTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCt                                        <  2:1272298/150‑1 (MQ=255)
                                                                                                                 cGACTCATCAACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGAt                                     <  1:3834641/150‑1 (MQ=255)
                                                                                                                  gACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGAtt                                    <  1:4163750/150‑1 (MQ=255)
                                                                                                                     tcatcaACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTa                                 <  1:1056447/150‑1 (MQ=255)
                                                                                                                     tcatcaACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTa                                 <  1:1720082/150‑1 (MQ=255)
                                                                                                                     tcatcaACAATCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTa                                 <  1:3508685/150‑1 (MQ=255)
                                                                                                                       atcaACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATc                               <  1:78008/150‑1 (MQ=255)
                                                                                                                          aacaaTCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGAc                            >  1:2481339/1‑150 (MQ=255)
                                                                                                                               tCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAAACAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGAt                       <  2:2044257/150‑1 (MQ=255)
                                                                                                                               tCAGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGAt                       <  1:4036761/150‑1 (MQ=255)
                                                                                                                                cAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATa                      >  2:3332459/1‑150 (MQ=255)
                                                                                                                                 aGCTCTGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATAc                     <  1:3980685/150‑1 (MQ=255)
                                                                                                                                  gCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCCGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACg                    <  1:3273095/150‑1 (MQ=255)
                                                                                                                                    tctGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGgc                  <  1:632746/150‑1 (MQ=255)
                                                                                                                                      tGGATGGTCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGgcgc                >  2:3094169/1‑150 (MQ=255)
                                                                                                                                          tGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCg            <  2:2888128/150‑1 (MQ=255)
                                                                                                                                             tCGTATTTATTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACa         <  1:1173916/150‑1 (MQ=255)
                                                                                                                                              cGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACaa        >  1:3002901/1‑150 (MQ=255)
                                                                                                                                                  tttCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACAAAGTa    >  2:1009789/1‑150 (MQ=255)
                                                                                                                                                  tttCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACAAAGTa    >  2:149705/1‑150 (MQ=255)
                                                                                                                                                    tcttctCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACAAAGTAAg  <  2:1754909/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAGTTTGTGCGCTGATTTTATCGACCTGATACGGCGCATCGACAAAGTAAG  >  CP009273/2950541‑2950838

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: