Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,950,431 T→A 9.5% N605Y (AAT→TAT)  ptrA ← protease III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,950,4310TA 9.5% 763.7 / 56.9 253N605Y (AAT→TAT) ptrAprotease III
Reads supporting (aligned to +/- strand):  ref base T (117/112);  new base A (11/13);  total (128/125)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 260 total aligned reads displayed.

AACGCTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATT  >  CP009273/2950283‑2950578
                                                                                                                                                    |                                                                                                                                                   
 aCGCTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTa                                                                                                                                                    >  2:364653/1‑149 (MQ=255)
 aCGCTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTa                                                                                                                                                    <  1:2662986/149‑1 (MQ=255)
 aCGCTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTa                                                                                                                                                    <  1:598187/149‑1 (MQ=255)
  cGCTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTaa                                                                                                                                                   <  1:3429478/149‑1 (MQ=255)
  cGCTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTaa                                                                                                                                                   <  1:559553/149‑1 (MQ=255)
     tttACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCa                                                                                                                                                <  1:2874575/149‑1 (MQ=255)
      ttACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCAt                                                                                                                                               >  1:3236598/1‑149 (MQ=255)
         cccTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAg                                                                                                                                            >  1:2145295/1‑149 (MQ=255)
           cTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGc                                                                                                                                          >  2:1673699/1‑149 (MQ=255)
           cTTTTCTGCGGAATACATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGc                                                                                                                                          <  1:2845555/149‑1 (MQ=255)
             tttCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCg                                                                                                                                        >  2:1910978/1‑149 (MQ=255)
              ttCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCgt                                                                                                                                       <  2:3610611/149‑1 (MQ=255)
               tCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCgtt                                                                                                                                      <  1:2231125/149‑1 (MQ=255)
                 tGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCgttgt                                                                                                                                    <  2:1899021/149‑1 (MQ=255)
                 tGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCgttgt                                                                                                                                    >  2:3159269/1‑149 (MQ=255)
                 tGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGCAACCATTAGCATTAACCATAAGGCCgttgt                                                                                                                                    >  2:3298717/1‑149 (MQ=255)
                   cGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTa                                                                                                                                  >  1:2079380/1‑149 (MQ=255)
                    ggAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAg                                                                                                                                 >  1:1734231/1‑149 (MQ=255)
                    ggAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAg                                                                                                                                 >  2:1585574/1‑149 (MQ=255)
                       aTCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGt                                                                                                                              >  1:1293017/1‑149 (MQ=255)
                          catcatCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTtgg                                                                                                                           <  1:131806/149‑1 (MQ=255)
                           atcatcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTATCGTtggt                                                                                                                          >  1:1966242/1‑149 (MQ=255)
                           atcatcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTtggt                                                                                                                          <  2:3047158/149‑1 (MQ=255)
                             catcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTtggtgg                                                                                                                        <  1:1519343/149‑1 (MQ=255)
                               tcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGaa                                                                                                                      >  1:1113782/1‑149 (MQ=255)
                               tcTGGTTATACCAGGACTTCGCCTGCTCAAGATGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGATGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGaa                                                                                                                      >  1:2977835/1‑149 (MQ=255)
                                  ggTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGaaaaa                                                                                                                   >  2:2467553/1‑149 (MQ=255)
                                   gTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAAc                                                                                                                  <  2:210687/149‑1 (MQ=255)
                                   gTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCCAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAAc                                                                                                                  >  1:3484957/1‑149 (MQ=255)
                                      ataCCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTa                                                                                                               <  2:2274908/149‑1 (MQ=255)
                                      ataCCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTa                                                                                                               <  2:3539164/149‑1 (MQ=255)
                                             gACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGccacc                                                                                                        <  2:4101089/149‑1 (MQ=255)
                                                 tCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAAcc                                                                                                    <  1:1215160/149‑1 (MQ=255)
                                                   gCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGa                                                                                                  <  2:2862518/149‑1 (MQ=255)
                                                   gCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGa                                                                                                  <  2:3005716/149‑1 (MQ=255)
                                                    ccTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGAc                                                                                                 >  1:1454280/1‑149 (MQ=255)
                                                       gCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGcc                                                                                              <  1:3101098/149‑1 (MQ=255)
                                                         tCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTg                                                                                            <  1:1005856/149‑1 (MQ=255)
                                                           aaGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGt                                                                                          >  2:3654474/1‑149 (MQ=255)
                                                            aGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGtt                                                                                         >  2:2356641/1‑149 (MQ=255)
                                                                gATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCtt                                                                                     <  1:1836495/149‑1 (MQ=255)
                                                                 aTCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTa                                                                                    >  2:3357726/1‑149 (MQ=255)
                                                                     tCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTg                                                                                <  1:2598317/149‑1 (MQ=255)
                                                                          aGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCaa                                                                           >  2:1888777/1‑149 (MQ=255)
                                                                          aGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCaa                                                                           >  2:4023522/1‑149 (MQ=255)
                                                                             ggTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAgcg                                                                        >  1:2942129/1‑149 (MQ=255)
                                                                              gTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAgcgc                                                                       <  1:468632/149‑1 (MQ=255)
                                                                               tataGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGcc                                                                      >  1:3328637/1‑149 (MQ=255)
                                                                                 taGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAg                                                                    >  1:511877/1‑149 (MQ=255)
                                                                                       aaGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGc                                                              <  1:2290394/149‑1 (MQ=255)
                                                                                         gTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCga                                                            <  2:1296395/149‑1 (MQ=255)
                                                                                          tACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCgag                                                           >  2:2516950/1‑149 (MQ=255)
                                                                                           aCCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGGTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCGGGTTGCTTAACTGATCAAGCGCCAGCCCTGCgaga                                                          >  2:1608643/1‑149 (MQ=255)
                                                                                            ccccTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGAt                                                         <  1:3862095/149‑1 (MQ=255)
                                                                                             cccTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATa                                                        <  2:1478807/149‑1 (MQ=255)
                                                                                              ccTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATaa                                                       <  1:2478493/149‑1 (MQ=255)
                                                                                              ccTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATaa                                                       <  1:407553/149‑1 (MQ=255)
                                                                                                  gAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCAt                                                   <  2:1826052/149‑1 (MQ=255)
                                                                                                  gAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAACCAt                                                   <  2:3865151/149‑1 (MQ=255)
                                                                                                   aGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCAtt                                                  >  2:1569682/1‑149 (MQ=255)
                                                                                                      aaTGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATACGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAg                                               <  1:364653/149‑1 (MQ=255)
                                                                                                      aaTGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAg                                               >  2:3775043/1‑149 (MQ=255)
                                                                                                      aaTGCCTGGAACAGCTGCGGCAGACGCTGGGGGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAg                                               >  1:4042470/1‑149 (MQ=255)
                                                                                                       aTGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAgc                                              >  1:3127509/1‑149 (MQ=255)
                                                                                                       aTGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAgc                                              >  2:3984013/1‑149 (MQ=255)
                                                                                                            tGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCaaa                                         >  2:114552/1‑149 (MQ=255)
                                                                                                            tGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCaaa                                         >  1:2629296/1‑149 (MQ=255)
                                                                                                              gAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACa                                       <  1:3400607/149‑1 (MQ=255)
                                                                                                              gAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACa                                       <  1:763488/149‑1 (MQ=255)
                                                                                                               aaCAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATCATCATTGAGCGCAAACAt                                      >  2:2578432/1‑149 (MQ=255)
                                                                                                                 cAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCa                                    >  1:82262/1‑149 (MQ=255)
                                                                                                                    cTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCt                                 <  2:901309/149‑1 (MQ=255)
                                                                                                                       cGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGAt                              <  1:2604863/149‑1 (MQ=255)
                                                                                                                         gCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTg                            >  2:425519/1‑149 (MQ=255)
                                                                                                                         gCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTg                            >  2:3719564/1‑149 (MQ=255)
                                                                                                                           aGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCg                          >  1:2435598/1‑149 (MQ=255)
                                                                                                                            gACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCgg                         >  2:1375741/1‑149 (MQ=255)
                                                                                                                            gACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCgg                         <  1:3983492/149‑1 (MQ=255)
                                                                                                                            gACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCgg                         >  1:3748718/1‑149 (MQ=255)
                                                                                                                              cGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGgc                       <  2:294122/149‑1 (MQ=255)
                                                                                                                              cGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGgc                       <  2:746806/149‑1 (MQ=255)
                                                                                                                                  gggTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTg                   >  2:2748109/1‑149 (MQ=255)
                                                                                                                                  gggTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTg                   <  2:2758067/149‑1 (MQ=255)
                                                                                                                                      gtAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGTAAAAACTTATGCCACCAACCGACGCCTGGTTTCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCa               <  2:1966242/149‑1 (MQ=255)
                                                                                                                                       tAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCAt              >  2:2659242/1‑149 (MQ=255)
                                                                                                                                        aaCCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATg             >  1:860898/1‑149 (MQ=255)
                                                                                                                                        aaCCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATg             <  1:3016127/149‑1 (MQ=255)
                                                                                                                                         aCCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATgg            >  1:2226101/1‑149 (MQ=255)
                                                                                                                                         aCCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATgg            >  2:3966638/1‑149 (MQ=255)
                                                                                                                                         aCCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATgg            <  1:3626187/149‑1 (MQ=255)
                                                                                                                                          ccATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGc           <  2:3340007/149‑1 (MQ=255)
                                                                                                                                            aTTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCtt         <  2:2070307/149‑1 (MQ=255)
                                                                                                                                               aGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCg      >  1:1766033/1‑149 (MQ=255)
                                                                                                                                               aGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCg      >  2:2039022/1‑149 (MQ=255)
                                                                                                                                               aGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCg      <  1:789402/149‑1 (MQ=255)
                                                                                                                                                  aTTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGAt   <  1:2352767/149‑1 (MQ=255)
                                                                                                                                                   ttAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGAtt  <  1:1634659/149‑1 (MQ=255)
                                                                                                                                                   ttAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGAtt  <  1:2723085/149‑1 (MQ=255)
                                                                                                                                                   ttAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGAtt  <  2:1734231/149‑1 (MQ=255)
                                                                                                                                                   ttAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGAtt  >  2:92697/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACGCTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATT  >  CP009273/2950283‑2950578

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: