Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,197,739 A→G S167P (TCC→CCC)  ymfE ← e14 prophage; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,7390AG98.9% 332.4 / NA 94S167P (TCC→CCC) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (56/37);  total (56/38)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

CTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTAC  >  NC_000913/1197643‑1197835
                                                                                                |                                                                                                
ctGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAt                                                                                               >  1:306256/1‑100 (MQ=255)
 tGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAtt                                                                                              >  1:46300/1‑100 (MQ=255)
 tGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGAAAAAAAAGATTATTAGGAAGGGAAtt                                                                                              <  1:1127917/100‑1 (MQ=255)
  gagaATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAtt                                                                                              <  2:937343/99‑1 (MQ=255)
  gagaATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTa                                                                                             >  1:844316/1‑100 (MQ=255)
      aTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGc                                                                                         >  1:1169235/1‑100 (MQ=255)
      aTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGc                                                                                         <  2:2031648/100‑1 (MQ=255)
       tGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCa                                                                                        >  2:1221655/1‑100 (MQ=255)
       tGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCa                                                                                        >  1:1670563/1‑100 (MQ=255)
         aTAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAAt                                                                                      >  1:1957414/1‑100 (MQ=255)
          tAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAAtt                                                                                     <  2:323805/100‑1 (MQ=255)
          tAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAAtt                                                                                     >  2:1603447/1‑100 (MQ=255)
          tAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAAtt                                                                                     >  1:1673980/1‑100 (MQ=255)
              aTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCg                                                                                 <  1:2426/100‑1 (MQ=255)
              aTTACAAAAAAGAATAAAGACAGCAAAGAACCGTAAACATGAAAAATATCCATTTTTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCg                                                                                 <  2:640493/100‑1 (MQ=255)
                tACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaa                                                                               <  1:1221661/100‑1 (MQ=255)
                   aaaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAt                                                                            >  2:1613724/1‑100 (MQ=255)
                   aaaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAt                                                                            >  1:1845565/1‑100 (MQ=255)
                    aaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAtt                                                                           >  2:1862835/1‑100 (MQ=255)
                    aaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAtt                                                                           >  1:519012/1‑100 (MQ=255)
                     aaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCAATAAAAAAAAAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTc                                                                          <  2:1989288/100‑1 (MQ=255)
                      aaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCa                                                                         <  2:306255/100‑1 (MQ=255)
                      aaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCa                                                                         >  2:557138/1‑100 (MQ=255)
                       aaGAATAAAGACATCAAATAACCTACAACATGAAAAATATCCAGTTGTGTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaa                                                                        <  2:2035040/100‑1 (MQ=255)
                        aGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaa                                                                       <  2:1937166/100‑1 (MQ=255)
                        aGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaa                                                                       >  1:1153953/1‑100 (MQ=255)
                         gAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaa                                                                      >  1:1060788/1‑100 (MQ=255)
                            tAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAat                                                                   <  1:863340/100‑1 (MQ=255)
                            tAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAat                                                                   >  1:197546/1‑100 (MQ=255)
                             aaaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAata                                                                  >  1:377506/1‑100 (MQ=255)
                             aaaGACAGAAAAGAAACTTTAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAata                                                                  <  2:1675894/100‑1 (MQ=255)
                              aaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAatat                                                                 >  2:688195/1‑100 (MQ=255)
                              aaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAatat                                                                 >  1:887702/1‑100 (MQ=255)
                              aaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAatat                                                                 >  2:338897/1‑100 (MQ=255)
                              aaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAatat                                                                 >  2:1332853/1‑100 (MQ=255)
                                gACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATcc                                                               <  2:1599157/100‑1 (MQ=255)
                                gACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATcc                                                               <  2:1819025/100‑1 (MQ=255)
                                gACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATcc                                                               <  1:239145/100‑1 (MQ=255)
                                  cAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaa                                                             <  1:192963/100‑1 (MQ=255)
                                    gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                           >  2:1613493/1‑100 (MQ=255)
                                    gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                           >  2:893426/1‑100 (MQ=255)
                                    gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                           >  2:292233/1‑100 (MQ=255)
                                    gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                           >  2:650055/1‑100 (MQ=255)
                                    gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                           >  2:162442/1‑100 (MQ=255)
                                    gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                           >  1:920007/1‑100 (MQ=255)
                                       aaGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATa                                                        >  2:1532948/1‑100 (MQ=255)
                                          aaCCTTAAACATGTAAAATATCCATTTAATTGAAAAAAAAGATTATTAGGAAAGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAg                                                       <  2:1687301/98‑1 (MQ=255)
                                             cTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATAc                                                  >  2:1048407/1‑100 (MQ=255)
                                                    aTGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTc                                           >  2:116456/1‑100 (MQ=255)
                                                     tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTcc                                          >  1:326405/1‑100 (MQ=255)
                                                     tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTcc                                          <  1:807478/100‑1 (MQ=255)
                                                      gAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCa                                         >  1:113924/1‑100 (MQ=255)
                                                       aaaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCag                                        >  2:488094/1‑100 (MQ=255)
                                                       aaaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCag                                        <  1:688197/100‑1 (MQ=255)
                                                       aaaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCag                                        >  1:722515/1‑100 (MQ=255)
                                                          aaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCagaag                                     <  1:650057/100‑1 (MQ=255)
                                                           atatCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCagaaga                                    >  1:1428599/1‑100 (MQ=255)
                                                             atCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGt                                  >  2:1063037/1‑100 (MQ=255)
                                                              tCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGtt                                 <  1:1048412/100‑1 (MQ=255)
                                                               ccATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTc                                <  2:1525896/100‑1 (MQ=255)
                                                               ccATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTc                                <  1:1603459/100‑1 (MQ=255)
                                                                cATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCa                               >  2:1353860/1‑100 (MQ=255)
                                                                 aTTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCaa                              <  1:1603618/100‑1 (MQ=255)
                                                                  tttgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAat                             <  1:84287/100‑1 (MQ=255)
                                                                    tgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAatat                           <  2:113923/100‑1 (MQ=255)
                                                                    tgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAatat                           >  1:1668348/1‑100 (MQ=255)
                                                                    tggttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCCat                             >  2:280724/4‑98 (MQ=255)
                                                                     gtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAAtata                          >  1:1068474/1‑100 (MQ=255)
                                                                     gtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAAtata                          >  1:201635/1‑100 (MQ=255)
                                                                     gtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAAtata                          >  1:1934442/1‑100 (MQ=255)
                                                                        tgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATg                       <  1:2055/100‑1 (MQ=255)
                                                                        tgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATg                       <  1:1691938/100‑1 (MQ=255)
                                                                        tacAAAAAAAGATTATTAGGAAGCGAAGAAATTAAATTATCGAAAATTCAAAAAATATCCAAAAACATTATACTTTATCTCAGAAGAGTTCAAtata                          <  2:880725/95‑1 (MQ=255)
                                                                         gCAAAAAAAGATTATTAGGAAGGGAATTAATGCCATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                      >  2:1980745/1‑100 (MQ=255)
                                                                         gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                      >  1:240514/1‑100 (MQ=255)
                                                                         gCAAAAAAAGATTATTAGGAAGGGAATTAATACAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                      >  2:235237/1‑100 (MQ=255)
                                                                          cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                     >  2:1807204/1‑100 (MQ=255)
                                                                            aaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                   <  1:292234/100‑1 (MQ=255)
                                                                            aaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                   <  1:526153/100‑1 (MQ=255)
                                                                            aaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                   <  2:1169229/100‑1 (MQ=255)
                                                                              aaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTc                 <  1:488095/100‑1 (MQ=255)
                                                                              aaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAATAAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTc                 <  2:1934424/100‑1 (MQ=255)
                                                                                    ttattaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCagt           >  2:1243786/1‑98 (MQ=255)
                                                                                     tattaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAtt          >  2:1654977/1‑100 (MQ=255)
                                                                                      attaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttt         >  2:1808753/1‑100 (MQ=255)
                                                                                      attaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttt         >  1:1740072/1‑100 (MQ=255)
                                                                                       ttaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAtttt        >  1:106094/1‑100 (MQ=255)
                                                                                       ttaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAtttt        >  1:2005562/1‑100 (MQ=255)
                                                                                       ttaGGAAGGGAATTAATGCAATTATCAAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAtttt        <  2:844313/100‑1 (MQ=255)
                                                                                         aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc      <  2:377505/100‑1 (MQ=255)
                                                                                         aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc      <  1:816583/100‑1 (MQ=255)
                                                                                         aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc      <  2:1673967/100‑1 (MQ=255)
                                                                                         aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc      <  2:1957396/100‑1 (MQ=255)
                                                                                          ggaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCt     <  1:961913/100‑1 (MQ=255)
                                                                                          ggaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCt     >  2:761682/1‑100 (MQ=255)
                                                                                           gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt    <  1:1793663/100‑1 (MQ=255)
                                                                                           gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt    <  2:354832/100‑1 (MQ=255)
                                                                                           gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt    <  1:772711/100‑1 (MQ=255)
                                                                                           gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt    <  1:162443/100‑1 (MQ=255)
                                                                                             agggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTAc  <  1:971718/100‑1 (MQ=255)
                                                                                                |                                                                                                
CTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTAC  >  NC_000913/1197643‑1197835

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