Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,652 G→T 9.5% pseudogene (2004/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,461,6520GT 9.5% 1001.8 / 63.3 290pseudogene (2004/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (115/147);  new base T (16/12);  total (131/159)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 293 total aligned reads displayed.

TCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGA  >  CP009273/1461503‑1461801
                                                                                                                                                     |                                                                                                                                                     
 cGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGa                                                                                                                                                      >  1:1703659/1‑150 (MQ=255)
 cGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGa                                                                                                                                                      >  2:1208411/1‑150 (MQ=255)
    tGGCAGGCAGTATGCAGGTTGGTGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAata                                                                                                                                                   >  2:4615441/1‑150 (MQ=255)
    tGGCAGGCAGTATGCAGGTTGGTGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAata                                                                                                                                                   >  2:2011364/1‑150 (MQ=255)
    tGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAata                                                                                                                                                   <  2:1358707/150‑1 (MQ=255)
      gcaggcagTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATc                                                                                                                                                 <  1:4252305/150‑1 (MQ=255)
      gcaggcagTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATc                                                                                                                                                 <  1:4795515/150‑1 (MQ=255)
        aggcagTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAATATCAc                                                                                                                                               >  2:4848778/1‑150 (MQ=255)
         ggcagTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACt                                                                                                                                              >  1:3546187/1‑150 (MQ=255)
         ggcagTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACt                                                                                                                                              >  2:3747308/1‑150 (MQ=255)
           cagTATGCAGGTTGGCGATTTCTCGACCGGGGTGGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTgg                                                                                                                                            <  2:4453495/150‑1 (MQ=255)
             gTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTa                                                                                                                                          <  2:908052/150‑1 (MQ=255)
               aTGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAAc                                                                                                                                        <  1:2475991/150‑1 (MQ=255)
                 gCAGGTTGGCGATTTCTCGACCGGGGTAGCTCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGt                                                                                                                                      <  2:635934/150‑1 (MQ=255)
                  cAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGtt                                                                                                                                     <  1:58853/150‑1 (MQ=255)
                    ggTTGTCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCt                                                                                                                                   <  2:4154750/150‑1 (MQ=255)
                        ggCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGACTATCACTGGTAACGTTCTGGtt                                                                                                                               >  2:5145082/1‑150 (MQ=255)
                        ggCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGtt                                                                                                                               <  1:1169755/150‑1 (MQ=255)
                        ggCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGtt                                                                                                                               <  1:2710804/150‑1 (MQ=255)
                           gATTTCTCGACCGGGGTAGTTCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGAt                                                                                                                            <  2:4304830/150‑1 (MQ=255)
                           gATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGAt                                                                                                                            >  1:866533/1‑150 (MQ=255)
                           gATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGAACTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGAt                                                                                                                            <  2:4346830/150‑1 (MQ=255)
                            aTTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATa                                                                                                                           >  1:2824252/1‑150 (MQ=255)
                              ttCTCGACCGTGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAACAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAg                                                                                                                         <  2:3199426/150‑1 (MQ=255)
                               tctcGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAgg                                                                                                                        <  2:2039833/150‑1 (MQ=255)
                                 tcGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGAt                                                                                                                      <  1:1562920/150‑1 (MQ=255)
                                      cGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAAc                                                                                                                 >  1:265193/1‑150 (MQ=255)
                                       ggggTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAAcc                                                                                                                >  1:1849316/1‑150 (MQ=255)
                                       ggggTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAAcc                                                                                                                <  2:3290869/150‑1 (MQ=255)
                                       ggggTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGGTCGGGTTGATAGGGCTAAAAcc                                                                                                                >  1:4472137/1‑150 (MQ=255)
                                          gTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGAACGAAAAGTGTTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCa                                                                                                             >  2:5008001/1‑150 (MQ=255)
                                           tAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAg                                                                                                            >  1:1512451/1‑150 (MQ=255)
                                                cTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAAt                                                                                                       >  1:3563739/1‑150 (MQ=255)
                                                       gCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGcggcgg                                                                                                <  1:2191156/150‑1 (MQ=255)
                                                        caacaaTAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGa                                                                                               >  1:2016735/1‑150 (MQ=255)
                                                           caaTAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAata                                                                                            <  1:2012076/150‑1 (MQ=255)
                                                           caaTAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAata                                                                                            >  2:2520472/1‑150 (MQ=255)
                                                            aaTAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAatat                                                                                           <  2:4093321/150‑1 (MQ=255)
                                                                aCAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATttct                                                                                       <  1:794976/150‑1 (MQ=255)
                                                                  aGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATttcttc                                                                                     <  1:1731571/150‑1 (MQ=255)
                                                                  aGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATttcttc                                                                                     >  2:1707827/1‑150 (MQ=255)
                                                                     gTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGat                                                                                  <  2:5136108/150‑1 (MQ=255)
                                                                      tGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTCCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGatc                                                                                 <  1:5148749/150‑1 (MQ=255)
                                                                            tGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCg                                                                           <  2:2880395/150‑1 (MQ=255)
                                                                             ggcggcAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTATTCGATCCATCCGt                                                                          >  1:4207538/1‑150 (MQ=255)
                                                                              gcggcAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTg                                                                         <  1:987668/150‑1 (MQ=255)
                                                                               cggcAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTgg                                                                        >  2:2399208/1‑150 (MQ=255)
                                                                                   aaaaGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTAt                                                                    >  1:2007003/1‑150 (MQ=255)
                                                                                   aaaaGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTAt                                                                    <  1:1296707/150‑1 (MQ=255)
                                                                                     aaGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCa                                                                  >  2:3971164/1‑150 (MQ=255)
                                                                                       gATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAAc                                                                >  1:2510972/1‑150 (MQ=255)
                                                                                        aTCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACg                                                               >  2:1574126/1‑150 (MQ=255)
                                                                                          cTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGtt                                                             <  1:3894942/150‑1 (MQ=255)
                                                                                            aaaaGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTa                                                           >  2:1794330/1‑150 (MQ=255)
                                                                                             aaaGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTAc                                                          <  1:3752138/150‑1 (MQ=255)
                                                                                               aGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACgg                                                        >  2:2314388/1‑150 (MQ=255)
                                                                                                 tGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCa                                                      <  2:2166600/150‑1 (MQ=255)
                                                                                                    tCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTg                                                   <  1:3146814/150‑1 (MQ=255)
                                                                                                       ggCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGata                                                >  2:1179271/1‑150 (MQ=255)
                                                                                                         cAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGataat                                              >  2:3443172/1‑150 (MQ=255)
                                                                                                           gAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGAtaataa                                            <  1:725346/150‑1 (MQ=255)
                                                                                                           gAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGAtaataa                                            <  2:3643061/150‑1 (MQ=255)
                                                                                                             aaGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACg                                          >  2:4502617/1‑150 (MQ=255)
                                                                                                             aaGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACg                                          >  2:4002994/1‑150 (MQ=255)
                                                                                                             aaGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACg                                          >  2:3893729/1‑150 (MQ=255)
                                                                                                               gCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTg                                        <  1:5239412/150‑1 (MQ=255)
                                                                                                                 gACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGAc                                      <  1:4924417/150‑1 (MQ=255)
                                                                                                                 gACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGAc                                      <  1:3633698/150‑1 (MQ=255)
                                                                                                                 gACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGGGGGTATCAACGTTTACGGCAGTGATAATAACGTGAc                                      >  2:712399/1‑150 (MQ=255)
                                                                                                                        aTAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCtggat                               >  1:902263/1‑150 (MQ=255)
                                                                                                                        aTAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCtggat                               <  1:2929925/150‑1 (MQ=255)
                                                                                                                        aTAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCtggat                               >  2:3555004/1‑150 (MQ=255)
                                                                                                                        aTAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCtggat                               <  2:3926952/150‑1 (MQ=255)
                                                                                                                         tAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCtggatg                              >  2:1267473/1‑150 (MQ=255)
                                                                                                                          aaaCGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCtggatgg                             >  2:1344840/1‑150 (MQ=255)
                                                                                                                           aaCGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCtggatgga                            <  2:2999169/150‑1 (MQ=255)
                                                                                                                            aCGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGaa                           <  1:3388923/150‑1 (MQ=255)
                                                                                                                             cGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGaaa                          <  2:3765549/150‑1 (MQ=255)
                                                                                                                               tttCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGt                        <  2:2646975/150‑1 (MQ=255)
                                                                                                                                ttCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGtt                       <  2:4803281/150‑1 (MQ=255)
                                                                                                                                 tCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAACACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTa                      <  1:3801781/150‑1 (MQ=255)
                                                                                                                                 tCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTa                      <  1:1029284/150‑1 (MQ=255)
                                                                                                                                  cTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTaa                     <  2:3321271/150‑1 (MQ=255)
                                                                                                                                    gtcGATGCGAATACAGTGAATATCCCTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACt                   <  2:2704318/148‑1 (MQ=255)
                                                                                                                                    ggCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACt                   >  2:3042860/1‑150 (MQ=255)
                                                                                                                                    ggCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACt                   <  2:877292/150‑1 (MQ=255)
                                                                                                                                       gATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACtgtt                >  1:2395752/1‑150 (MQ=255)
                                                                                                                                         tGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACtgttgt              <  2:3721911/150‑1 (MQ=255)
                                                                                                                                           cGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTAt            <  1:1908578/150‑1 (MQ=255)
                                                                                                                                            gAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGACAGTTAACtgtttgatc           >  2:4925513/1‑145 (MQ=255)
                                                                                                                                            gAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATc           >  1:3576180/1‑150 (MQ=255)
                                                                                                                                             aaTACAGTGAATATCACTTGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCa          <  2:1356143/150‑1 (MQ=255)
                                                                                                                                             aaTACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCa          <  2:1889827/150‑1 (MQ=255)
                                                                                                                                             aaTACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCa          >  2:1926984/1‑150 (MQ=255)
                                                                                                                                             aaTACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCa          >  1:1677420/1‑150 (MQ=255)
                                                                                                                                             aaTACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCa          >  1:2804164/1‑150 (MQ=255)
                                                                                                                                              aTACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAg         >  1:2255747/1‑150 (MQ=255)
                                                                                                                                                  aGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAg     <  2:1512451/150‑1 (MQ=255)
                                                                                                                                                   gTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGt    >  1:1323727/1‑150 (MQ=255)
                                                                                                                                                    tGAATATCACTGGTAACGTTCTGGGTGATAAGGATAAAACCGCAGACAATGCGGCGGAATCTTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTg   >  2:1093550/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGA  >  CP009273/1461503‑1461801

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: