Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,758,577 G→T 100% N8K (AAC→AAA yfjH ← CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,758,5770GT98.5% 430.9 / ‑1.8 131N8K (AAC→AAAyfjHCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (65/64);  total (66/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 133 total aligned reads displayed.

CTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAATACTA  >  NC_000913/2758482‑2758673
                                                                                               |                                                                                                
cTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGtttt                                                                                               >  2:2937052/1‑99 (MQ=255)
cTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttt                                                                                              <  2:2043657/100‑1 (MQ=255)
cTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttt                                                                                              <  1:469862/100‑1 (MQ=255)
   cAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttttt                                                                                            >  2:1427633/1‑99 (MQ=255)
   cAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGGTTTTTTg                                                                                           >  1:1774311/1‑100 (MQ=255)
     aaTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGcc                                                                                         <  1:866477/100‑1 (MQ=255)
     aaTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGcc                                                                                         <  1:1703039/100‑1 (MQ=255)
     aaTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGc                                                                                          >  2:892640/1‑99 (MQ=255)
     aaTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGc                                                                                          >  1:2530309/1‑99 (MQ=255)
      aTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCt                                                                                        <  2:620794/100‑1 (MQ=255)
           ttCATAAAAATCATTCATATAAATCATAAAACTATCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAg                                                                                   <  1:182505/100‑1 (MQ=255)
           ttCATAAAAATCATTCAGAGAAATCATAAAACCTCCATAATTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAg                                                                                   <  1:529164/100‑1 (MQ=255)
           ttCATAAAAATCATTCAGAGAAATCATAAAACCACTATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAg                                                                                   <  1:686820/100‑1 (MQ=255)
           ttCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAg                                                                                   <  1:1264480/100‑1 (MQ=255)
           ttCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAg                                                                                   >  1:435298/1‑100 (MQ=255)
           ttCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAg                                                                                   <  1:2103174/100‑1 (MQ=255)
           ttCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTCGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAg                                                                                   <  2:2490934/100‑1 (MQ=255)
            tCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGa                                                                                  >  1:1091900/1‑100 (MQ=255)
            tCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGa                                                                                  >  1:2127776/1‑100 (MQ=255)
             cATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTTTCTTa                                                                                      >  2:843166/1‑95 (MQ=255)
             cATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGAt                                                                                 >  2:1663404/1‑100 (MQ=255)
               tAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATa                                                                                <  2:172818/99‑1 (MQ=255)
               tAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATAc                                                                               >  2:740795/1‑100 (MQ=255)
                  aaaTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCa                                                                            <  2:1057913/100‑1 (MQ=255)
                  aaaTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCa                                                                            >  1:2141919/1‑100 (MQ=255)
                  aaaTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATAGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCa                                                                            >  1:971795/1‑100 (MQ=255)
                    aTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCAtt                                                                          <  2:2626468/100‑1 (MQ=255)
                      cATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTaa                                                                        <  1:2424842/100‑1 (MQ=255)
                       aTTCAGAGAAATCATAAAACCACCATAACTCAATGTATTAGAGGATGTAAAAATTGTATTTTTAATACATTGTTTTTTTGCCTTAAAGATACCCATTAAt                                                                       >  1:834244/1‑100 (MQ=255)
                        ttCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTTTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtc                                                                      <  1:2956379/100‑1 (MQ=255)
                        ttCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtc                                                                      <  1:167417/100‑1 (MQ=255)
                        ttCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtc                                                                      <  2:1714845/100‑1 (MQ=255)
                         tCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtct                                                                     >  1:1780793/1‑100 (MQ=255)
                          cAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTTTTTTTTATTTACCCATTAAtct                                                                     >  2:796047/1‑99 (MQ=255)
                          cAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtctc                                                                    <  1:310112/100‑1 (MQ=255)
                          cAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtctc                                                                    >  1:2028842/1‑100 (MQ=255)
                          cAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtctc                                                                    >  1:2770490/1‑100 (MQ=255)
                          cAGAGAAATCATAAAACCACCATAAATCAATGGATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtctc                                                                    >  2:167214/1‑100 (MQ=255)
                           agagAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTTCCTTAAAGATACCCATTAAtct                                                                     >  1:2678564/1‑98 (MQ=255)
                           agagAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtctct                                                                   >  1:2816160/1‑100 (MQ=255)
                              gaAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCAATTAATctctct                                                                 >  1:2925162/1‑99 (MQ=255)
                                aaTCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTCCATGGTTTTTCTGCCTTAAAGATACCCATTAATCTCTCTgta                                                              <  1:2048329/100‑3 (MQ=255)
                                aaTCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCa                                                              <  1:1101738/100‑1 (MQ=255)
                                 aTCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCaa                                                             <  2:1881344/100‑1 (MQ=255)
                                 aTCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCaa                                                             >  1:1335902/1‑100 (MQ=255)
                                 aTCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCaa                                                             >  1:2602835/1‑100 (MQ=255)
                                  tCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAc                                                            <  2:172674/100‑1 (MQ=255)
                                  tCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAc                                                            <  1:1259378/100‑1 (MQ=255)
                                  tCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAc                                                            >  1:701561/1‑100 (MQ=255)
                                   cATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAcc                                                           >  2:1692471/1‑100 (MQ=255)
                                   cATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAcc                                                           >  1:684343/1‑100 (MQ=255)
                                   cATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAcc                                                           >  1:2031497/1‑100 (MQ=255)
                                   cATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAcc                                                           >  1:1595756/1‑100 (MQ=255)
                                     tAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCaa                                                         >  1:1548670/1‑100 (MQ=255)
                                      aaaaCCACCATCATTCAATGTATAAAAAGAAGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAcc                                                           <  2:1394279/97‑1 (MQ=255)
                                       aaaCCACCATAAATCAATGTATTATAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCaaa                                                        <  1:843170/99‑1 (MQ=255)
                                          ccaccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTg                                                     <  1:22018/99‑1 (MQ=255)
                                          ccaccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGa                                                    <  2:304165/100‑1 (MQ=255)
                                           caccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGaa                                                   >  1:2398681/1‑100 (MQ=255)
                                           caccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTATTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGaa                                                   >  1:2399971/1‑100 (MQ=255)
                                           caccaTAAATCAATGGATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGaa                                                   <  2:1462323/100‑1 (MQ=255)
                                            accaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAAc                                                  >  1:1581042/1‑100 (MQ=255)
                                             ccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAAc                                                  <  2:2094538/99‑1 (MQ=255)
                                              caTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAAc                                                  <  1:360118/98‑1 (MQ=255)
                                                tAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTTCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAAt                                              >  2:1196411/1‑100 (MQ=255)
                                                 aaaTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATg                                             <  2:1210224/100‑1 (MQ=255)
                                                 aaaTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATg                                             <  1:2240004/100‑1 (MQ=255)
                                                 aaaTCAATGTATTAGAAGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACACATTAATCTCTCTGCAACCAAAGTGAACCaa                                               >  2:705933/1‑98 (MQ=255)
                                                   aTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATgag                                           <  2:2850241/100‑1 (MQ=255)
                                                        tGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTTTTTTTATGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGca                                      >  2:2056803/1‑100 (MQ=255)
                                                        tGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGagag                                         <  2:2418548/97‑1 (MQ=255)
                                                        tGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGca                                      <  1:392382/100‑1 (MQ=255)
                                                           aTTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGcaaca                                   >  1:1141565/1‑100 (MQ=255)
                                                                 ggTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATg                             >  1:292567/1‑100 (MQ=255)
                                                                 ggTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAAGCAACAAGAATg                             >  1:608455/1‑100 (MQ=255)
                                                                  gTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGATCCAATGAGAGGCAACAAGAATGa                            <  1:2577260/100‑1 (MQ=255)
                                                                   ttGTAAAAATTGTATTTTTAGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCGATCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGa                            <  2:2591079/99‑1 (MQ=255)
                                                                    tGTAAAAATTGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAAc                          <  2:1774291/100‑1 (MQ=255)
                                                                     gTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACa                         <  1:1489953/100‑1 (MQ=255)
                                                                        aaaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatc                      >  1:3024530/1‑100 (MQ=255)
                                                                        aaaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatc                      <  1:494304/100‑1 (MQ=255)
                                                                        aaaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATAAATCTATCTCCAACCAAAGTGAACCAATGAGAGGCAACAAGAATAAACAatc                      >  2:1311475/1‑100 (MQ=255)
                                                                         aaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcg                     >  2:2407549/1‑100 (MQ=255)
                                                                            ttGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcgatc                  <  1:1238923/100‑1 (MQ=255)
                                                                            ttGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcgatc                  <  2:336731/100‑1 (MQ=255)
                                                                            tacaaTTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcgatc                  >  2:2884266/5‑100 (MQ=255)
                                                                             tGTATTTTTGGTACATGGTTTTTTTTTTTTAAAGAAACCCATTAAGCTCTCTGCAACCAAAATGAACCAATTAGAGGCAACAAGAATGAACAatcgaacg                 >  2:889049/1‑100 (MQ=255)
                                                                             tGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGTCAACAAGAATGAACAatcgatcg                 >  1:2132267/1‑100 (MQ=255)
                                                                               tATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGc               <  2:361328/100‑1 (MQ=255)
                                                                                aTTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGcc              >  1:1406863/1‑100 (MQ=255)
                                                                                 tttttGGTACATGGTTTTTTTTCCTTAAAAATACCCATTAATCTCTCTGCAACCAAAGTGAACAAATAAGAGACAACAAGAATGAACAaacgatc                  >  2:1761446/1‑95 (MQ=255)
                                                                                    ttGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTa          <  1:1638671/100‑1 (MQ=255)
                                                                                     tGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAg         <  2:499489/100‑1 (MQ=255)
                                                                                       gTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGaa       <  2:1604852/100‑1 (MQ=255)
                                                                                         ncATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAAta     >  2:2591870/2‑100 (MQ=255)
                                                                                         aCATGGTTTTTTTTCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAAta     >  2:1443268/1‑100 (MQ=255)
                                                                                           aTGGTTTTTTTGCCTTAAAGATACAAATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACCAGACTGAACAATCAAACGGCAGTTAGAAtact   >  2:609088/1‑100 (MQ=255)
                                                                                            tGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAAtacta  >  2:1554134/1‑100 (MQ=255)
                                                                                            tGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAAtacta  >  2:1554240/1‑100 (MQ=255)
                                                                                            tGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATAAGAGGCAACAAGAAAGAACAATAGATCGACCGTTAGAAtacta  >  1:106242/1‑100 (MQ=255)
                                                                                               |                                                                                                
CTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAATACTA  >  NC_000913/2758482‑2758673

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: