Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 274,735 A→T 93.1% pseudogene (138/348 nt) afuB ← pseudogene, CP4‑6 prophage; ABC transporter permease family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273274,7350AT93.1% 928.0 / 42.3 287pseudogene (138/348 nt)afuBpseudogene, CP4‑6 prophage; ABC transporter permease family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (13/6);  new base T (153/114);  total (167/120)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 287 total aligned reads displayed.

GAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATG  >  CP009273/274588‑274880
                                                                                                                                                   |                                                                                                                                                 
gAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTAt                                                                                                                                                 <  1:2917899/150‑1 (MQ=255)
 aaaaTAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATg                                                                                                                                                <  2:1250123/150‑1 (MQ=255)
 aaaaTAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATg                                                                                                                                                >  2:3944070/1‑150 (MQ=255)
   aataatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgc                                                                                                                                              <  1:4100124/150‑1 (MQ=255)
    ataatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgcg                                                                                                                                             >  2:2016707/1‑150 (MQ=255)
     taatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgcgc                                                                                                                                            >  2:3199348/1‑150 (MQ=255)
      aatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgcgcg                                                                                                                                           >  1:1039720/1‑150 (MQ=255)
      aatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgcgcg                                                                                                                                           >  1:1602972/1‑150 (MQ=255)
          gCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGcaca                                                                                                                                       >  1:2096966/1‑150 (MQ=255)
           ccAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGACCACAATGGCGCTGACGGTGGTTATGGCGCGCACaa                                                                                                                                      >  1:2515992/1‑150 (MQ=255)
           ccAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACaa                                                                                                                                      >  1:1975292/1‑150 (MQ=255)
            cAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACaaa                                                                                                                                     >  2:3994647/1‑150 (MQ=255)
              gCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGc                                                                                                                                   <  2:61191/150‑1 (MQ=255)
                aTCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTg                                                                                                                                 >  2:3934888/1‑150 (MQ=255)
                      acGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATc                                                                                                                           >  2:903139/1‑150 (MQ=255)
                       cGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCa                                                                                                                          >  1:18967/1‑150 (MQ=255)
                        gATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAg                                                                                                                         >  2:2668391/1‑150 (MQ=255)
                        gATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAg                                                                                                                         <  2:4521790/150‑1 (MQ=255)
                        gATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAg                                                                                                                         >  2:4578563/1‑150 (MQ=255)
                         aTCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgc                                                                                                                        >  1:2504693/1‑150 (MQ=255)
                         aTCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgc                                                                                                                        >  1:4556343/1‑150 (MQ=255)
                          tCAGAATGGAACCGTAGGCAATCGCTACACCGTTTTCGCCATCTTCCACGCTGTTGAGGATGTTAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACTTTGGCGCTGACGGTGGTTATGGCGAGCACAAAGCTGTAAATCAgcg                                                                                                                       <  2:1330075/150‑1 (MQ=255)
                          tCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgcg                                                                                                                       >  1:83762/1‑150 (MQ=255)
                          tCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgcg                                                                                                                       >  1:3673655/1‑150 (MQ=255)
                                   aaCCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGa                                                                                                              >  1:1467738/1‑150 (MQ=255)
                                   aaCCGTAGGCAATCGCTACACCGTATTCGCCATCGTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGa                                                                                                              <  1:1116093/150‑1 (MQ=255)
                                     ccGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATc                                                                                                            <  2:2697474/150‑1 (MQ=255)
                                      cGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCg                                                                                                           >  2:3598006/1‑150 (MQ=255)
                                               tCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCgca                                                                                                  >  1:4448972/1‑150 (MQ=255)
                                                  cTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCgcagca                                                                                               <  1:647200/150‑1 (MQ=255)
                                                   tACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGcagcag                                                                                              >  2:204195/1‑150 (MQ=255)
                                                   tACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGcagcag                                                                                              >  2:1785856/1‑150 (MQ=255)
                                                    acacCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCagcagc                                                                                             <  2:2796213/150‑1 (MQ=255)
                                                     cacCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCg                                                                                            <  1:2554157/150‑1 (MQ=255)
                                                      acCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCgg                                                                                           >  2:1146313/1‑150 (MQ=255)
                                                      acCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCgg                                                                                           >  2:2967758/1‑150 (MQ=255)
                                                      acCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCgg                                                                                           >  2:4519648/1‑150 (MQ=255)
                                                        cGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgca                                                                                         >  2:442016/1‑150 (MQ=255)
                                                         gTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgcag                                                                                        >  1:1289793/1‑150 (MQ=255)
                                                           aTTCGCCATCTTCCACGCGGAAGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgcagca                                                                                      >  2:2853096/1‑150 (MQ=255)
                                                             tCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAgg                                                                                    <  1:3699148/150‑1 (MQ=255)
                                                              cGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGa                                                                                   >  2:3498485/1‑150 (MQ=255)
                                                              cGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGa                                                                                   >  2:1084096/1‑150 (MQ=255)
                                                              cGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGa                                                                                   >  2:512796/1‑150 (MQ=255)
                                                               gCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAt                                                                                  >  2:1609410/1‑150 (MQ=255)
                                                                 cATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAtgt                                                                                >  1:40026/1‑150 (MQ=255)
                                                                 cATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAtgt                                                                                >  1:539983/1‑150 (MQ=255)
                                                                   tCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGc                                                                              <  2:1469083/150‑1 (MQ=255)
                                                                     ttCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGt                                                                            >  2:1770152/1‑150 (MQ=255)
                                                                     ttCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGt                                                                            >  2:2888757/1‑150 (MQ=255)
                                                                      tCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTa                                                                           <  2:3440485/150‑1 (MQ=255)
                                                                           gcgGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGt                                                                      <  2:2809797/150‑1 (MQ=255)
                                                                            cgGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTa                                                                     <  1:1370356/150‑1 (MQ=255)
                                                                            cgGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTa                                                                     <  2:4325671/150‑1 (MQ=255)
                                                                             ggTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTAc                                                                    <  1:49289/150‑1 (MQ=255)
                                                                               ttGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGc                                                                  >  1:2640300/1‑150 (MQ=255)
                                                                               ttGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGc                                                                  >  1:1705739/1‑150 (MQ=255)
                                                                                   ggATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAggg                                                              >  1:3606872/1‑150 (MQ=255)
                                                                                       gTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAAcc                                                          <  1:2389833/150‑1 (MQ=255)
                                                                                       gTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAAcc                                                          <  2:1699273/150‑1 (MQ=255)
                                                                                        tAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCg                                                         <  2:1914751/150‑1 (MQ=255)
                                                                                         aaGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGACCCgg                                                        >  1:2486788/1‑150 (MQ=255)
                                                                                              gTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGgcgcgc                                                   >  2:3480906/1‑150 (MQ=255)
                                                                                              gTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGgcgcgc                                                   >  2:3687786/1‑150 (MQ=255)
                                                                                               tCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCa                                                  <  2:3853236/150‑1 (MQ=255)
                                                                                                   tACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCt                                              <  2:2758725/150‑1 (MQ=255)
                                                                                                   tACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCt                                              <  1:916937/150‑1 (MQ=255)
                                                                                                    aCGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtg                                             >  1:202164/1‑150 (MQ=255)
                                                                                                    aCGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtg                                             >  1:469486/1‑150 (MQ=255)
                                                                                                       cgGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtgagt                                          <  2:4051003/150‑1 (MQ=255)
                                                                                                       cgGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtgagt                                          >  1:3572278/1‑150 (MQ=255)
                                                                                                            aTCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGc                                     >  2:1437654/1‑150 (MQ=255)
                                                                                                            aTCCGGCGTGAAGATGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACATAGCTGTTTATCCGCGCCGAGAGGATCGACGGGCGCAGCAGCGGCAGGAGGATGTGCGTAATGGTACGCTGGGTACCGGCGCGCAGGCTGTGTGAGGc                                     <  2:2175636/150‑1 (MQ=255)
                                                                                                               cGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTc                                  <  2:2132291/150‑1 (MQ=255)
                                                                                                                ggCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCa                                 >  1:2209903/1‑150 (MQ=255)
                                                                                                                 gCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCAt                                <  2:742923/150‑1 (MQ=255)
                                                                                                                  cGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATc                               >  1:322384/1‑150 (MQ=255)
                                                                                                                     gACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAg                            <  1:3651027/150‑1 (MQ=255)
                                                                                                                     gACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAg                            <  1:3651107/150‑1 (MQ=255)
                                                                                                                     gACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAg                            >  1:3763776/1‑150 (MQ=255)
                                                                                                                     gACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAg                            <  1:726430/150‑1 (MQ=255)
                                                                                                                      aCGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGc                           <  1:4296462/150‑1 (MQ=255)
                                                                                                                        gAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGa                         >  2:724721/1‑150 (MQ=255)
                                                                                                                        gAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGa                         >  2:4203151/1‑150 (MQ=255)
                                                                                                                         aGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGAt                        >  2:1101056/1‑150 (MQ=255)
                                                                                                                            aaCACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTa                     >  2:1488417/1‑150 (MQ=255)
                                                                                                                             acacAAGGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTAt                    <  2:3473254/150‑1 (MQ=255)
                                                                                                                                caATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGa                 >  2:3346716/1‑150 (MQ=255)
                                                                                                                                caATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGa                 <  1:3091219/150‑1 (MQ=255)
                                                                                                                                 aaTGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGAt                >  2:434778/1‑150 (MQ=255)
                                                                                                                                 aaTGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGAt                >  2:242019/1‑150 (MQ=255)
                                                                                                                                   tGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGAttt              <  2:3885229/150‑1 (MQ=255)
                                                                                                                                     gcgcTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTgg            <  2:3708041/150‑1 (MQ=255)
                                                                                                                                      cgcTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGAGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCTATTTGGc           >  1:2758891/1‑150 (MQ=255)
                                                                                                                                         tGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCg        <  1:3293423/150‑1 (MQ=255)
                                                                                                                                              gtggtTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGACCCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAAt   <  1:176252/150‑1 (MQ=255)
                                                                                                                                              gtggtTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAAt   <  1:4237194/150‑1 (MQ=255)
                                                                                                                                              gtggtTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAAt   <  1:3459944/150‑1 (MQ=255)
                                                                                                                                               tggtTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATg  <  2:1809818/150‑1 (MQ=255)
                                                                                                                                               tggtTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATg  <  2:2535791/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATG  >  CP009273/274588‑274880

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: