Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,315,800 A→G 15.4% intergenic (‑30/+56) 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,8000AG15.4% 335.9 / 50.0 149intergenic (‑30/+56)yciF/yciGputative rubrerythrin/ferritin‑like metal‑binding protein/KGG family protein
Reads supporting (aligned to +/- strand):  ref base A (69/57);  new base G (16/7);  total (85/64)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

Only 100 of 153 total aligned reads displayed.

ATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCT  >  NC_000913/1315701‑1315898
                                                                                                   |                                                                                                  
aTTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCa                                                                                                    >  1:342281/1‑100 (MQ=255)
aTTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCa                                                                                                    >  1:3324948/1‑100 (MQ=255)
aTTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCa                                                                                                    >  1:3194980/1‑100 (MQ=255)
aTTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCa                                                                                                    <  1:2845512/100‑1 (MQ=255)
 ttGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCAt                                                                                                   <  2:2025049/100‑1 (MQ=255)
  tGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATa                                                                                                  >  1:2582472/1‑100 (MQ=255)
  tGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATa                                                                                                  >  1:1218678/1‑100 (MQ=255)
  tGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATa                                                                                                  >  2:3154725/1‑100 (MQ=255)
  tGTTTTTCTGCGCTGTAGGCATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATa                                                                                                  >  1:624416/1‑100 (MQ=255)
   gTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATAc                                                                                                 <  1:3104441/100‑1 (MQ=255)
     ttttCTGCGCTGTAGGTATCTGAAAGCAGGTGAATACATACATCTTCTATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGGCATACgg                                                                                               <  2:3648323/100‑1 (MQ=255)
     ttttCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACgg                                                                                               <  1:3220201/100‑1 (MQ=255)
          tGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGc                                                                                          >  2:2817253/1‑100 (MQ=255)
           gcgcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCa                                                                                         <  1:635465/100‑1 (MQ=255)
            cgcTGTAGGTATCTGAACGCAGGTGAATAAATACATCTTCAATGGCCTTCATACTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGCCCCGc                                                                                          >  2:1725269/1‑98 (MQ=255)
            cgcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAg                                                                                        <  1:2079211/100‑1 (MQ=255)
            cgcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGGCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGcc                                                                                             >  2:165595/1‑95 (MQ=255)
             gcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTATTCATCTTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGa                                                                                       >  1:242620/1‑100 (MQ=255)
             gcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTATTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGc                                                                                          >  2:299063/1‑97 (MQ=255)
              cTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtg                                                                                      <  1:1010593/100‑3 (MQ=255)
              cTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtg                                                                                      <  1:133307/100‑3 (MQ=255)
              cTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgt                                                                                       <  2:2547389/99‑2 (MQ=255)
              cTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGAt                                                                                      <  1:1622957/100‑1 (MQ=255)
                 tAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtgatt                                                                                   <  1:834641/100‑6 (MQ=38)
                 tAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtgatt                                                                                   <  1:1778610/100‑6 (MQ=38)
                   ggTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtgatttc                                                                                 <  1:1155199/100‑8 (MQ=37)
                   ggTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAAt                                                                                 <  1:914286/100‑1 (MQ=255)
                   ggTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAAt                                                                                 <  2:525176/100‑1 (MQ=255)
                   ggTATCTGAAAGCAGGTGAATAAATACATATTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAAt                                                                                 <  1:2244595/100‑1 (MQ=255)
                     tATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTa                                                                               <  2:2980571/100‑1 (MQ=255)
                     tATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTa                                                                               <  2:1917664/100‑1 (MQ=255)
                        cTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCt                                                                            <  1:3604506/100‑1 (MQ=255)
                         tGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGTGCCATACGGCCCGCAGATTTAATTATCtt                                                                           >  2:3785194/1‑100 (MQ=255)
                          gAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                          >  2:843315/1‑100 (MQ=255)
                            aaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTAt                                                                        <  1:1815324/100‑1 (MQ=255)
                             aGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATc                                                                       <  1:1063711/100‑1 (MQ=255)
                             aGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATc                                                                       <  1:1292351/100‑1 (MQ=255)
                             aGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATc                                                                       <  1:2047528/100‑1 (MQ=255)
                                  gTGAATAAATATATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTTAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCt                                                                   <  2:1319259/99‑1 (MQ=255)
                                  gTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTg                                                                  <  1:1818858/100‑1 (MQ=255)
                                   tGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGa                                                                 <  1:304939/100‑1 (MQ=255)
                                   tGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGa                                                                 <  2:3689189/100‑1 (MQ=255)
                                    gAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGat                                                                >  2:1921509/1‑100 (MQ=255)
                                     aataaataCATCATCAAGGTGCTTCATATTCATATGTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCt                                                                   <  2:2462612/96‑1 (MQ=255)
                                       taaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATAtt                                                             <  2:2853344/100‑1 (MQ=255)
                                       taaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATAtt                                                             <  1:3744010/100‑1 (MQ=255)
                                           taCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaata                                                         <  2:2878486/100‑1 (MQ=255)
                                                  tCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACg                                                  >  1:2354273/1‑100 (MQ=255)
                                                    aaTGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGt                                                <  1:390512/100‑1 (MQ=255)
                                                         atctcatattcatatTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATc                                           <  2:2575167/97‑1 (MQ=255)
                                                          cTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCa                                          >  1:459532/1‑100 (MQ=255)
                                                          cTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCa                                          >  2:1480956/1‑100 (MQ=255)
                                                           ttcatattcataTTTTTCTCCAGTGAAATCACTGCGAACCATACGGTCCTCAGATGGAATTATCTTTATCAGCTGATATTAATAATTAACGGTGTATCa                                          >  2:3442844/1‑99 (MQ=255)
                                                            tcatattcatatTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATCCAATTATCTTTATTAGCTGATATTAATAATTAACGGTGAATCAg                                         >  1:561157/1‑99 (MQ=255)
                                                             cgtattcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAgat                                       >  1:2575186/3‑98 (MQ=255)
                                                              tctttcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAt                                      <  2:769370/97‑1 (MQ=255)
                                                              atattcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAt                                      >  2:109328/1‑100 (MQ=255)
                                                              atattcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAt                                      >  1:1460562/1‑100 (MQ=255)
                                                               tattcattttTTTCTCCAGTGAAATCACTGCGAGCCATAAGGCCCGAAGATGCAATTATCTTTATCTGCTGATATTAATAATTAACGGTGAATCAGGAtt                                     <  2:201803/100‑1 (MQ=255)
                                                               tattcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAtt                                     >  2:2877938/1‑100 (MQ=255)
                                                               tattcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTAAATCa                                          >  2:866446/1‑95 (MQ=255)
                                                                 ttcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAtaga                                   >  1:769378/1‑97 (MQ=255)
                                                                 ttcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTg                                   >  2:2826951/1‑100 (MQ=255)
                                                                 ttcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTg                                   >  2:273076/1‑100 (MQ=255)
                                                                 ttcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTg                                   >  2:273060/1‑100 (MQ=255)
                                                                  tcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGc                                  <  2:1600238/100‑1 (MQ=255)
                                                                   catattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGcc                                 <  1:2822177/100‑1 (MQ=255)
                                                                    atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCt                                >  1:1706323/1‑100 (MQ=255)
                                                                    atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCt                                >  2:3257950/1‑100 (MQ=255)
                                                                    atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCACAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCt                                >  2:3322427/1‑100 (MQ=255)
                                                                     tattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTg                               >  2:1314303/1‑100 (MQ=255)
                                                                           tctcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGGGAATCAGGATTTGCCTGATttat                         <  2:1658199/100‑1 (MQ=255)
                                                                            ctcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatt                        >  1:3629229/1‑100 (MQ=255)
                                                                             tcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatta                       >  1:2593547/1‑100 (MQ=255)
                                                                              ccAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTac                      >  1:3776819/1‑100 (MQ=255)
                                                                                aGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTacca                    >  2:2451540/1‑100 (MQ=255)
                                                                                 gTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTaccac                   >  2:491447/1‑100 (MQ=255)
                                                                                   gAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCg                 >  2:2145346/1‑100 (MQ=255)
                                                                                   gAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCg                 >  2:1002518/1‑100 (MQ=255)
                                                                                      aTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCtt              >  2:2816623/1‑100 (MQ=255)
                                                                                      aTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCtt              >  2:1844146/1‑100 (MQ=255)
                                                                                      aTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCtt              >  1:1255005/1‑100 (MQ=255)
                                                                                      aTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCtt              >  2:794992/1‑100 (MQ=255)
                                                                                      aTCACTGCGAGCCATACGGCCCGCAGAGGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCtt              <  1:299065/100‑1 (MQ=255)
                                                                                        cACTGCGAGCCATACGGTCCGTAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCTGGTTTTGCCTTAGTTATTACCACCGCTttg            <  2:45575/100‑1 (MQ=255)
                                                                                        cACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttg            <  2:459527/100‑1 (MQ=255)
                                                                                        atcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttg            >  1:2716091/3‑100 (MQ=255)
                                                                                         tcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgt           >  2:2197118/2‑100 (MQ=255)
                                                                                         tcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttg            >  1:2976428/2‑99 (MQ=255)
                                                                                          cTGCGGGCCGTATGGCTCGCAGATGCAATTATATTAATCAGCTGATATTGATAATTAACGGTGAAGCAGGCTTTGCCTGATTAATTACCACAACTttgtt          >  2:1258249/1‑100 (MQ=255)
                                                                                          cTGCGGGCCGTATGGCGCGCAGATGCAATTATCTTTATCAGGTGATATTACTAATAAACGGCGAAGCAGGATTTGACTGAGTTATTACCACCGCTtcgtt          >  2:953080/1‑100 (MQ=255)
                                                                                           tGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTATCACCGCTttgttg         >  2:661967/1‑100 (MQ=255)
                                                                                           tGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTAGCACCGCTttgtt          >  2:1197935/1‑99 (MQ=255)
                                                                                           tGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgttg         <  2:57025/100‑1 (MQ=255)
                                                                                           tGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgttg         <  1:3437383/100‑1 (MQ=255)
                                                                                            gCGAGCCATACGGCCCGCAGAAGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATATGCCTGATTTATTACCACCGCTTTGTTGa        >  2:231537/1‑100 (MQ=255)
                                                                                              gAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGAcc      <  1:1542620/100‑1 (MQ=255)
                                                                                               aGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCg     >  1:297210/1‑100 (MQ=255)
                                                                                                  cATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCt  >  1:3802097/1‑100 (MQ=255)
                                                                                                  cATACGGCCCGCAGATGCAATTATCTTTATCAGATGATATTAATAATTAAAGGTGAATGAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCt  >  1:3132934/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
ATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCT  >  NC_000913/1315701‑1315898

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: