Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 632,017 (G)7→8 intergenic (+18/‑165) cstA → / → ybdD carbon starvation protein A/conserved protein YbdD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913632,0101.G98.4% 228.8 / ‑3.2 63intergenic (+11/‑172)cstA/ybdDcarbon starvation protein A/conserved protein YbdD
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (25/37);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-01

GCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCAT  >  NC_000913/631871‑632143
                                                                                                                                            |                                                                                                                                     
gCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGa                                                                                                                                        >  1:846633/1‑140 (MQ=255)
     gTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAggggg                                                                                                                                  >  1:1196783/1‑140 (MQ=255)
     gTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAggggg                                                                                                                                  >  1:120533/1‑140 (MQ=255)
      ttCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAgggggg                                                                                                                                 >  1:937108/1‑140 (MQ=255)
      ttCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAgggggg                                                                                                                                 >  1:89424/1‑140 (MQ=255)
      ttCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAgggggg                                                                                                                                 >  1:206631/1‑140 (MQ=255)
          cGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAggggggggat                                                                                                                             <  2:1277224/140‑4 (MQ=255)
                aGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGc                                                                                                                        >  2:959685/1‑140 (MQ=255)
                aGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGc                                                                                                                        >  2:957143/1‑140 (MQ=255)
                aGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGc                                                                                                                        >  1:147736/1‑140 (MQ=255)
                aGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGCGAGGGGGGGGATGTTGGc                                                                                                                        <  2:883053/140‑1 (MQ=255)
                 gACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCg                                                                                                                       <  1:183054/140‑1 (MQ=255)
                     gCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATg                                                                                                                   >  2:448421/1‑140 (MQ=255)
                      cACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGt                                                                                                                  <  1:957143/140‑1 (MQ=255)
                        cTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTg                                                                                                                <  1:270217/140‑1 (MQ=255)
                           gCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCt                                                                                                             <  2:1144517/140‑1 (MQ=255)
                           gCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCtt                                                                                                            <  2:120533/140‑1 (MQ=255)
                             ggCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTa                                                                                                           <  2:793064/140‑1 (MQ=255)
                             ggCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTa                                                                                                           <  2:147736/140‑1 (MQ=255)
                               cATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGt                                                                                                         >  2:807824/1‑140 (MQ=255)
                                 ttGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGc                                                                                                       <  1:1303843/140‑1 (MQ=255)
                                   gAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGccc                                                                                                     >  2:1106209/1‑140 (MQ=255)
                                   gAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGccc                                                                                                     <  1:379615/140‑1 (MQ=255)
                                    aaaGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCa                                                                                                    <  2:206631/140‑1 (MQ=255)
                                        aTCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGt                                                                                                <  2:937108/140‑1 (MQ=255)
                                         tCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGtt                                                                                               <  1:448421/140‑1 (MQ=255)
                                          ccGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTc                                                                                              <  2:513297/140‑1 (MQ=255)
                                             aaGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCct                                                                                           <  2:846633/140‑1 (MQ=255)
                                              aGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctc                                                                                          >  2:140955/1‑140 (MQ=255)
                                                 cAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCa                                                                                       <  2:157354/140‑1 (MQ=255)
                                                 cAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCa                                                                                       <  2:13019/140‑1 (MQ=255)
                                                  aaCGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTAGCTTAGTGCCCAGGGTTCCCTCTCAc                                                                                      >  1:601140/1‑140 (MQ=255)
                                                   aCGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAcc                                                                                     >  1:951677/1‑140 (MQ=255)
                                                   aCGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAcc                                                                                     >  1:1270084/1‑140 (MQ=255)
                                                    cGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAg                                                                                                   <  2:549281/125‑1 (MQ=255)
                                                    cGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAg                                                                                                   >  1:549281/1‑125 (MQ=255)
                                                    cGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGTTTCCCTCTCAccc                                                                                    >  1:1090822/1‑140 (MQ=255)
                                                         aaagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGTTTCCCTCTCACCCTAAcc                                                                               <  2:1090822/140‑1 (MQ=255)
                                                         aaagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAcc                                                                               <  2:357285/140‑1 (MQ=255)
                                                           agaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCct                                                                             <  2:1081044/140‑1 (MQ=255)
                                                                cGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTccc                                                                        <  1:959685/140‑1 (MQ=255)
                                                                      aTGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCTAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                                                                  <  1:262337/140‑1 (MQ=255)
                                                                      aTGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                                                                  <  2:623188/140‑1 (MQ=255)
                                                                      aTGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                                                                  <  2:1270084/140‑1 (MQ=255)
                                                                          aCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCg                                                              <  2:918612/140‑1 (MQ=255)
                                                                               tGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgggg                                                         >  2:87406/1‑140 (MQ=255)
                                                                                gCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGa                                                        <  2:1190172/140‑1 (MQ=255)
                                                                                    gAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGa                                                    >  1:649003/1‑140 (MQ=255)
                                                                                    gAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGa                                                    >  2:618696/1‑140 (MQ=255)
                                                                                        aTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGa                                                <  1:968955/140‑1 (MQ=255)
                                                                                         tGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAg                                               >  1:966522/1‑140 (MQ=255)
                                                                                           tCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAgcg                                             <  1:618696/140‑1 (MQ=255)
                                                                                             gaggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGcgcg                                           <  2:612576/140‑1 (MQ=255)
                                                                                              aggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGt                                          >  2:1034078/1‑140 (MQ=255)
                                                                                                  gaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGAt                                      <  2:172278/140‑1 (MQ=255)
                                                                                                   aTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATa                                     >  2:1418969/1‑140 (MQ=255)
                                                                                                    tCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAg                                    <  2:649003/140‑1 (MQ=255)
                                                                                                    tCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAg                                    <  2:1196783/140‑1 (MQ=255)
                                                                                                    tCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAg                                    <  1:1200403/140‑1 (MQ=255)
                                                                                                    tCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAg                                    <  2:1333703/140‑1 (MQ=255)
                                                                                                        ggCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCAtt                                >  1:1299529/1‑140 (MQ=255)
                                                                                                          cacAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTg                              >  2:422063/3‑140 (MQ=255)
                                                                                                              ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAg                                    <  2:282484/130‑1 (MQ=255)
                                                                                                              ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAg                                    >  1:282484/1‑130 (MQ=255)
                                                                                                              ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGACGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAgg                          >  2:1132420/1‑140 (MQ=255)
                                                                                                               gCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGc                         <  1:314326/140‑1 (MQ=255)
                                                                                                                         cacaCTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGc               >  2:566413/1‑140 (MQ=255)
                                                                                                                         cacaCTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGc               >  2:432477/1‑140 (MQ=255)
                                                                                                                          acacTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCg              <  2:1317074/140‑1 (MQ=255)
                                                                                                                                      gagTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCAt  >  2:237411/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                     
GCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCAT  >  NC_000913/631871‑632143

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