Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,835,154 G→A 11.6% intergenic (‑55/+156) nepI ← / ← yicN putative transporter/DUF1198 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,835,1540GA11.6% 925.1 / 84.2 293intergenic (‑55/+156)nepI/yicNputative transporter/DUF1198 family protein
Reads supporting (aligned to +/- strand):  ref base G (110/148);  new base A (18/16);  total (129/164)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 305 total aligned reads displayed.

AGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTGATAAGACG  >  CP009273/3835005‑3835301
                                                                                                                                                     |                                                                                                                                                   
aGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAAtg                                                                                                                                                     >  2:3060317/1‑150 (MQ=255)
    aGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAAtgtgtg                                                                                                                                                 >  2:1282508/1‑150 (MQ=255)
    aGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAAtgtgtg                                                                                                                                                 >  1:309085/1‑150 (MQ=255)
     ggCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGa                                                                                                                                                >  2:1176215/1‑150 (MQ=255)
      gCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGAc                                                                                                                                               <  1:2935175/150‑1 (MQ=255)
          cacaAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCAt                                                                                                                                           <  2:3523648/150‑1 (MQ=255)
            caAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCAttt                                                                                                                                         <  1:29840/150‑1 (MQ=255)
            caAAACGCCACCGAGAAAACGGCTGAACAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCAttt                                                                                                                                         <  1:28699/150‑1 (MQ=255)
             aaaaCGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAACTTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTa                                                                                                                                        <  2:1702881/150‑1 (MQ=255)
             aaaaCGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGACATGTGTGACGCATTTa                                                                                                                                        >  2:2686567/1‑150 (MQ=255)
             aaaaCGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATATGTGACGCATTTa                                                                                                                                        >  1:1238808/1‑150 (MQ=255)
                 cGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGt                                                                                                                                    <  1:1781634/150‑1 (MQ=255)
                  gCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGtt                                                                                                                                   <  2:776050/150‑1 (MQ=255)
                  gCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGtt                                                                                                                                   <  2:725362/150‑1 (MQ=255)
                   ccACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGttt                                                                                                                                  <  1:2314885/150‑1 (MQ=255)
                    cACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGtttt                                                                                                                                 >  2:260761/1‑150 (MQ=255)
                        gagaAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCa                                                                                                                             <  1:3237389/150‑1 (MQ=255)
                        gagaAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATATGTGACGCATTTAACGTTTTTGCa                                                                                                                             <  2:287774/150‑1 (MQ=255)
                          gaAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAAc                                                                                                                           <  2:1660554/150‑1 (MQ=255)
                          gaAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAAc                                                                                                                           <  2:2349500/150‑1 (MQ=255)
                           aaaaCGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACt                                                                                                                          <  2:2928550/150‑1 (MQ=255)
                             aaCGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACttt                                                                                                                        <  1:452567/150‑1 (MQ=255)
                                   tGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATACATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGgt                                                                                                                  <  2:3779942/150‑1 (MQ=255)
                                   tGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGgt                                                                                                                  >  1:1310134/1‑150 (MQ=255)
                                   tGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGgt                                                                                                                  <  1:2447533/150‑1 (MQ=255)
                                    gACCAGTGCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGTCAATAAAGTCACTCCTGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACTCATTGAACGTTTTTGCAACTTTTACGgtg                                                                                                                 <  2:4070821/150‑1 (MQ=255)
                                     aCCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGgtgt                                                                                                                <  1:460445/150‑1 (MQ=255)
                                      ccAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGtt                                                                                                               <  1:4112685/150‑1 (MQ=255)
                                      ccAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGtt                                                                                                               <  1:727248/150‑1 (MQ=255)
                                       cAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTg                                                                                                              <  1:3822057/150‑1 (MQ=255)
                                        aGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGc                                                                                                             <  2:4342320/150‑1 (MQ=255)
                                         gTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCAGTTAACGTTTTTGCAACTTTTACGGTGTTGCt                                                                                                            >  2:3203848/1‑150 (MQ=255)
                                                 gggTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTgg                                                                                                    >  1:3360897/1‑150 (MQ=255)
                                                     gATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGtt                                                                                                >  2:1643646/1‑150 (MQ=255)
                                                     gATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGtt                                                                                                <  2:32343/150‑1 (MQ=255)
                                                      aTGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTATAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTa                                                                                               >  1:1434684/1‑150 (MQ=255)
                                                          cATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGAt                                                                                           <  2:2392155/150‑1 (MQ=255)
                                                          cATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGAt                                                                                           <  2:911554/150‑1 (MQ=255)
                                                               gcgccgcgGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGaaa                                                                                      <  2:1252646/150‑1 (MQ=255)
                                                               gcgccgcgGTTTTCGGCAAGAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGaaa                                                                                      <  2:136781/150‑1 (MQ=255)
                                                                 gccgcgGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACg                                                                                    >  1:1042284/1‑150 (MQ=255)
                                                                   cgcgGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGa                                                                                  <  2:4019584/150‑1 (MQ=255)
                                                                   cgcgGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGa                                                                                  >  1:1615604/1‑150 (MQ=255)
                                                                      ggTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACa                                                                               >  1:2539653/1‑150 (MQ=255)
                                                                             ggCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATATGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCaaa                                                                        >  2:1447657/1‑150 (MQ=255)
                                                                              gCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAg                                                                       >  2:2991299/1‑150 (MQ=255)
                                                                                aaTAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGAt                                                                     >  1:2847342/1‑150 (MQ=255)
                                                                                  tAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTa                                                                   >  2:797939/1‑150 (MQ=255)
                                                                                   aaaTTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTAc                                                                  >  1:645635/1‑150 (MQ=255)
                                                                                    aaTTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTAcc                                                                 >  2:926808/1‑150 (MQ=255)
                                                                                     aTTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATATGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCg                                                                >  2:2483780/1‑150 (MQ=255)
                                                                                      ttCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCgg                                                               >  2:2756955/1‑150 (MQ=255)
                                                                                         aCTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATg                                                            <  1:1167082/150‑1 (MQ=255)
                                                                                         aCTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATg                                                            <  1:3017745/150‑1 (MQ=255)
                                                                                          cTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATgc                                                           <  1:1915646/150‑1 (MQ=255)
                                                                                          cTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATgc                                                           <  1:7683/150‑1 (MQ=255)
                                                                                            cATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATgcgg                                                         <  2:2279400/150‑1 (MQ=255)
                                                                                               ggTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCTCTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGt                                                      >  2:3519450/1‑150 (MQ=255)
                                                                                               ggTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGt                                                      >  1:4300509/1‑150 (MQ=255)
                                                                                                gtgtTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTg                                                     <  2:4204568/150‑1 (MQ=255)
                                                                                                gtgtTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATATGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTg                                                     >  2:2799888/1‑150 (MQ=255)
                                                                                                 tgtTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGa                                                    <  1:2679916/150‑1 (MQ=255)
                                                                                                  gtTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGaa                                                   <  1:2153148/150‑1 (MQ=255)
                                                                                                   ttACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAAc                                                  >  2:3701085/1‑150 (MQ=255)
                                                                                                    tACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACg                                                 <  1:2030583/150‑1 (MQ=255)
                                                                                                     aCCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGt                                                <  2:1590906/150‑1 (MQ=255)
                                                                                                       ccGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCt                                              <  2:3864659/150‑1 (MQ=255)
                                                                                                       ccGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCt                                              <  2:3128575/150‑1 (MQ=255)
                                                                                                        cGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCtt                                             <  1:3697705/150‑1 (MQ=255)
                                                                                                         gTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTa                                            <  1:3434367/150‑1 (MQ=255)
                                                                                                          tGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTAt                                           <  1:981813/150‑1 (MQ=255)
                                                                                                           ggCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATATGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATc                                          >  1:4067193/1‑150 (MQ=255)
                                                                                                             cATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCa                                        >  2:1013701/1‑150 (MQ=255)
                                                                                                              aTGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCAc                                       >  2:2019396/1‑150 (MQ=255)
                                                                                                                gTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCAccc                                     <  1:3565158/150‑1 (MQ=255)
                                                                                                                 tttCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCg                                    <  2:1712355/150‑1 (MQ=255)
                                                                                                                     aGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAg                                <  2:254785/150‑1 (MQ=255)
                                                                                                                      ggggAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAgg                               >  2:513255/1‑150 (MQ=255)
                                                                                                                      ggggAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAgg                               >  2:392828/1‑150 (MQ=255)
                                                                                                                        ggAAAAGCCGACCACTTTAACATTGAAATATGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTc                             <  2:3320625/150‑1 (MQ=255)
                                                                                                                          aaaaGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAg                           <  2:2905308/150‑1 (MQ=255)
                                                                                                                             aGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGgcg                        >  1:963285/1‑150 (MQ=255)
                                                                                                                                cGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCgg                     >  1:3921258/1‑150 (MQ=255)
                                                                                                                                  aCCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTc                   <  2:1073943/150‑1 (MQ=255)
                                                                                                                                   ccACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCt                  <  2:2090172/150‑1 (MQ=255)
                                                                                                                                    cACTTTACCCTTGAACTGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTg                 <  2:3283822/150‑1 (MQ=255)
                                                                                                                                      cTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCCGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTa               >  2:2362648/1‑150 (MQ=255)
                                                                                                                                      cTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCCGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTa               >  2:855008/1‑150 (MQ=255)
                                                                                                                                        ttAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAgg             <  2:121239/150‑1 (MQ=255)
                                                                                                                                           aCATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTg          <  1:2389754/150‑1 (MQ=255)
                                                                                                                                            cATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTGa         >  1:4049262/1‑150 (MQ=255)
                                                                                                                                            cATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTGa         <  1:1566973/150‑1 (MQ=255)
                                                                                                                                             aTTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTGAt        <  1:2270082/150‑1 (MQ=255)
                                                                                                                                             aTTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTGAt        <  2:1047110/150‑1 (MQ=255)
                                                                                                                                             aTTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTGAt        >  1:840626/1‑150 (MQ=255)
                                                                                                                                              ttGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTGATa       <  2:1625253/150‑1 (MQ=255)
                                                                                                                                                gAAATATGTGACGCATTTAACTTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTGATAAg     <  1:3290112/150‑1 (MQ=255)
                                                                                                                                                  aaTGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTGATAAGAc   >  2:1793512/1‑150 (MQ=255)
                                                                                                                                                  aaTGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTGATAAGAc   >  2:3441927/1‑150 (MQ=255)
                                                                                                                                                   aTGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTGATAAGACg  <  1:2738713/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTACCCGTGGCATGTTTCAGGGGAAAAGCCGACCACTTTAACATTGAAATGTGTGACGCATTTAACGTTTTTGCAACTTTTACGGTGTTGCTTGAGCTGGAGTTAAGATCGAAACGGAACAGGGCAAAGATTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTGATAAGACG  >  CP009273/3835005‑3835301

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: