Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,758,577 G→T 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,5770GT93.4% 352.1 / 15.6 106N8K (AAC→AAA) yfjHCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (2/5);  new base T (50/49);  total (52/54)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

Only 100 of 109 total aligned reads displayed.

CTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAATACTA  >  NC_000913/2758482‑2758673
                                                                                               |                                                                                                
cTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttt                                                                                              >  2:814347/1‑100 (MQ=255)
cTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttt                                                                                              >  1:1916221/1‑100 (MQ=255)
cTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttt                                                                                              >  1:1788681/1‑100 (MQ=255)
cTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttt                                                                                              >  1:1480596/1‑100 (MQ=255)
 tACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGtttttt                                                                                             <  1:1007822/100‑1 (MQ=255)
 tACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGtttttt                                                                                             <  2:1906696/100‑1 (MQ=255)
 tACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGtttttt                                                                                             <  2:1713661/100‑1 (MQ=255)
  aCAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttttt                                                                                            <  1:1115311/100‑1 (MQ=255)
  aCAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttttt                                                                                            <  1:1449147/100‑1 (MQ=255)
   cAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTg                                                                                           <  2:921445/100‑1 (MQ=255)
   cAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTg                                                                                           <  2:1896376/100‑1 (MQ=255)
   cAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTg                                                                                           >  2:1002943/1‑100 (MQ=255)
   cAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTg                                                                                           >  1:380005/1‑100 (MQ=255)
   cAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTg                                                                                           <  1:1885264/100‑1 (MQ=255)
    aaaTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGc                                                                                          <  1:201956/100‑1 (MQ=255)
     aaTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGcc                                                                                         <  1:655865/100‑1 (MQ=255)
      aTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCt                                                                                        <  1:1255242/100‑1 (MQ=255)
      aTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCt                                                                                        <  2:50265/100‑1 (MQ=255)
      aTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCt                                                                                        <  1:1854307/100‑1 (MQ=255)
      aTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCt                                                                                        <  2:426275/100‑1 (MQ=255)
        ttGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTa                                                                                      <  2:1295864/100‑1 (MQ=255)
          gTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTaaa                                                                                    >  2:477612/1‑100 (MQ=255)
          gTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTaaa                                                                                    >  1:1853187/1‑100 (MQ=255)
           ttCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAg                                                                                   <  2:118221/100‑1 (MQ=255)
           ttCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAg                                                                                   <  2:204742/100‑1 (MQ=255)
            tCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGa                                                                                  >  2:61856/1‑100 (MQ=255)
             cATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGAt                                                                                 >  2:71218/1‑100 (MQ=255)
              aTAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATa                                                                                >  2:1939715/1‑100 (MQ=255)
               tAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATAc                                                                               >  1:152648/1‑100 (MQ=255)
               tAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATAc                                                                               >  2:1537903/1‑100 (MQ=255)
               tAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATAc                                                                               <  2:944665/100‑1 (MQ=255)
                aaaaaTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATAc                                                                               >  1:1265542/1‑99 (MQ=255)
                  aaaTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCa                                                                            <  2:1521871/100‑1 (MQ=255)
                  aaaTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCa                                                                            <  2:1265546/100‑1 (MQ=255)
                     tCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTa                                                                         >  1:995648/1‑100 (MQ=255)
                       aTTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAt                                                                       >  1:190595/1‑100 (MQ=255)
                       aTTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAt                                                                       >  2:1596620/1‑100 (MQ=255)
                        ttCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtc                                                                      <  1:1835227/100‑1 (MQ=255)
                        ttCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtc                                                                      <  1:1760907/100‑1 (MQ=255)
                          cAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCTTTAAAGCTACCCATCAAtcgc                                                                    >  1:707577/1‑98 (MQ=255)
                          cAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtctc                                                                    <  2:1548076/100‑1 (MQ=255)
                             agaAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATctctct                                                                 <  1:1240896/100‑1 (MQ=255)
                             agaAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATctctct                                                                 >  1:1476814/1‑100 (MQ=255)
                              gaAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTg                                                                <  1:536265/100‑1 (MQ=255)
                               aaaTCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtctctc                                                                  <  2:599893/97‑1 (MQ=255)
                                aaTCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCa                                                              <  1:1030248/100‑1 (MQ=255)
                                aaTCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCa                                                              <  1:1002942/100‑1 (MQ=255)
                                 aTCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCaa                                                             >  1:591726/1‑100 (MQ=255)
                                  tCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAATAGTGTATTTTTGGTACATGGGTTTTTTGCCTTATAGATACCCATTAATCTCTCTGAAAc                                                            <  2:647294/100‑1 (MQ=255)
                                  tCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAc                                                            <  1:862993/100‑1 (MQ=255)
                                  tCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAc                                                            <  2:241425/100‑1 (MQ=255)
                                  tCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAc                                                            <  2:500232/100‑1 (MQ=255)
                                   cATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAcc                                                           <  1:201640/100‑1 (MQ=255)
                                   cATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAcc                                                           <  1:1596609/100‑1 (MQ=255)
                                   cATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAcc                                                           <  1:448132/100‑1 (MQ=255)
                                   cATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAcc                                                           >  1:1805684/1‑100 (MQ=255)
                                   cATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAc                                                            <  2:563217/99‑1 (MQ=255)
                                    aTAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTAGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATctctct                                                                 <  1:1483668/93‑1 (MQ=255)
                                       aaaCCACCATAAATCAATGTATTAGAGGTTGTAATAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTTTGCAACCAAAg                                                       <  1:1025613/100‑1 (MQ=255)
                                       aaaCCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAg                                                       >  1:69884/1‑100 (MQ=255)
                                       aaaCCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAg                                                       <  2:867589/100‑1 (MQ=255)
                                           caccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGaa                                                   >  1:1754414/1‑100 (MQ=255)
                                            accaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAAc                                                  <  1:134389/100‑1 (MQ=255)
                                              caTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTTTTTTTTTTATACCCATTAATCTCTCTGCAAACAAAGTGTAcaa                                                >  1:252289/1‑98 (MQ=255)
                                              caTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGTCTTATAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCa                                                >  1:1795838/1‑100 (MQ=255)
                                                tAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAAt                                              >  2:500457/1‑100 (MQ=255)
                                                  aaTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATga                                            >  1:163526/1‑100 (MQ=255)
                                                   aTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATgag                                           <  1:677287/100‑1 (MQ=255)
                                                    tCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATgaga                                          <  2:410475/100‑1 (MQ=255)
                                                    tCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATgaga                                          >  2:418799/1‑100 (MQ=255)
                                                     cAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGagag                                         >  2:772728/1‑100 (MQ=255)
                                                     cAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGagag                                         <  1:1939702/100‑1 (MQ=255)
                                                       acgTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGc                                       <  1:1182109/98‑1 (MQ=255)
                                                        tGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGca                                      <  1:77356/100‑1 (MQ=255)
                                                           aTTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTACCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGcaaca                                   >  1:1524798/1‑100 (MQ=255)
                                                             tAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGcaacaa                                  >  1:840671/1‑99 (MQ=255)
                                                             tAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAg                                 >  2:1717976/1‑100 (MQ=255)
                                                              agagGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGa                                >  1:829434/1‑100 (MQ=255)
                                                               gagGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTACCTTAAAGATACCCATTAATCTCTCTGCAATCAATGTGAACCAATGAGAGGCAACAATaa                               >  1:1717040/1‑100 (MQ=255)
                                                                agGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAAt                              <  1:1914758/100‑1 (MQ=255)
                                                                 ggTTGTAAAAATTGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATg                             <  2:1395979/100‑1 (MQ=255)
                                                                   ttGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGCTACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGAAACAAGGATGaa                           >  1:267859/1‑100 (MQ=255)
                                                                    tGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAAc                          <  1:506464/100‑1 (MQ=255)
                                                                    tGTAAAAATTGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAAc                          >  1:805171/1‑100 (MQ=255)
                                                                     gTAAAAATTGTATTTTTGGTACATGGTTTTTTTTTTTTATATATACCCATTAATCTCTCTGCAACTAAAGTGAACCAATGAGAGGCAACAAGAATGaa                           >  1:1614346/1‑98 (MQ=255)
                                                                      tAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACaa                        >  1:1464826/1‑100 (MQ=255)
                                                                       aaaaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAat                       >  1:3286/1‑100 (MQ=255)
                                                                       aaaaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAat                       <  1:639238/100‑1 (MQ=255)
                                                                       aaaaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAat                       <  2:1447701/100‑1 (MQ=255)
                                                                        aaaaTTGTATTTTTGGTACATGGTTTTTTTTTTTTTTTTATACCCATTAATCTCTCTGCAACCAAAGTGAACTAATGAGAGGCAACAAGAATGAACAatc                      >  1:1857005/1‑100 (MQ=255)
                                                                        aaaaTTGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatc                      <  1:1023005/100‑1 (MQ=255)
                                                                         aaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcg                     >  1:1861481/1‑100 (MQ=255)
                                                                            ttGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcgatc                  <  2:1045616/100‑1 (MQ=255)
                                                                            ttGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcgatc                  >  2:70519/1‑100 (MQ=255)
                                                                                aTTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGcc              <  2:172610/100‑1 (MQ=255)
                                                                                   tttGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGtt           >  2:1678476/1‑100 (MQ=255)
                                                                                         aCATGGTTTTTTTGTCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCCACAAGAATGAACAATCGATCGGCCGTTAGtata     >  1:1293408/1‑100 (MQ=255)
                                                                                         aCATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAAta     >  1:504687/1‑100 (MQ=255)
                                                                                          cATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAAtac    >  2:1736564/1‑100 (MQ=255)
                                                                                            tGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAAtacta  >  2:485422/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 < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑