Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,125,785 C→A 5.4% A179S (GCC→TCC)  flgA ← assembly protein for flagellar basal‑body periplasmic P ring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,125,7850CA 5.4% 1317.2 / 28.5 334A179S (GCC→TCC) flgAassembly protein for flagellar basal‑body periplasmic P ring
Reads supporting (aligned to +/- strand):  ref base C (162/154);  new base A (11/7);  total (173/161)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 348 total aligned reads displayed.

AGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCAC  >  CP009273/1125637‑1125932
                                                                                                                                                    |                                                                                                                                                   
aGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACTCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGgc                                                                                                                                                     >  2:2059971/1‑149 (MQ=255)
aGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGgc                                                                                                                                                     >  1:2521228/1‑149 (MQ=255)
aGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGgc                                                                                                                                                     >  2:2301386/1‑149 (MQ=255)
 gggACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGgcg                                                                                                                                                    >  2:1220934/1‑149 (MQ=255)
   gACGGGTAATCTTTAACAGCTTCCAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCt                                                                                                                                                  >  2:1756089/1‑149 (MQ=255)
   gACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCACCAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCt                                                                                                                                                  <  2:1659593/149‑1 (MQ=255)
   gACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCt                                                                                                                                                  >  2:1418217/1‑149 (MQ=255)
      gggTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTaaa                                                                                                                                               <  1:2031730/149‑1 (MQ=255)
          aaTCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCa                                                                                                                                           >  1:3674600/1‑149 (MQ=255)
          aaTCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCa                                                                                                                                           <  2:3397515/149‑1 (MQ=255)
                 aaCAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGc                                                                                                                                    >  1:1452653/1‑149 (MQ=255)
                 aaCAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGc                                                                                                                                    >  2:4650131/1‑149 (MQ=255)
                 aaCAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGc                                                                                                                                    >  1:877945/1‑149 (MQ=255)
                  aCAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGAGCTAAACCCATCACCGCt                                                                                                                                   >  1:3175942/1‑149 (MQ=255)
                       ttACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGa                                                                                                                              >  2:725469/1‑149 (MQ=255)
                        tACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGAt                                                                                                                             >  1:4498620/1‑149 (MQ=255)
                             gTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACAt                                                                                                                        >  2:62624/1‑149 (MQ=255)
                             gTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACAt                                                                                                                        <  1:1973190/149‑1 (MQ=255)
                                tataAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGa                                                                                                                     >  1:2615505/1‑149 (MQ=255)
                                 ataAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGAc                                                                                                                    >  1:4071138/1‑149 (MQ=255)
                                 ataAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGAc                                                                                                                    >  1:4071481/1‑149 (MQ=255)
                                   aaGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGAcgc                                                                                                                  >  1:1064527/1‑149 (MQ=255)
                                         atTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTc                                                                                                            >  1:1639235/1‑149 (MQ=255)
                                          ttCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTcc                                                                                                           >  1:3476773/1‑149 (MQ=255)
                                          ttCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGAACACATTGACGCGTTGTcc                                                                                                           <  2:245777/149‑1 (MQ=255)
                                          ttCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACACATCACCGCTGGCGATCACATTGACGCGTTGTcc                                                                                                           >  1:2012924/1‑149 (MQ=255)
                                           tCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTccc                                                                                                          >  1:307602/1‑149 (MQ=255)
                                            cccATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCg                                                                                                         >  1:1402623/1‑149 (MQ=255)
                                             ccATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGc                                                                                                        >  2:1455238/1‑149 (MQ=255)
                                              cATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCt                                                                                                       >  2:800730/1‑149 (MQ=255)
                                                tCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCttt                                                                                                     <  1:50523/149‑1 (MQ=255)
                                                  tGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTa                                                                                                   >  1:2203472/1‑149 (MQ=255)
                                                    cATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTAcc                                                                                                 <  2:3778061/149‑1 (MQ=255)
                                                     aTCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTAccc                                                                                                <  1:1130255/149‑1 (MQ=255)
                                                      tCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCg                                                                                               >  1:2064283/1‑149 (MQ=255)
                                                         acaacGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCa                                                                                            <  1:731308/149‑1 (MQ=255)
                                                         acaacGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCa                                                                                            >  2:2384041/1‑149 (MQ=255)
                                                         acaacGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCa                                                                                            <  1:4230644/149‑1 (MQ=255)
                                                         acaacGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCACATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCa                                                                                            >  1:4393300/1‑149 (MQ=255)
                                                          caacGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCAt                                                                                           <  2:4344941/149‑1 (MQ=255)
                                                           aacGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATg                                                                                          <  1:2271794/149‑1 (MQ=255)
                                                           aacGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATg                                                                                          >  1:3620139/1‑149 (MQ=255)
                                                             cgccgcTGTCTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCTCCATGcc                                                                                        <  2:1279594/149‑1 (MQ=255)
                                                             cgccgcTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGAGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGcc                                                                                        <  1:3005798/149‑1 (MQ=255)
                                                             cgccgcTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGAGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGcc                                                                                        <  1:4544587/149‑1 (MQ=255)
                                                              gccgcTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCt                                                                                       >  2:3589291/1‑149 (MQ=255)
                                                              gccgcTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCt                                                                                       <  1:4632099/149‑1 (MQ=255)
                                                                    gACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGAGCTAAACCCATCACCGCTGGCGATCCCATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGaa                                                                                 >  2:3652520/1‑149 (MQ=255)
                                                                      cTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAAc                                                                               <  1:3307868/149‑1 (MQ=255)
                                                                          cTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTggg                                                                           >  1:556120/1‑149 (MQ=255)
                                                                           tGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGt                                                                          <  1:2483963/149‑1 (MQ=255)
                                                                               ccGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAAc                                                                      <  2:4444138/149‑1 (MQ=255)
                                                                               ccGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGGTGTCCAGCTTTTACCCGCCATGCCTGGCGAAACTGGGGTAAc                                                                      >  2:828998/1‑149 (MQ=255)
                                                                                 gATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTg                                                                    >  2:217207/1‑149 (MQ=255)
                                                                                 gATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTg                                                                    >  1:2165692/1‑149 (MQ=255)
                                                                                      cATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAg                                                               <  2:3208690/149‑1 (MQ=255)
                                                                                      cATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAg                                                               >  2:4582763/1‑149 (MQ=255)
                                                                                      cATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAg                                                               >  1:3894212/1‑149 (MQ=255)
                                                                                      cATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAg                                                               >  2:1177168/1‑149 (MQ=255)
                                                                                       aTGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAgg                                                              <  1:3446012/149‑1 (MQ=255)
                                                                                        tGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGt                                                             <  2:2924936/149‑1 (MQ=255)
                                                                                        tGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGt                                                             <  1:596003/149‑1 (MQ=255)
                                                                                         gcgcACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGtt                                                            <  1:1223428/149‑1 (MQ=255)
                                                                                          cgcACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTg                                                           <  1:1002259/149‑1 (MQ=255)
                                                                                            cACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTCACTGGATAGGTTGAt                                                         >  2:785304/1‑149 (MQ=255)
                                                                                             aCCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATc                                                        >  1:3359616/1‑149 (MQ=255)
                                                                                              cccGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGAGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCg                                                       >  2:1671352/1‑149 (MQ=255)
                                                                                                cgcgCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCggg                                                     >  2:4073808/1‑149 (MQ=255)
                                                                                                    cATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGAt                                                 <  1:2132974/149‑1 (MQ=255)
                                                                                                      ttCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAg                                               >  2:2029812/1‑149 (MQ=255)
                                                                                                       tCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGa                                              <  1:4251918/149‑1 (MQ=255)
                                                                                                        cTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGAt                                             <  2:3897273/149‑1 (MQ=255)
                                                                                                           tgCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATcgc                                          <  2:2810298/149‑1 (MQ=255)
                                                                                                             cGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCgcgc                                        <  2:3507743/149‑1 (MQ=255)
                                                                                                               aCGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAg                                      >  1:4687155/1‑149 (MQ=255)
                                                                                                                  gCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCt                                   >  1:1496285/1‑149 (MQ=255)
                                                                                                                     gCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAAt                                >  1:1399518/1‑149 (MQ=255)
                                                                                                                      cATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATg                               <  2:4071358/149‑1 (MQ=255)
                                                                                                                       aTTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATgg                              >  1:3107303/1‑149 (MQ=255)
                                                                                                                       aTTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATgg                              <  2:1343038/149‑1 (MQ=255)
                                                                                                                        ttGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGc                             <  2:4160596/149‑1 (MQ=255)
                                                                                                                           ttCAGCGCCTGACCGTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACAATGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTTGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATc                          <  1:231103/149‑1 (MQ=255)
                                                                                                                             cAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATcaa                        <  2:877945/149‑1 (MQ=255)
                                                                                                                             cAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATcaa                        >  1:1585671/1‑149 (MQ=255)
                                                                                                                             cAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATcaa                        <  2:3455448/149‑1 (MQ=255)
                                                                                                                                cgcCTGACCTTCTGCGTTGGAGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATcaacaa                     <  1:1285996/149‑1 (MQ=255)
                                                                                                                                  ccTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGt                   <  1:2424970/149‑1 (MQ=255)
                                                                                                                                    tGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttg                 >  1:3406857/1‑149 (MQ=255)
                                                                                                                                    tGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttg                 >  1:560973/1‑149 (MQ=255)
                                                                                                                                     gACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttga                <  1:3844232/149‑1 (MQ=255)
                                                                                                                                      aCCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttgat               <  2:1890911/149‑1 (MQ=255)
                                                                                                                                       ccTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttgatt              <  2:371377/149‑1 (MQ=255)
                                                                                                                                       ccTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttgatt              >  2:770970/1‑149 (MQ=255)
                                                                                                                                        cTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttgattg             >  1:3204388/1‑149 (MQ=255)
                                                                                                                                         ttCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttgattga            <  1:952688/149‑1 (MQ=255)
                                                                                                                                            tGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGatat         >  1:2087287/1‑149 (MQ=255)
                                                                                                                                               gTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCa      <  2:3203520/149‑1 (MQ=255)
                                                                                                                                                ttGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAg     >  1:3040169/1‑149 (MQ=255)
                                                                                                                                                  ggCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCa   <  2:1064527/149‑1 (MQ=255)
                                                                                                                                                   gcgcTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCAc  >  2:2610693/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCAC  >  CP009273/1125637‑1125932

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: