Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,984,987 C→A 100% T38T (ACG→ACT yecA ← UPF0149 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,984,9870CA99.4% 1048.9 / ‑1.8 307T38T (ACG→ACTyecAUPF0149 family protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (141/163);  total (143/164)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 308 total aligned reads displayed.

GCGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGC  >  CP009273/1984841‑1985132
                                                                                                                                                  |                                                                                                                                                 
gCGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGt                                                                                                                                                 <  2:2758109/149‑1 (MQ=255)
  ggTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCa                                                                                                                                               >  1:52557/1‑149 (MQ=255)
   gTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCaa                                                                                                                                              >  2:3863466/1‑149 (MQ=255)
   gTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCaa                                                                                                                                              >  1:3446169/1‑149 (MQ=255)
     aTCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATa                                                                                                                                            <  1:2666497/149‑1 (MQ=255)
       cGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATaaa                                                                                                                                          <  1:1558645/149‑1 (MQ=255)
       cGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATaaa                                                                                                                                          <  1:3947929/149‑1 (MQ=255)
         gCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAAcc                                                                                                                                        <  2:885971/149‑1 (MQ=255)
         gCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAAcc                                                                                                                                        >  1:2378487/1‑149 (MQ=255)
          ccATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCg                                                                                                                                       >  2:4275334/1‑149 (MQ=255)
          ccATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCg                                                                                                                                       <  2:3650935/149‑1 (MQ=255)
           cATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGt                                                                                                                                      >  1:2033349/1‑149 (MQ=255)
            atatGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTc                                                                                                                                     <  2:3190174/149‑1 (MQ=255)
             tatGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTcc                                                                                                                                    <  1:1044296/149‑1 (MQ=255)
              atGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCa                                                                                                                                   <  1:1850327/149‑1 (MQ=255)
               tGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAg                                                                                                                                  >  1:3260892/1‑149 (MQ=255)
               tGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAg                                                                                                                                  >  1:3030183/1‑149 (MQ=255)
                 ttGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCt                                                                                                                                >  2:1285103/1‑149 (MQ=255)
                   gAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTcc                                                                                                                              <  1:1928698/149‑1 (MQ=255)
                     aaaGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGc                                                                                                                            >  2:661682/1‑149 (MQ=255)
                     aaaGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGc                                                                                                                            >  1:760924/1‑149 (MQ=255)
                      aaGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGcc                                                                                                                           >  2:1914584/1‑149 (MQ=255)
                       aGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCca                                                                                                                          >  1:2336527/1‑149 (MQ=255)
                       aGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCca                                                                                                                          <  1:39861/149‑1 (MQ=255)
                          cAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATACACCGTCCAGCTCCGCCACAt                                                                                                                       >  1:3095219/1‑149 (MQ=255)
                              ttCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACTCGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCATCACCGCAGTCAATAAACCGTCCAGCTTCGCCACATCAAg                                                                                                                   <  1:4158665/149‑1 (MQ=255)
                              ttCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAg                                                                                                                   <  2:4179447/149‑1 (MQ=255)
                              ttCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAg                                                                                                                   <  1:1258228/149‑1 (MQ=255)
                                cATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGa                                                                                                                 <  1:2616755/149‑1 (MQ=255)
                                 aTAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGAt                                                                                                                <  1:2928742/149‑1 (MQ=255)
                                 aTAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGAt                                                                                                                <  2:4592216/149‑1 (MQ=255)
                                 aTAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGAt                                                                                                                >  2:4062861/1‑149 (MQ=255)
                                 aTAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGAt                                                                                                                <  1:1725202/149‑1 (MQ=255)
                                 aTAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGAt                                                                                                                <  2:3419398/149‑1 (MQ=255)
                                  tAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATg                                                                                                               >  1:3718566/1‑149 (MQ=255)
                                  tAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATg                                                                                                               <  1:4427236/149‑1 (MQ=255)
                                    aaGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGc                                                                                                             >  2:590378/1‑149 (MQ=255)
                                     aGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCg                                                                                                            >  1:3572503/1‑149 (MQ=255)
                                      gcgcGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGt                                                                                                           <  1:2510847/149‑1 (MQ=255)
                                          gTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTc                                                                                                       >  2:3172989/1‑149 (MQ=255)
                                             aTCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgt                                                                                                    >  1:1031918/1‑149 (MQ=255)
                                             aTCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgt                                                                                                    <  1:910476/149‑1 (MQ=255)
                                             aTCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgt                                                                                                    <  1:4314267/149‑1 (MQ=255)
                                             aTCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgt                                                                                                    >  1:378391/1‑149 (MQ=255)
                                              tctcTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgtg                                                                                                   >  1:2953995/1‑149 (MQ=255)
                                              tctcTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgtg                                                                                                   >  2:1038721/1‑149 (MQ=255)
                                              tctcTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgtg                                                                                                   >  1:561760/1‑149 (MQ=255)
                                                tctTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgtt                                                                                                 >  2:2338621/1‑149 (MQ=255)
                                                  tttCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgttgt                                                                                               <  1:1159002/149‑1 (MQ=255)
                                                   ttCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTa                                                                                              >  1:4459513/1‑149 (MQ=255)
                                                     ctctGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTagt                                                                                            >  2:1502460/1‑147 (MQ=255)
                                                     ctctGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTAtt                                                                                            >  2:3921512/1‑149 (MQ=255)
                                                      tctGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTAttt                                                                                           >  2:3979841/1‑149 (MQ=255)
                                                       ctGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTg                                                                                          <  2:2320359/149‑1 (MQ=255)
                                                        tGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTTCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTgg                                                                                         <  2:3221512/149‑1 (MQ=255)
                                                        tGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTgg                                                                                         <  1:1450225/149‑1 (MQ=255)
                                                        tGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTgg                                                                                         >  2:2643094/1‑149 (MQ=255)
                                                          aCGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTc                                                                                       <  1:4519280/149‑1 (MQ=255)
                                                               cAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAat                                                                                  >  2:4697799/1‑149 (MQ=255)
                                                                     ggCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTc                                                                            <  2:2894694/149‑1 (MQ=255)
                                                                     ggCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTc                                                                            <  2:2894009/149‑1 (MQ=255)
                                                                     ggCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTAGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTc                                                                            >  2:3461929/1‑149 (MQ=255)
                                                                      gCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTcc                                                                           <  2:123060/149‑1 (MQ=255)
                                                                      gCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTcc                                                                           >  2:1945425/1‑149 (MQ=255)
                                                                         caTAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGc                                                                        >  1:2113376/1‑149 (MQ=255)
                                                                         caTAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCGTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGc                                                                        <  2:2463708/149‑1 (MQ=255)
                                                                            aGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTTCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCAt                                                                     <  1:4328904/149‑1 (MQ=255)
                                                                             gTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCAtt                                                                    <  2:1892030/149‑1 (MQ=255)
                                                                             gTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCAtt                                                                    <  2:775761/149‑1 (MQ=255)
                                                                                  ccccaccccaCACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTTCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAAc                                                               <  1:2191224/149‑1 (MQ=255)
                                                                                      accccacACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCAc                                                           >  2:4496412/1‑149 (MQ=255)
                                                                                       ccccacACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACt                                                          <  1:959613/149‑1 (MQ=255)
                                                                                            cacGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGt                                                     >  1:1727208/1‑149 (MQ=255)
                                                                                              cGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGttt                                                   <  2:2915508/149‑1 (MQ=255)
                                                                                              cGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGttt                                                   <  1:1340729/149‑1 (MQ=255)
                                                                                                gCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTaa                                                 <  1:99333/149‑1 (MQ=255)
                                                                                                 cAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAAc                                                <  1:3392847/149‑1 (MQ=255)
                                                                                                 cAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAAc                                                <  1:3646047/149‑1 (MQ=255)
                                                                                                   aCCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACgg                                              >  1:4673897/1‑149 (MQ=255)
                                                                                                     cAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTc                                            <  1:428630/149‑1 (MQ=255)
                                                                                                      aGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTTCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTcc                                           >  2:2378656/1‑149 (MQ=255)
                                                                                                             gTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTc                                    <  1:616838/149‑1 (MQ=255)
                                                                                                                ccGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATa                                 <  2:3269379/149‑1 (MQ=255)
                                                                                                                 cGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATACACCGTCCAGCGTCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATAc                                <  2:2142554/149‑1 (MQ=255)
                                                                                                                 cGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATAc                                <  1:2469068/149‑1 (MQ=255)
                                                                                                                  ggTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACg                               >  2:267542/1‑149 (MQ=255)
                                                                                                                   gTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGc                              >  2:1766522/1‑149 (MQ=255)
                                                                                                                    ttCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCt                             <  1:1610624/149‑1 (MQ=255)
                                                                                                                     tCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCtt                            >  2:3561817/1‑149 (MQ=255)
                                                                                                                     tCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCtt                            <  1:946628/149‑1 (MQ=255)
                                                                                                                        aTCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTc                         <  2:505164/149‑1 (MQ=255)
                                                                                                                           tctTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCtt                      <  1:2746443/149‑1 (MQ=255)
                                                                                                                                 ggAGAACTCAACACCGCAGTCAATAAACCGTCCAGCGTCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAt                <  2:4305801/149‑1 (MQ=255)
                                                                                                                                   agaACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtct              <  1:200625/149‑1 (MQ=255)
                                                                                                                                     aaCTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCa            >  2:4420950/1‑149 (MQ=255)
                                                                                                                                       cTCAACACCGCAGTCACTAACCCGTCCAGCTTCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAAc          <  2:3524542/149‑1 (MQ=255)
                                                                                                                                       cTCAACACCGCAGTCAATAAACCGTCCAGCTTCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAAc          <  2:1057952/149‑1 (MQ=255)
                                                                                                                                       cTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAAc          <  1:4438255/149‑1 (MQ=255)
                                                                                                                                        tCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACt         >  1:2733597/1‑149 (MQ=255)
                                                                                                                                           acacCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGc      >  1:1136741/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGC  >  CP009273/1984841‑1985132

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: