Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,315,796 A→G 6.6% intergenic (‑26/+60) yciF ← / ← yciG putative rubrerythrin/ferritin‑like metal‑binding protein/KGG family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,315,7960AG 6.6% 492.3 / 15.7 148intergenic (‑26/+60)yciF/yciGputative rubrerythrin/ferritin‑like metal‑binding protein/KGG family protein
Reads supporting (aligned to +/- strand):  ref base A (65/73);  new base G (2/8);  total (67/81)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 152 total aligned reads displayed.

ACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAC  >  NC_000913/1315648‑1315944
                                                                                                                                                    |                                                                                                                                                    
aCTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGa                                                                                                                                                      <  3:121211/149‑1 (MQ=255)
aCTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGa                                                                                                                                                      <  3:1330870/149‑1 (MQ=255)
aCTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGa                                                                                                                                                      <  2:204582/149‑1 (MQ=255)
  ttAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGc                                                                                                                                                    >  3:1265602/1‑149 (MQ=255)
   tAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGcc                                                                                                                                                   <  2:422366/149‑1 (MQ=255)
    aaTTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCa                                                                                                                                                  <  2:122322/149‑1 (MQ=255)
      tttttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATa                                                                                                                                                <  2:182058/149‑1 (MQ=255)
      tttttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATa                                                                                                                                                >  1:848872/1‑149 (MQ=255)
       ttttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATAc                                                                                                                                               <  4:541495/149‑1 (MQ=255)
        tttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGTGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACg                                                                                                                                              >  1:327271/1‑149 (MQ=255)
             tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTATTCTCCAGTGAAATCACTGCGAGCCATACGGCCCg                                                                                                                                         >  4:952366/1‑149 (MQ=255)
             tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCACGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGGCATCATATTCATATTTTTCTCCAGTGAAATCACAGCGAGCCATACGGCCCg                                                                                                                                         >  4:912898/1‑149 (MQ=255)
                   gTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtga                                                                                                                                   <  3:1269856/149‑4 (MQ=255)
                    ttGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGc                                                                                                                                  >  2:900121/1‑149 (MQ=255)
                    ttGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGc                                                                                                                                  <  1:1231505/149‑1 (MQ=255)
                          ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCaa                                                                                                                                <  1:506301/145‑1 (MQ=255)
                          ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCaa                                                                                                                                >  2:506301/1‑145 (MQ=255)
                          ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTAt                                                                                                                            <  2:1322305/149‑1 (MQ=255)
                          ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTAt                                                                                                                            >  2:704970/1‑149 (MQ=255)
                           tGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCACTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATc                                                                                                                           <  3:952366/149‑1 (MQ=255)
                              gAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                        <  3:256711/149‑1 (MQ=255)
                              gAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                        >  4:377833/1‑149 (MQ=255)
                              gAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                        >  3:1340266/1‑149 (MQ=255)
                              gAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                        >  2:231175/1‑149 (MQ=255)
                              gAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                        >  4:146038/1‑149 (MQ=255)
                                gTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTAt                                                                                                                      <  2:1368380/149‑1 (MQ=255)
                                gTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTAt                                                                                                                      <  1:723431/149‑1 (MQ=255)
                                  tttGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCa                                                                                                                    >  1:948388/1‑149 (MQ=255)
                                  tttGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCa                                                                                                                    <  3:773313/149‑1 (MQ=255)
                                      ccAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTg                                                                                                                >  3:721927/1‑149 (MQ=255)
                                         gTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgt                                                                                                                                     <  2:893456/125‑2 (MQ=255)
                                         gTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGata                                                                                                             <  3:929805/149‑1 (MQ=255)
                                         gTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGata                                                                                                             <  3:688107/149‑1 (MQ=255)
                                         gTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGata                                                                                                             >  1:356893/1‑149 (MQ=255)
                                              cGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaat                                                                                                        <  2:446809/149‑1 (MQ=255)
                                               gggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTTATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaata                                                                                                       >  3:926873/1‑149 (MQ=255)
                                                ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaataa                                                                                                      <  3:1191450/149‑1 (MQ=255)
                                                ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaataa                                                                                                      >  4:382571/1‑149 (MQ=255)
                                                ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaataa                                                                                                      <  1:851752/149‑1 (MQ=255)
                                                     aTTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAAc                                                                                                 <  1:221247/149‑1 (MQ=255)
                                                     aTTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAAc                                                                                                 <  4:1340266/149‑1 (MQ=255)
                                                     aTTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAAc                                                                                                 <  3:108983/149‑1 (MQ=255)
                                                       tGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACgg                                                                                               <  3:1402773/149‑1 (MQ=255)
                                                        gTTTTTCTGCGCTGTAGGTATCTGAGAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGt                                                                                              <  3:955403/149‑1 (MQ=255)
                                                        gTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGt                                                                                              <  3:227598/149‑1 (MQ=255)
                                                         tttttCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTg                                                                                             <  1:1363574/149‑1 (MQ=255)
                                                           tttCTGCGCTGTAGGTATCTGAAAGCGGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGaa                                                                                           <  2:1150463/149‑1 (MQ=255)
                                                                  gcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAt                                                                                    <  3:382571/149‑1 (MQ=255)
                                                                      tAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGc                                                                                >  1:1112364/1‑149 (MQ=255)
                                                                        ggTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCt                                                                              >  3:1424536/1‑149 (MQ=255)
                                                                           aTCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAt                                                                           <  1:356697/149‑1 (MQ=255)
                                                                            tCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAtt                                                                          <  3:1077421/149‑1 (MQ=255)
                                                                            tCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAtt                                                                          >  1:518005/1‑149 (MQ=255)
                                                                            tCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAtt                                                                          <  3:657403/149‑1 (MQ=255)
                                                                                aaaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatt                                                                      <  2:337153/149‑1 (MQ=255)
                                                                                aaaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatt                                                                      <  2:1475338/149‑1 (MQ=255)
                                                                                 aaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatta                                                                     >  1:1236849/1‑149 (MQ=255)
                                                                                 aaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatta                                                                     <  4:119271/149‑1 (MQ=255)
                                                                                  aGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTac                                                                    <  4:347523/149‑1 (MQ=255)
                                                                                   gCAGGTGAATAAATACATCTTCAATGGCCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTacc                                                                   >  2:374085/1‑149 (MQ=255)
                                                                                    cAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTacca                                                                  <  2:968238/149‑1 (MQ=255)
                                                                                           ataaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCttt                                                           <  1:1053421/149‑1 (MQ=255)
                                                                                           ataaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCttt                                                           <  4:810902/149‑1 (MQ=255)
                                                                                            taaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttg                                                          <  4:914589/149‑1 (MQ=255)
                                                                                              aataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTCGCCTGATTTATTACCACCGCTttgtt                                                        <  4:49944/149‑1 (MQ=255)
                                                                                               ataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                        >  2:1245915/1‑84 (MQ=255)
                                                                                                 aCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGAc                                                     <  1:900706/149‑1 (MQ=255)
                                                                                                   aTCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCg                                                   >  2:519626/1‑149 (MQ=255)
                                                                                                   aTCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCg                                                   >  3:471809/1‑149 (MQ=255)
                                                                                                      ttCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCTTGATTTATTACCACCGCTTGGTTGACCGCCt                                                >  1:301401/1‑149 (MQ=255)
                                                                                                      ttCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatta                                                                     <  2:208835/128‑1 (MQ=255)
                                                                                                      ttCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatta                                                                     >  1:208835/1‑128 (MQ=255)
                                                                                                      ttCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCt                                                >  1:1307902/1‑149 (MQ=255)
                                                                                                      ttCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCAGATTTATTACCACCGCTTTGTTGACCGCCt                                                >  1:959487/1‑149 (MQ=255)
                                                                                                         aaTGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttt                                             <  1:742409/149‑1 (MQ=255)
                                                                                                          aTGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCttttt                                            <  4:471809/149‑1 (MQ=255)
                                                                                                             gTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAt                                                                           <  4:578276/115‑1 (MQ=255)
                                                                                                             gTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAt                                                                           >  3:578276/1‑115 (MQ=255)
                                                                                                                  catattcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCt                                    <  4:171071/149‑1 (MQ=255)
                                                                                                                   atattcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCtt                                   >  3:1442521/1‑149 (MQ=255)
                                                                                                                     attcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                 >  3:828283/1‑149 (MQ=255)
                                                                                                                         atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATg                             >  3:1358223/1‑149 (MQ=255)
                                                                                                                          tattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgc                            >  2:828315/1‑149 (MQ=255)
                                                                                                                          tattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgc                            >  1:684183/1‑149 (MQ=255)
                                                                                                                            tttttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgcgc                          <  1:224011/149‑1 (MQ=255)
                                                                                                                                 ctcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGc                     <  4:721927/149‑1 (MQ=255)
                                                                                                                                 ctcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGc                     >  4:271506/1‑149 (MQ=255)
                                                                                                                                     aGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGAt                 >  3:16124/1‑149 (MQ=255)
                                                                                                                                      gTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGa                  >  4:11204/1‑147 (MQ=255)
                                                                                                                                         aaaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAtt             <  4:16124/149‑1 (MQ=255)
                                                                                                                                         aaaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAtt             <  2:301035/149‑1 (MQ=255)
                                                                                                                                          aaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttt            >  1:184044/1‑149 (MQ=255)
                                                                                                                                          aaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttt            <  3:599762/149‑1 (MQ=255)
                                                                                                                                            tCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttt                                             <  4:317275/114‑1 (MQ=255)
                                                                                                                                            tCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCtt                                   >  1:432561/1‑124 (MQ=255)
                                                                                                                                            tCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCtt                                   <  2:432561/124‑1 (MQ=255)
                                                                                                                                             cACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa         <  2:48135/149‑1 (MQ=255)
                                                                                                                                              tcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                                                        >  3:1358860/2‑101 (MQ=255)
                                                                                                                                              tcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                                                        <  4:1358860/100‑1 (MQ=255)
                                                                                                                                                  cGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAtt    >  2:1162405/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAC  >  NC_000913/1315648‑1315944

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