Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,197,739 A→G 100% 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,7390AG97.7% 527.0 / ‑1.3 164S167P (TCC→CCC) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/4);  new base G (87/72);  total (88/76)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

CACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTA  >  NC_000913/1197641‑1197834
                                                                                                  |                                                                                               
cactGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAggaaggaa                                                                                                >  2:3797970/1‑100 (MQ=255)
   tgcgaATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGATTGCAAAAAAAGATTATTAGGAAGGGAAtt                                                                                             >  2:3978158/4‑100 (MQ=255)
       aaTGAAAAAATTACAAAAAAGAATAAAGACAGCAAAGAAAAAAAAAAAAAAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAAt                                                                                          <  2:1894062/99‑1 (MQ=255)
        aTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATACATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGc                                                                                        >  2:2252480/1‑100 (MQ=255)
         tGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCa                                                                                       >  1:3237142/1‑100 (MQ=255)
         tGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACACGAAAAATATCCATTTGTTTGCAAAAAAAGATTATAAGGAAGGGAATTAATGCa                                                                                       >  2:2310327/1‑100 (MQ=255)
            tAAAATTACAAAGAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAAtt                                                                                    >  1:2742029/1‑100 (MQ=255)
             aaaaTAACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTa                                                                                   <  1:1961673/100‑1 (MQ=255)
                  tACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaa                                                                              <  1:428158/100‑1 (MQ=255)
                  tACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaa                                                                              >  2:568002/1‑100 (MQ=255)
                   aCAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaaa                                                                             <  2:3710043/100‑1 (MQ=255)
                     aaaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAt                                                                           <  1:1382263/100‑1 (MQ=255)
                      aaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAACAtt                                                                          >  1:1837115/1‑100 (MQ=255)
                      aaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAtt                                                                          >  2:1073914/1‑100 (MQ=255)
                       aaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTc                                                                         <  1:556577/100‑1 (MQ=255)
                       aaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTc                                                                         >  1:2526432/1‑100 (MQ=255)
                       aaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTc                                                                         >  1:3468903/1‑100 (MQ=255)
                        aaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCa                                                                        <  1:3346222/100‑1 (MQ=255)
                          aGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATCTGTTTGCAAAAAAAGATTATTAGAAAGGGAATTAATGCAATTATCGAAAATTCaaa                                                                      >  2:3566396/1‑100 (MQ=255)
                           gAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaa                                                                     >  1:1037353/1‑100 (MQ=255)
                            aaTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaaa                                                                    >  1:2257156/1‑100 (MQ=255)
                             aaaaaaACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaaaa                                                                   <  1:3043498/98‑1 (MQ=255)
                               aaaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTTTTAGGAAGGGTATTAATGCAATTATCGAGAATTCAAAAAata                                                                 <  1:2816599/100‑1 (MQ=255)
                               aaaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAata                                                                 <  2:2889039/100‑1 (MQ=255)
                                aaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAatat                                                                >  2:1170460/1‑100 (MQ=255)
                                aaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAatat                                                                >  1:2180433/1‑100 (MQ=255)
                                 aGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATc                                                               >  2:4041756/1‑100 (MQ=255)
                                 aGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATc                                                               >  1:997954/1‑100 (MQ=255)
                                  gACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATcc                                                              <  1:1626462/100‑1 (MQ=255)
                                   aCAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCa                                                             >  1:1324332/1‑100 (MQ=255)
                                    cAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAatat                                                                >  2:3012733/1‑96 (MQ=255)
                                      gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                          >  2:2542283/1‑100 (MQ=255)
                                      gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATCAAGAAGAGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                          >  2:3104351/1‑100 (MQ=255)
                                       cAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaaa                                                         >  1:633805/1‑100 (MQ=255)
                                        aaaGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAAt                                                        <  1:2586175/100‑1 (MQ=255)
                                         aaGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATa                                                       <  1:3145524/100‑1 (MQ=255)
                                         aaGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATa                                                       >  1:2582435/1‑100 (MQ=255)
                                          aGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAg                                                      <  1:2252494/100‑1 (MQ=255)
                                             aCCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGtat                                                   >  2:1364553/1‑100 (MQ=255)
                                             aCCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGta                                                    >  2:3537622/1‑99 (MQ=255)
                                              ccTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGtata                                                  >  1:342642/1‑100 (MQ=255)
                                              ccTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGtata                                                  >  1:1075705/1‑100 (MQ=255)
                                                ttCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACt                                                <  1:3658942/100‑1 (MQ=255)
                                                    aCATGAAAAATATCCATTTGTTTGAAAAAAAGGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTAt                                            <  1:1299128/100‑1 (MQ=255)
                                                      aTGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTc                                          <  1:885902/100‑1 (MQ=255)
                                                       tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTcc                                         <  2:1549669/100‑1 (MQ=255)
                                                       tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTcc                                         <  1:3782021/100‑1 (MQ=255)
                                                       tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTcc                                         <  1:1766697/100‑1 (MQ=255)
                                                         aaaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCag                                       >  2:3184351/1‑100 (MQ=255)
                                                         aaaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCag                                       <  1:18343/100‑1 (MQ=255)
                                                         aaaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCag                                       <  2:3643568/100‑1 (MQ=255)
                                                          aaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCaga                                      >  2:1447245/1‑100 (MQ=255)
                                                            aaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCagaag                                    <  2:2998662/100‑1 (MQ=255)
                                                                tCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGtt                                >  1:3149994/1‑100 (MQ=255)
                                                                 ccATTTAATTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTAAAAAAAAAAAAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTc                               <  2:450383/100‑1 (MQ=255)
                                                                  cATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCa                              >  1:89978/1‑100 (MQ=255)
                                                                  cATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCa                              >  2:2317341/1‑100 (MQ=255)
                                                                  cATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCa                              <  2:1739654/100‑1 (MQ=255)
                                                                  cATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCa                              <  1:3169932/100‑1 (MQ=255)
                                                                  cATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATCATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCa                              >  2:3614591/1‑100 (MQ=255)
                                                                   aTTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCaa                             <  1:3075563/100‑1 (MQ=255)
                                                                   aTTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCaa                             <  2:1231196/100‑1 (MQ=255)
                                                                      tgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAatat                          >  1:3086155/1‑100 (MQ=255)
                                                                      tgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAatat                          >  1:4038453/1‑100 (MQ=255)
                                                                      tgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAatat                          >  1:3804618/1‑100 (MQ=255)
                                                                       gtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAAtata                         >  2:709149/1‑100 (MQ=255)
                                                                       aaaaaaaaaaaaaGATAATTAGGAAGGAAAGTAATGCAATTATCGAAAATACAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGtt                                <  2:3386308/87‑1 (MQ=255)
                                                                         ttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAAt                       >  2:75946/1‑100 (MQ=255)
                                                                         ttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAAt                       <  1:2283083/100‑1 (MQ=255)
                                                                         ttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAAt                       <  1:1876724/100‑1 (MQ=255)
                                                                          tgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATACAAAAATAGTATACTTTATTCCAGACGAGTTCACTCTAAt                       <  2:3481558/99‑1 (MQ=255)
                                                                           gCAAAAGAAGATTATTAGGAAGGGAACTAATGCCATTATCGACAATTAAAAAAATATCCAAAAATAGAATACTTTATTCCAGAATATTTCAATaaaa                        >  2:1625617/1‑97 (MQ=255)
                                                                           gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATg                      >  2:3674319/1‑99 (MQ=255)
                                                                           gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                     <  2:3541224/100‑1 (MQ=255)
                                                                           gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                     >  1:762777/1‑100 (MQ=255)
                                                                           gCAAAAAAAGATTATTAGGAAGGGAATTAATGAAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                     <  2:344090/100‑1 (MQ=255)
                                                                            cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                    >  1:1055778/1‑100 (MQ=255)
                                                                            cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                    >  1:3844834/1‑100 (MQ=255)
                                                                            cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                     >  2:4087787/1‑99 (MQ=255)
                                                                            cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAAAAAAAAAAAATAGTATACTTTATTCCAGAAGAGTTCAAtata                         <  2:3121653/95‑1 (MQ=255)
                                                                             aaaaaaaGATTATTAGGAAGAAAATTAATGAAATTATCGAAAATACAAAAAAAAAAAAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGttt                   <  2:517631/100‑1 (MQ=255)
                                                                             aaaaaaaGATTATAAGGAAGGGAATTAATGCAATTATCAAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATAAAATGttt                   >  2:2575870/1‑100 (MQ=255)
                                                                              aaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                  <  2:176981/100‑1 (MQ=255)
                                                                              aaaaaaGATTATTAGGAAGGGAATTAATGCAATTAAAAAAAAAAAAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                  <  2:3691074/100‑1 (MQ=255)
                                                                                aaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTc                <  2:2868766/100‑1 (MQ=255)
                                                                                   aGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTc             <  1:3592620/100‑1 (MQ=255)
                                                                                   aGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTc             <  1:817886/100‑1 (MQ=255)
                                                                                   aGATTATTAAGAAGGGAATTAATGAAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTc             <  2:1249528/100‑1 (MQ=255)
                                                                                    gATTATTAGGAAGGGAATTAATGCAATTATCGGAAATTCAAAAAATATCCAAAAATAGTATACATAATTCCAGAAGAGTTCAATATAATGTTTGTCTTCa            >  2:1322255/1‑100 (MQ=255)
                                                                                     attattAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCaa           >  1:3149342/1‑100 (MQ=255)
                                                                                     attattAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCaa           <  1:2905654/100‑1 (MQ=255)
                                                                                     attattAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCaa           >  1:2671211/1‑100 (MQ=255)
                                                                                     attattAGGAAGGGAATTAATGCAATTATCGAAAAATAAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCaa           <  2:3091571/100‑1 (MQ=255)
                                                                                        actaGGAAGGGAAATAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttt        >  1:2807584/3‑100 (MQ=255)
                                                                                          tgggaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttttt      >  1:2834500/3‑100 (MQ=255)
                                                                                          taGGAAGGGAATTAATGCAATTATCGAAAATTCAAAACATATCCAAAAATAGTATACTTTATTCCAGAAGCGTTCAAAATAATGTTTGTCTTCTAttttt      >  2:3371253/1‑100 (MQ=255)
                                                                                          taGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttttt      <  2:2113634/100‑1 (MQ=255)
                                                                                            ggaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCt    >  2:839399/1‑100 (MQ=255)
                                                                                             gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt   <  2:1457402/100‑1 (MQ=255)
                                                                                             gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt   <  2:3281/100‑1 (MQ=255)
                                                                                                  |                                                                                               
CACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTA  >  NC_000913/1197641‑1197834

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