Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,132,494 T→A 15.7% I255I (ATT→ATA flgI → putative flagellar basal body protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,132,4940TA15.7% 830.5 / 132.8 292I255I (ATT→ATAflgIputative flagellar basal body protein
Reads supporting (aligned to +/- strand):  ref base T (141/105);  new base A (25/21);  total (166/126)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 296 total aligned reads displayed.

CGTGCGTGGATATGGCAGCGCCACCGCGTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTA  >  CP009273/1132347‑1132642
                                                                                                                                                   |                                                                                                                                                    
cgtgcgtgGATATGGCAGCGCCACCGCGTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATaaa                                                                                                                                                    <  1:2963162/150‑1 (MQ=255)
cgtgcgtgGATATGGCAGCGCCACCGCGTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTaa                                                                                                                                                    <  1:2553845/150‑1 (MQ=255)
cgtgcgtgGATATGGCAGCGCCACCGCGTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTaa                                                                                                                                                    >  1:2203974/1‑150 (MQ=255)
    cgtgGATATGGCAGCGCCACCGCGTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTcg                                                                                                                                                >  2:943137/1‑150 (MQ=255)
     gtgGATATGGCAGCGCCACCGCGTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTcgc                                                                                                                                               >  2:1102673/1‑150 (MQ=255)
            tGGCAGCGCCACCGCGTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCgg                                                                                                                                        >  1:856904/1‑150 (MQ=255)
             ggCAGCGCCACCGCGTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGt                                                                                                                                       >  2:3704420/1‑150 (MQ=255)
                aGCGCCACCGCGTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCg                                                                                                                                    >  2:1944061/1‑150 (MQ=255)
                aGCGCCACCGCGTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCg                                                                                                                                    >  2:3058424/1‑150 (MQ=255)
                  cgcCACCGCGTTAGATGCGCTGACGATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCggt                                                                                                                                  <  2:2729535/150‑1 (MQ=255)
                  cgcCACCGCGTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCggt                                                                                                                                  <  2:3409627/150‑1 (MQ=255)
                    ccACCGCGTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCggtgg                                                                                                                                >  1:1396667/1‑150 (MQ=255)
                     cACCGCGTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCggtggt                                                                                                                               >  2:2014391/1‑150 (MQ=255)
                       ccGCGTTAGATGCGCGGACTAGTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGtga                                                                                                                             <  2:28444/150‑1 (MQ=255)
                           gTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGaa                                                                                                                         >  2:3205323/1‑150 (MQ=255)
                           gTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGaa                                                                                                                         >  2:1480209/1‑150 (MQ=255)
                            ttAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAAt                                                                                                                        >  2:4214675/1‑150 (MQ=255)
                            ttAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAAt                                                                                                                        <  1:2887835/150‑1 (MQ=255)
                             tAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATc                                                                                                                       >  2:3086345/1‑150 (MQ=255)
                              aGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATcg                                                                                                                      <  2:4269070/150‑1 (MQ=255)
                                aTGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGACTcgcg                                                                                                                    >  1:1827790/1‑150 (MQ=255)
                                      ggACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGa                                                                                                              >  1:2586648/1‑150 (MQ=255)
                                      ggACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGa                                                                                                              >  1:411913/1‑150 (MQ=255)
                                          tATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCt                                                                                                          >  1:1051729/1‑150 (MQ=255)
                                           aTTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTc                                                                                                         >  1:334581/1‑150 (MQ=255)
                                           aTTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTc                                                                                                         >  2:775000/1‑150 (MQ=255)
                                           aTTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTc                                                                                                         <  1:1643088/150‑1 (MQ=255)
                                                ggTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAg                                                                                                    >  2:4097533/1‑150 (MQ=255)
                                                ggTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAg                                                                                                    <  1:3274263/150‑1 (MQ=255)
                                                   gcgcGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTg                                                                                                 <  2:2472319/150‑1 (MQ=255)
                                                    cgcgTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTgc                                                                                                <  1:3582925/150‑1 (MQ=255)
                                                       gTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGcgcg                                                                                             >  1:3064911/1‑150 (MQ=255)
                                                             aGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAgcg                                                                                       <  2:2551640/150‑1 (MQ=255)
                                                               tGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCa                                                                                     <  2:637858/150‑1 (MQ=255)
                                                               tGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCa                                                                                     <  2:2092196/150‑1 (MQ=255)
                                                                 gCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAgg                                                                                   >  1:633150/1‑150 (MQ=255)
                                                                   aaCAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAgggg                                                                                 >  2:578259/1‑150 (MQ=255)
                                                                     cAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGaa                                                                               <  1:923684/150‑1 (MQ=255)
                                                                     cAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGaa                                                                               >  2:2221603/1‑150 (MQ=255)
                                                                      aGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAAt                                                                              <  2:1645050/150‑1 (MQ=255)
                                                                           ccAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAAtctctc                                                                         >  1:2638272/1‑150 (MQ=255)
                                                                               gTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTa                                                                     >  2:1183549/1‑150 (MQ=255)
                                                                                 ccGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCGAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAAc                                                                   >  1:1727951/1‑150 (MQ=255)
                                                                                  cGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACa                                                                  >  1:2371489/1‑150 (MQ=255)
                                                                                  cGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACa                                                                  <  2:1294030/150‑1 (MQ=255)
                                                                                  cGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACa                                                                  <  2:50974/150‑1 (MQ=255)
                                                                                   gCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAg                                                                 >  1:833091/1‑150 (MQ=255)
                                                                                    cTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGt                                                                >  1:3324338/1‑150 (MQ=255)
                                                                                     gtCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGtt                                                               >  1:973720/2‑150 (MQ=255)
                                                                                       ccTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTaa                                                             >  1:590097/1‑150 (MQ=255)
                                                                                             cGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCa                                                       >  1:1962437/1‑150 (MQ=255)
                                                                                             cGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCa                                                       <  2:211264/150‑1 (MQ=255)
                                                                                              gATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAg                                                      <  1:1164627/150‑1 (MQ=255)
                                                                                               atatTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAgg                                                     <  2:835170/150‑1 (MQ=255)
                                                                                                tatTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGc                                                    >  1:497822/1‑150 (MQ=255)
                                                                                                  ttCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCa                                                  <  2:3159838/150‑1 (MQ=255)
                                                                                                  ttCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCa                                                  <  2:1743832/150‑1 (MQ=255)
                                                                                                   tCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCaa                                                 >  1:1313270/1‑150 (MQ=255)
                                                                                                    cAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAAt                                                >  2:3865165/1‑150 (MQ=255)
                                                                                                     aGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATg                                               >  1:3745376/1‑150 (MQ=255)
                                                                                                      gAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGt                                              >  2:1803461/1‑150 (MQ=255)
                                                                                                      gAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGt                                              >  2:4119733/1‑150 (MQ=255)
                                                                                                        atatGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTca                                            <  1:800054/150‑1 (MQ=255)
                                                                                                          atGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTcagc                                          <  1:1763482/150‑1 (MQ=255)
                                                                                                            gCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTcagcca                                        <  1:2762323/150‑1 (MQ=255)
                                                                                                                    aTGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATa                                >  2:1886332/1‑150 (MQ=255)
                                                                                                                      gTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATaca                              <  1:1480209/150‑1 (MQ=255)
                                                                                                                         aCCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCg                           >  2:2677707/1‑150 (MQ=255)
                                                                                                                         aCCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCg                           >  2:4035665/1‑150 (MQ=255)
                                                                                                                          ccccGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGt                          <  2:734282/150‑1 (MQ=255)
                                                                                                                          ccccGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGt                          <  1:4287683/150‑1 (MQ=255)
                                                                                                                          ccccGCAGGACGCTAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGt                          >  2:2313310/1‑150 (MQ=255)
                                                                                                                           cccGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGtt                         >  2:2666606/1‑150 (MQ=255)
                                                                                                                              gCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTtgg                      <  1:3179489/150‑1 (MQ=255)
                                                                                                                               cAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTtggt                     <  2:2203974/150‑1 (MQ=255)
                                                                                                                               cAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTtggt                     <  2:1500645/150‑1 (MQ=255)
                                                                                                                                aGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTtggtg                    <  2:435833/150‑1 (MQ=255)
                                                                                                                                 ggACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTtggtgg                   <  1:1102673/150‑1 (MQ=255)
                                                                                                                                  gACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGGTtggtggt                  >  1:360822/1‑150 (MQ=255)
                                                                                                                                   aCGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTtggtggtg                 >  2:935050/1‑150 (MQ=255)
                                                                                                                                   aCGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTtggtggtg                 <  2:4159590/150‑1 (MQ=255)
                                                                                                                                      cTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGAc              >  1:884467/1‑150 (MQ=255)
                                                                                                                                      cTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGAc              >  2:356524/1‑150 (MQ=255)
                                                                                                                                      cTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGAc              >  1:3157943/1‑150 (MQ=255)
                                                                                                                                       tAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACa             >  1:4029883/1‑150 (MQ=255)
                                                                                                                                       tAAAGTAGTGATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCCGATACACCGTTTGGTGGTGGACa             >  2:1964558/1‑150 (MQ=255)
                                                                                                                                          agtagtGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGAc          >  2:3670005/1‑150 (MQ=255)
                                                                                                                                          agtagtGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGAc          <  1:3434805/150‑1 (MQ=255)
                                                                                                                                          agtagtGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGAc          <  2:870276/150‑1 (MQ=255)
                                                                                                                                          agtagtGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGAc          <  1:3907816/150‑1 (MQ=255)
                                                                                                                                           gtagtGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCCGATACCCCGTTTGGTGGTGGACAGACt         >  2:2966903/1‑150 (MQ=255)
                                                                                                                                             agtGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACtgt       >  1:1585152/1‑150 (MQ=255)
                                                                                                                                             agtGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACtgt       >  1:4091836/1‑150 (MQ=255)
                                                                                                                                             agtGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACtgt       >  1:3876492/1‑150 (MQ=255)
                                                                                                                                             agtGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACtgt       <  1:1944372/150‑1 (MQ=255)
                                                                                                                                                gATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGt    >  1:720470/1‑150 (MQ=255)
                                                                                                                                                gATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGt    <  1:1695091/150‑1 (MQ=255)
                                                                                                                                                gATAAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGt    >  1:2451568/1‑150 (MQ=255)
                                                                                                                                                 aTTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGtt   >  1:2793219/1‑150 (MQ=255)
                                                                                                                                                  ttAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTa  <  1:1423206/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGTGCGTGGATATGGCAGCGCCACCGCGTTAGATGCGCGGACTATTCAGGTGCGCGTACCGAGTGGCAACAGTTCCCAGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTA  >  CP009273/1132347‑1132642

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: