Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 631,241 A→G intergenic (+19/‑164) cstA → / → ybdD carbon starvation protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913631,2410AG99.3% 428.1 / ‑3.7 139intergenic (+19/‑164)cstA/ybdDcarbon starvation protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (74/64);  minor base . (0/1);  total (74/65)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.52e-01

Only 100 of 148 total aligned reads displayed.

AAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGAT  >  NC_000913/631151‑631330
                                                                                          |                                                                                         
aaagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTg                                                                                         >  1:4347559/1‑93 (MQ=255)
aaagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTg                                                                                         >  2:4013079/1‑93 (MQ=255)
 aagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGgtgt                                                                                        >  1:1791072/1‑93 (MQ=255)
   gaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTg                                                                                      <  2:826512/93‑1 (MQ=255)
   gaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTg                                                                                      <  2:1543531/93‑1 (MQ=255)
   gaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTg                                                                                      >  1:172635/1‑93 (MQ=255)
    aaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTgg                                                                                     >  1:3096377/1‑93 (MQ=255)
    aaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTgg                                                                                     <  1:3627206/93‑1 (MQ=255)
    aaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTgg                                                                                     >  2:1875402/1‑93 (MQ=255)
     aaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGc                                                                                    <  1:1322791/93‑1 (MQ=255)
        gCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGaa                                                                                 <  1:3364087/93‑1 (MQ=255)
          cGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATg                                                                               <  1:3101162/93‑1 (MQ=255)
          cGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGGGTTGGCGAATg                                                                               >  2:1632640/1‑93 (MQ=255)
            tATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGtt                                                                             <  2:3376939/93‑1 (MQ=255)
            tATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGtt                                                                             >  2:734111/1‑93 (MQ=255)
            tATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGtg                                                                             >  2:1018104/1‑92 (MQ=255)
            tATGAACCAATGCCGGAAAAGCCCTAGGAGAGCGCGGCGCAGGCAAAAGCCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGCATGtt                                                                             <  2:4077588/93‑1 (MQ=255)
             aTGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTg                                                                            >  2:1482041/1‑93 (MQ=255)
             aTGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTg                                                                            >  2:989906/1‑93 (MQ=255)
              tGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGGATGTTgg                                                                           >  1:3479882/1‑93 (MQ=255)
               gAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGc                                                                          <  1:1368237/93‑1 (MQ=255)
                aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCt                                                                         >  1:2060589/1‑93 (MQ=255)
                  ccAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGGTTa                                                                       >  1:1812161/1‑93 (MQ=255)
                  ccAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTa                                                                       >  1:1846458/1‑93 (MQ=255)
                  ccAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTa                                                                       <  1:2320845/93‑1 (MQ=255)
                   cAATGCGGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAg                                                                      >  2:3902296/1‑93 (MQ=255)
                   cAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGGTTat                                                                      >  1:3944606/1‑92 (MQ=255)
                      tGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGc                                                                   >  2:4000721/1‑93 (MQ=255)
                       gCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGcc                                                                  >  2:1950536/1‑93 (MQ=255)
                       gCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGGGTTGGCGAATGTTGGCTTAGTGcc                                                                  >  2:2562403/1‑93 (MQ=255)
                        ccGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGccc                                                                 <  1:3110159/93‑1 (MQ=255)
                          ggAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAg                                                               >  2:361602/1‑93 (MQ=255)
                          ggAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGCGTTGGCGAATGTTGGCTTAGTGCCGAg                                                               >  2:1340172/1‑93 (MQ=255)
                           gAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAgg                                                              >  2:135100/1‑93 (MQ=255)
                           gAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAgg                                                              >  1:3975950/1‑93 (MQ=255)
                            aaaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAggg                                                             >  1:1903382/1‑93 (MQ=255)
                            aaaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAggg                                                             <  1:3819508/93‑1 (MQ=255)
                            aaaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAggg                                                             <  1:3640798/93‑1 (MQ=255)
                             aaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGt                                                            >  2:1618922/1‑93 (MQ=255)
                              aaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGCGAGGGGGGGGAGTTGGCGAATGTTGGCTTAGGGCCCAggggt                                                           >  2:1030845/1‑91 (MQ=255)
                               aTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTc                                                          >  1:1351052/1‑93 (MQ=255)
                                  tCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCct                                                       >  1:4610321/1‑93 (MQ=255)
                                   cGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctc                                                      >  2:391752/1‑93 (MQ=255)
                                    gaggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctct                                                     >  2:4054322/1‑93 (MQ=255)
                                     aggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCtctc                                                    >  2:3708103/1‑93 (MQ=255)
                                     aggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCtctc                                                    <  2:1509845/93‑1 (MQ=255)
                                      ggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCtctcc                                                   >  2:1787699/1‑92 (MQ=255)
                                      ggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGGGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCa                                                   >  2:597157/1‑93 (MQ=255)
                                         gaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAccc                                                >  2:3203985/1‑93 (MQ=255)
                                         gaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAccc                                                >  1:4263637/1‑93 (MQ=255)
                                         gaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAccc                                                >  1:2278933/1‑93 (MQ=255)
                                          aTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCt                                               >  1:1196215/1‑93 (MQ=255)
                                          aTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGCTCCCTCTCCCCCt                                               >  2:2612946/1‑93 (MQ=255)
                                             gTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAc                                            <  2:2811359/93‑1 (MQ=255)
                                               ggCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCct                                         <  2:2121586/93‑1 (MQ=255)
                                                gcgcAGGCAAAATGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCct                                         >  1:642967/1‑93 (MQ=255)
                                                gcgcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCct                                         >  1:1770933/1‑93 (MQ=255)
                                                 cgcAGGCAAAAGTCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGCCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCctc                                        <  1:2084968/93‑1 (MQ=255)
                                                     ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTccc                                    <  2:4181103/93‑1 (MQ=255)
                                                     ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTccc                                    >  1:2554275/1‑93 (MQ=255)
                                                      gCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTcccc                                   <  2:2903495/93‑1 (MQ=255)
                                                       cAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCg                                  <  2:251979/93‑1 (MQ=255)
                                                       cAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCg                                  >  2:4392205/1‑93 (MQ=255)
                                                        aaaaGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCgg                                 <  2:473506/93‑1 (MQ=255)
                                                         aaaGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGCTGGCGAATGTTGGCTTAGTGTCCAGGGTTCCCTCTCACCCTAACCCTCTCCTCggg                                >  1:2429013/1‑92 (MQ=255)
                                                         aaaGGCGCACACTAAAGTCAGAGTGAGGGGGGGGCGTTGGCGAATGTTGGCTTATTGCCCAGGGTTACCTCTCACCCTAACCCTCTCCCCgcg                                >  2:3675432/1‑91 (MQ=255)
                                                           aGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                              <  1:1652382/93‑1 (MQ=255)
                                                           aGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                              <  1:4003889/93‑1 (MQ=255)
                                                           aGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                              <  2:3701989/93‑1 (MQ=255)
                                                              cgcACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGc                           <  2:4669162/93‑1 (MQ=255)
                                                               gcACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCg                          >  2:734899/1‑93 (MQ=255)
                                                                cacaCTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGa                         <  2:4138785/93‑1 (MQ=255)
                                                                 acacTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAg                        >  1:2169715/1‑93 (MQ=255)
                                                                 acacTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAg                        >  2:2839822/1‑93 (MQ=255)
                                                                  cacTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg                       <  1:1637127/93‑1 (MQ=255)
                                                                  cacTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg                       <  2:1206457/93‑1 (MQ=255)
                                                                   acTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAggg                      <  2:2661743/93‑1 (MQ=255)
                                                                    cTAAAGTCAGAGCGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgggg                     >  2:3329038/1‑93 (MQ=255)
                                                                     tAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGa                    <  1:2828580/93‑1 (MQ=255)
                                                                     tAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGa                    >  1:4474838/1‑93 (MQ=255)
                                                                     tAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGa                    >  1:4665742/1‑93 (MQ=255)
                                                                      aaaGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCCCCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc                   <  2:3500038/93‑1 (MQ=255)
                                                                      aaaGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc                   <  1:3583527/93‑1 (MQ=255)
                                                                      aaaGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc                   <  1:294308/93‑1 (MQ=255)
                                                                       aaGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACt                  <  1:2533776/93‑1 (MQ=255)
                                                                       aaGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACt                  >  1:653370/1‑93 (MQ=255)
                                                                       aaGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTACCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACt                  <  1:401350/93‑1 (MQ=255)
                                                                         gTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGa                <  1:1329188/93‑1 (MQ=255)
                                                                          tCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAc               <  2:3197228/93‑1 (MQ=255)
                                                                          tCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAc               <  2:2694551/93‑1 (MQ=255)
                                                                          tCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAc               <  1:4488737/93‑1 (MQ=255)
                                                                          gcAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAc               <  1:356268/92‑1 (MQ=255)
                                                                             gagTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGa            <  2:4635295/93‑1 (MQ=255)
                                                                              agTGAGGGGGGGGTGTTTGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAg           <  2:2230897/93‑1 (MQ=255)
                                                                              agTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAg           <  2:1539660/93‑1 (MQ=255)
                                                                                    ggggggGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGtt     >  2:4426014/1‑93 (MQ=255)
                                                                                    ggggggGTGTTGGCGAATGTTGGCTGAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGtt     <  1:2615110/93‑1 (MQ=255)
                                                                                      ggggGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGa   <  2:2979718/93‑1 (MQ=255)
                                                                                       ggggTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGAt  >  1:2955972/1‑93 (MQ=255)
                                                                                       ggggTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGAt  >  1:3802678/1‑93 (MQ=255)
                                                                                          |                                                                                         
AAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGAT  >  NC_000913/631151‑631330

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