Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,315,807 C→T 16.0% intergenic (‑37/+49) 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,8070CT16.0% 310.5 / 43.8 139intergenic (‑37/+49)yciF/yciGputative rubrerythrin/ferritin‑like metal‑binding protein/KGG family protein
Reads supporting (aligned to +/- strand):  ref base C (64/52);  new base T (15/7);  total (80/59)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

Only 100 of 143 total aligned reads displayed.

CTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACC  >  NC_000913/1315710‑1315906
                                                                                                 |                                                                                                   
 tGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGc                                                                                                  >  2:2817253/1‑100 (MQ=255)
 tGCGCTGTAGGTAACTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGc                                                                                                  >  1:2590167/1‑100 (MQ=255)
 tGCGCTGTAAGTATCTGAAAGTAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGc                                                                                                  <  2:1848709/100‑1 (MQ=255)
  gcgcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCa                                                                                                 <  1:635465/100‑1 (MQ=255)
   cgcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAg                                                                                                <  1:2079211/100‑1 (MQ=255)
   cgcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGGCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGcc                                                                                                     >  2:165595/1‑95 (MQ=255)
     cTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtg                                                                                              <  1:133307/100‑3 (MQ=255)
     cTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtg                                                                                              <  1:1010593/100‑3 (MQ=255)
     cTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgt                                                                                               <  2:2547389/99‑2 (MQ=255)
     cTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGAt                                                                                              <  1:1622957/100‑1 (MQ=255)
        tAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtgatt                                                                                           <  1:1778610/100‑6 (MQ=38)
         aGGTATCTGAAAGCAGGTGAATAAACACATCTCCAATGTTCTTCATCTTCATAGTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCaa                                                                                          <  2:2110922/100‑1 (MQ=255)
          ggTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAAt                                                                                         <  2:525176/100‑1 (MQ=255)
            tATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTa                                                                                       <  2:2980571/100‑1 (MQ=255)
             aTCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTAt                                                                                      >  2:2864856/1‑100 (MQ=255)
               cTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCACTTATCt                                                                                    >  1:1251460/1‑100 (MQ=255)
               cTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCt                                                                                    <  1:3604506/100‑1 (MQ=255)
               cTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCt                                                                                    >  2:53035/1‑100 (MQ=255)
               cTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATGCATATTTTTCTCCAGTGAAATTACTGCGAGCCATACGTCCCGCCGATGCAATTATCt                                                                                    >  2:889561/1‑100 (MQ=255)
                tGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGTGCCATACGGCCCGCAGATTTAATTATCtt                                                                                   >  2:3785194/1‑100 (MQ=255)
                 gAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                  >  2:480295/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:2047528/100‑1 (MQ=255)
                     gCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCa                                                                              <  1:878866/100‑1 (MQ=255)
                         gTGAATAAATATATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTTAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCt                                                                           <  2:1319259/99‑1 (MQ=255)
                         gTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTg                                                                          <  1:1818858/100‑1 (MQ=255)
                         gTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTg                                                                          <  2:689577/100‑1 (MQ=255)
                          tGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGa                                                                         <  1:304939/100‑1 (MQ=255)
                           gAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGat                                                                        >  2:1921509/1‑100 (MQ=255)
                            aataaataCATCATCAAGGTGCTTCATATTCATATGTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCt                                                                           <  2:2462612/96‑1 (MQ=255)
                             ataaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGatat                                                                      >  2:1010105/1‑100 (MQ=255)
                              taaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATAtt                                                                     <  2:2853344/100‑1 (MQ=255)
                              taaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATAtt                                                                     <  1:3744010/100‑1 (MQ=255)
                                aataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaa                                                                   >  2:2247350/1‑100 (MQ=255)
                                 ataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaat                                                                  <  1:2753250/100‑1 (MQ=255)
                                  taCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaata                                                                 <  1:3418225/100‑1 (MQ=255)
                                    cATCTTCAATGGTCTTCATATTCATATTTTTCTCAAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTaataat                                                               <  2:3067054/100‑1 (MQ=255)
                                         tCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACg                                                          >  1:2354273/1‑100 (MQ=255)
                                           aaTGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGt                                                        <  1:390512/100‑1 (MQ=255)
                                                tCTTCATATTCATATTTTTCTCTAGTGTATTCACTGCGAACCATACGACCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTTAAt                                                    >  2:2894596/1‑99 (MQ=255)
                                                atctcatattcatatTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATc                                                   <  2:2575167/97‑1 (MQ=255)
                                                 cTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCa                                                  >  2:1480956/1‑100 (MQ=255)
                                                  ttcatattcataTTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAg                                                 <  1:2084351/100‑1 (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)
                                                     atattcttattTTTCTCCAGTGAAATCATTGCTAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGa                                               <  1:2341825/99‑1 (MQ=255)
                                                     atattcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAt                                              >  2:109328/1‑100 (MQ=255)
                                                      tattcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAtt                                             >  2:2877938/1‑100 (MQ=255)
                                                        ttcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAtaga                                           >  1:769378/1‑97 (MQ=255)
                                                        ttcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTg                                           >  2:2626738/1‑100 (MQ=255)
                                                        ttcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTg                                           >  2:273076/1‑100 (MQ=255)
                                                         tcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGc                                          <  2:1600238/100‑1 (MQ=255)
                                                          catattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGcc                                         <  1:2822177/100‑1 (MQ=255)
                                                           atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCt                                        >  2:3257950/1‑100 (MQ=255)
                                                           atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCt                                        >  1:1706323/1‑100 (MQ=255)
                                                           atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCACAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCt                                        >  2:3322427/1‑100 (MQ=255)
                                                            tattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTg                                       >  2:1314303/1‑100 (MQ=255)
                                                                 ttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATtta                                  <  1:1488530/100‑1 (MQ=255)
                                                                  tctcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGGGAATCAGGATTTGCCTGATttat                                 <  2:1658199/100‑1 (MQ=255)
                                                                   ctcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatt                                <  1:3154749/100‑1 (MQ=255)
                                                                   ctcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatt                                >  1:3629229/1‑100 (MQ=255)
                                                                     ccAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTac                              >  1:3776819/1‑100 (MQ=255)
                                                                     ccAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTac                              >  2:3795311/1‑100 (MQ=255)
                                                                       aGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTacca                            >  2:2451540/1‑100 (MQ=255)
                                                                          gAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCg                         >  2:1002518/1‑100 (MQ=255)
                                                                          gAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCg                         >  2:2145346/1‑100 (MQ=255)
                                                                            tgcatcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCt                       >  1:1324079/6‑100 (MQ=38)
                                                                             aTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCtt                      >  2:1844146/1‑100 (MQ=255)
                                                                             aTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCtt                      >  1:653178/1‑100 (MQ=255)
                                                                             aTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCtt                      >  1:1255005/1‑100 (MQ=255)
                                                                             aTCACTGCGAGCCATACGGCCCGCAGAGGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCtt                      <  1:299065/100‑1 (MQ=255)
                                                                              catcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCctt                     >  2:2819952/4‑100 (MQ=255)
                                                                               cACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttg                    <  2:459527/100‑1 (MQ=255)
                                                                               atcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttg                    >  1:1103508/3‑100 (MQ=255)
                                                                               atcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttg                    >  1:2716091/3‑100 (MQ=255)
                                                                                tcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATCTTAATAATTAACGGTGAATCAGGATTTGCTTGATTTATTACCACCGCTttgt                   >  2:200534/2‑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)
                                                                                 cTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                  <  2:2811591/100‑1 (MQ=255)
                                                                                  tGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTAGCACCGCTttgtt                  >  2:1197935/1‑99 (MQ=255)
                                                                                  tGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgttg                 >  2:333445/1‑100 (MQ=255)
                                                                                  tGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgttg                 >  1:451082/1‑100 (MQ=255)
                                                                                  tGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgttg                 >  2:796557/1‑100 (MQ=255)
                                                                                  tGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgttg                 <  2:57025/100‑1 (MQ=255)
                                                                                  tGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgttg                 <  1:3437383/100‑1 (MQ=255)
                                                                                   gCGAGCCATACGGCCCGCAGAAGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATATGCCTGATTTATTACCACCGCTTTGTTGa                >  2:231537/1‑100 (MQ=255)
                                                                                    cGAGCCATACGGCCCGCAGATGCAATTATCTGTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGAc               <  1:3788254/100‑1 (MQ=255)
                                                                                     gAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGAcc              <  1:1542620/100‑1 (MQ=255)
                                                                                      aGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCg             <  2:1676061/100‑1 (MQ=255)
                                                                                      aGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCg             >  1:297210/1‑100 (MQ=255)
                                                                                         cATACGGCCCGCAGATGCAATTATCTTTATCAGATGATATTAATAATTAAAGGTGAATGAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCt          >  1:3132934/1‑100 (MQ=255)
                                                                                            aCGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttt       >  1:872284/1‑100 (MQ=255)
                                                                                                 ccGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTAcc  <  1:530536/100‑1 (MQ=255)
                                                                                                 ccGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTAcc  >  1:3232695/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
CTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACC  >  NC_000913/1315710‑1315906

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