Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 1,915,165 T→C intergenic (+54/+14) RPPX_08535 → / ← RPPX_08540 ferrioxamine B receptor/5S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099741,915,1650TC99.6% 643.4 / ‑3.2 222intergenic (+54/+14)RPPX_08535/RPPX_08540ferrioxamine B receptor/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (92/129);  minor base . (0/1);  total (92/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 226 total aligned reads displayed.

CCAATGAAAGCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGT  >  CP009974/1915017‑1915303
                                                                                                                                                    |                                                                                                                                          
ccAATGAAAGCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGctt                                                                                                                                          >  2:5822861/1‑151 (MQ=255)
 cAATGAAATCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGcttt                                                                                                                                         <  2:5518771/151‑1 (MQ=255)
  aaTGAAAGCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTg                                                                                                                                        <  2:1082926/151‑1 (MQ=255)
    aaaaaGCTCCATTACCTCATGTGCCAGCCTCAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTACCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTggg                                                                                                                                      <  2:879130/149‑3 (MQ=255)
         gCTACATTACCTCATGTGCCAGCCTGAACTACTGTTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTgggtttga                                                                                                                                 <  1:5320385/151‑8 (MQ=255)
         gCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTgggtttga                                                                                                                                 <  1:2066725/151‑8 (MQ=255)
                   cTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGa                                                                                                                      >  1:3877566/1‑151 (MQ=255)
                    tCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGAt                                                                                                                     >  1:764299/1‑151 (MQ=255)
                       tgtgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGAc                                                                                                                  >  2:3783672/1‑151 (MQ=255)
                       tgtgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGAc                                                                                                                  <  1:6461780/151‑1 (MQ=255)
                       tgtgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGAc                                                                                                                  >  1:3728235/1‑151 (MQ=255)
                        gtgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGAcc                                                                                                                 <  2:199157/151‑1 (MQ=255)
                        gtgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGAcc                                                                                                                 >  1:2322556/1‑151 (MQ=255)
                           ccAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTAc                                                                                                              <  1:1055153/151‑1 (MQ=255)
                           ccAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTAc                                                                                                              >  1:3267636/1‑151 (MQ=255)
                            cAGCGTGACCTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAACCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTAct                                                                                                             <  2:4248729/151‑1 (MQ=255)
                            cAGCCTGAACTACGGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTAct                                                                                                             <  2:6366618/151‑1 (MQ=255)
                              gCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTActct                                                                                                           >  2:5244365/1‑151 (MQ=255)
                              gCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTActct                                                                                                           >  2:5700128/1‑151 (MQ=255)
                              gCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTActct                                                                                                           >  1:4343739/1‑151 (MQ=255)
                                    actactGCTAGATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATg                                                                                                     <  2:2863024/151‑1 (MQ=255)
                                    actactGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATg                                                                                                     >  1:2305638/1‑151 (MQ=39)
                                    actactGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATg                                                                                                     <  2:2310255/151‑1 (MQ=255)
                                     ctactGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATgg                                                                                                    >  2:1784643/1‑151 (MQ=255)
                                      tactGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATggg                                                                                                   >  1:5175107/1‑151 (MQ=39)
                                      tactGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATggg                                                                                                   <  1:2500266/151‑1 (MQ=255)
                                       actGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATgggg                                                                                                  >  2:1917891/1‑151 (MQ=39)
                                       actGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATgggg                                                                                                  <  1:5883831/151‑1 (MQ=255)
                                              cATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAAcccc                                                                                           <  1:2190191/151‑1 (MQ=38)
                                              cATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAAcccc                                                                                           <  2:342012/151‑1 (MQ=38)
                                              cATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAAcccc                                                                                           >  1:5971423/1‑151 (MQ=38)
                                               aTGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCca                                                                                          >  2:3478016/1‑151 (MQ=38)
                                               aTGGGTGAAGAGCACAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCca                                                                                          <  2:3977469/151‑1 (MQ=38)
                                                tGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCcac                                                                                         >  1:1613302/1‑151 (MQ=38)
                                                 gggTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCcaca                                                                                        >  1:2230223/1‑151 (MQ=38)
                                                    tGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAGACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTa                                                                                     >  2:6114591/1‑151 (MQ=37)
                                                          gcgcAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCg                                                                               <  2:1857896/151‑1 (MQ=37)
                                                          gcgcAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCg                                                                               <  1:4281448/151‑1 (MQ=37)
                                                           ctatgtgtgAGCGCCACGGTCAGCTACCAGTTCTACCCGGTATTGAAACCCGTGAAAGACAACGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCgg                                                                              <  2:841416/146‑1 (MQ=21)
                                                           cgcAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCgg                                                                              >  2:6236564/1‑151 (MQ=37)
                                                            gcAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGc                                                                             >  1:3692711/1‑151 (MQ=37)
                                                             cAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCg                                                                            <  2:4703516/151‑1 (MQ=37)
                                                             cAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCg                                                                            <  2:4967376/151‑1 (MQ=37)
                                                               aTGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGAt                                                                          >  1:2029274/1‑151 (MQ=37)
                                                               aTGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGAt                                                                          <  1:4583526/151‑1 (MQ=37)
                                                                 gTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGc                                                                        >  1:841416/1‑151 (MQ=25)
                                                                    aGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATc                                                                     <  2:6430318/151‑1 (MQ=25)
                                                                      cgcTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGt                                                                   >  2:3033395/1‑151 (MQ=25)
                                                                      cgcTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGt                                                                   <  1:5665958/151‑1 (MQ=25)
                                                                       gcTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGtt                                                                  <  2:5293389/151‑1 (MQ=25)
                                                                        cTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGttt                                                                 <  1:1798325/151‑1 (MQ=25)
                                                                          aCGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGc                                                                        >  2:5317814/1‑142 (MQ=25)
                                                                            ggggaGCGATCAGTGCTGACCGATGTTGAAAGCCGTGAAAGACAGAGCCCCGGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCact                                                             <  2:4321465/147‑1 (MQ=18)
                                                                            ggTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCact                                                             >  2:148785/1‑151 (MQ=25)
                                                                              tgaggtATCAGTTCTACCCGATATTGCAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCactac                                                           <  2:392979/146‑1 (MQ=18)
                                                                                aGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCActgctg                                                         <  2:415695/151‑1 (MQ=21)
                                                                                  cTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAg                                                       <  1:2575458/151‑1 (MQ=21)
                                                                                     tCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTc                                                    >  1:4213362/1‑151 (MQ=21)
                                                                                      cAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCg                                                   >  1:6111055/1‑151 (MQ=21)
                                                                                      cAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCg                                                   >  1:3629810/1‑151 (MQ=21)
                                                                                      cAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCg                                                   >  2:350250/1‑151 (MQ=21)
                                                                                         ttCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATAGTTTCACTGCTGAGTTCggaa                                                >  2:4217545/1‑149 (MQ=18)
                                                                                          tCTAACCGATATTGAAACCCGTGAAAGACAACGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggat                                               <  2:4636445/151‑1 (MQ=18)
                                                                                          tCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggat                                               <  2:2322024/151‑1 (MQ=18)
                                                                                          tCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggat                                               >  2:2046247/1‑151 (MQ=18)
                                                                                           cTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGTCGATGCATCGTTTCACTGCTGTGTTCgggatg                                              >  1:2960551/1‑151 (MQ=18)
                                                                                           cTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggttg                                              >  1:2412721/1‑151 (MQ=18)
                                                                                           cTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatg                                              <  1:1178540/151‑1 (MQ=18)
                                                                                            tAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatgg                                             >  2:6077551/1‑151 (MQ=18)
                                                                                                   tatTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAgg                                      <  1:864815/151‑1 (MQ=18)
                                                                                                      tGAATCCCGTGAAAGACAACGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAggtgg                                   <  2:3281664/151‑1 (MQ=17)
                                                                                                       gAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggtggt                                  <  1:799489/151‑1 (MQ=17)
                                                                                                       gAAACCCGGGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggtggt                                  <  2:5070526/151‑1 (MQ=17)
                                                                                                        aataccGTGAGAGACAAAGCCCCTGTGTGCAGGGCAGACGGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGtt                                 <  2:6007521/147‑1 (MQ=14)
                                                                                                        aaaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGtt                                 <  2:504917/151‑1 (MQ=17)
                                                                                                         aaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTc                                >  1:2421529/1‑151 (MQ=17)
                                                                                                          aCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTcc                               <  1:3449864/151‑1 (MQ=17)
                                                                                                          aCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTcc                               <  1:6264891/151‑1 (MQ=17)
                                                                                                            ccGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCaa                             <  1:147716/151‑1 (MQ=17)
                                                                                                              gTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATg                           <  2:788174/151‑1 (MQ=14)
                                                                                                               tGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACCCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGCTGGGATCAGGTGGTTCCAATGc                          >  1:3345240/1‑151 (MQ=14)
                                                                                                                gAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGct                         >  1:4978716/1‑151 (MQ=14)
                                                                                                                   aGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTa                      <  1:1863991/151‑1 (MQ=14)
                                                                                                                   aGACAAAGCCCCTGTCAGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTa                      <  2:1831899/151‑1 (MQ=14)
                                                                                                                     aCAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATg                    <  1:573572/151‑1 (MQ=14)
                                                                                                                      cAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATgg                   <  2:5335596/151‑1 (MQ=14)
                                                                                                                        aaGCCCCTGTCTGCATCGCAGACAGGGGCTTGGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtc                 <  1:6435073/151‑1 (MQ=14)
                                                                                                                          gCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTCCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTCTGgtcgt               >  2:4015583/1‑151 (MQ=14)
                                                                                                                            cccTTTCTGCATCGCAGACAGGGGCTTTGGGTTTTAATCTTGACGATGACCTACTCTCACATGGGTAAACCCCACACTACCATCGGCGATGCATCGGTTCACTGCTGAGTTCGGGATGTGATCAGGTGGTTCCAATGCTCTATGATCGTCa             >  2:512003/1‑151 (MQ=11)
                                                                                                                            cccTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             >  1:5122179/1‑151 (MQ=14)
                                                                                                                            cccTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             >  2:216255/1‑151 (MQ=14)
                                                                                                                            cccTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             >  1:6020359/1‑151 (MQ=14)
                                                                                                                             ccTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCaa            >  2:5158717/1‑151 (MQ=14)
                                                                                                                              cTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAg           >  2:3667526/1‑151 (MQ=14)
                                                                                                                                  cTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAAt       <  2:6475439/151‑1 (MQ=11)
                                                                                                                                  cTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAAt       <  2:72938/151‑1 (MQ=11)
                                                                                                                                      aTCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTc     <  2:2327629/149‑1 (MQ=11)
                                                                                                                                       tCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGCCCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGt  <  2:438245/151‑1 (MQ=11)
                                                                                                                                       tCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGt  >  1:570204/1‑151 (MQ=11)
                                                                                                                                       tCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGt  <  2:2707024/151‑1 (MQ=11)
                                                                                                                                                    |                                                                                                                                          
CCAATGAAAGCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGT  >  CP009974/1915017‑1915303

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: