Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,844,375 T→G 5.7% E166A (GAG→GCG)  ynjI ← inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,844,3750TG 5.7% 385.7 / 10.9 105E166A (GAG→GCG) ynjIinner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (50/49);  new base G (2/4);  total (52/53)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 107 total aligned reads displayed.

CAACAAGGCGATAACCTGCTCTGGGGCTTCCTCCATGTTTGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGAATCGATAAGGTGCGATCAATTG  >  NC_000913/1844236‑1844520
                                                                                                                                           |                                                                                                                                                 
caacaaGGCGATAACCTGCTCTGGGGCTTCCTCCATGTTTGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCtt                                                                                                                                           <  1:428617/148‑1 (MQ=255)
 aacaaGGCGATAACCTGCTCTGGGGCTTCCTCCATGTTTGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAACGACGCGAGTCttt                                                                                                                                          <  1:845734/148‑1 (MQ=255)
    aaGGCGATAACCTGCTCTGGGGCTTCCTCCATGTTTGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACGCGAGTCtttaga                                                                                                                                       <  1:601905/148‑4 (MQ=255)
              ccTGCTCTGGGGCTTCCTCCATGTTTGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGat                                                                                                                                <  2:989284/145‑1 (MQ=255)
                tGCTCTGGGGCTTCCTCCATGTTTGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACt                                                                                                                                      >  1:1108792/1‑137 (MQ=255)
                tGCTCTGGGGCTTCCTCCATGTTTGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACt                                                                                                                                      <  2:1108792/137‑1 (MQ=255)
                                ccATGTTTGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTa                                                                                                                                        >  2:488936/1‑119 (MQ=255)
                                ccATGTTTGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTa                                                                                                                                        >  3:609/1‑119 (MQ=255)
                                ccATGTTTGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTa                                                                                                                                        <  1:488936/119‑1 (MQ=255)
                                ccATGTTTGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTa                                                                                                                                        <  4:609/119‑1 (MQ=255)
                                      ttGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATAACCCAATGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATAGAGTCAAAAGGACGCACATg                                                                                                     >  1:867351/1‑148 (MQ=255)
                                         cTTTAAAGGTATATGCTCCATGGTCGCCAGAGAGAAAATGCTCCATTAAGGCACAATGACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGata                                                                                                  <  2:867351/148‑1 (MQ=255)
                                           ttAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCtt                                                                                                                                           >  3:52428/1‑105 (MQ=255)
                                           ttAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCtt                                                                                                                                           <  4:52428/105‑1 (MQ=255)
                                           ttAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATatcga                                                                                                                           >  3:838283/1‑121 (MQ=255)
                                           ttAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATatcga                                                                                                                           <  4:838283/121‑1 (MQ=255)
                                                  tATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGt                                                                                                                                              >  1:1068388/1‑95 (MQ=255)
                                                  tATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGt                                                                                                                                              <  2:1068388/95‑1 (MQ=255)
                                                        cTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATg                                                                                                     <  1:499570/130‑1 (MQ=255)
                                                        cTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATg                                                                                                     >  2:499570/1‑130 (MQ=255)
                                                         tCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGt                                                                                                                                              <  4:283928/88‑1 (MQ=255)
                                                         tCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGt                                                                                                                                              >  3:283928/1‑88 (MQ=255)
                                                             tGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACt                                                                                                                                                   >  4:613137/1‑79 (MQ=255)
                                                             tGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAAGAACTTGCGCTATCTTCGATACCCCATTGATCCTCTAAAGACt                                                                                                                                                   <  3:613137/79‑1 (MQ=255)
                                                                   ccAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATatcga                                                                                                                           <  1:808975/97‑1 (MQ=255)
                                                                   ccAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATatcga                                                                                                                           >  2:808975/1‑97 (MQ=255)
                                                                    cAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATatcgatcg                                                                                                                        >  2:231693/1‑99 (MQ=255)
                                                                    cAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATatcgatcg                                                                                                                        <  1:231693/99‑1 (MQ=255)
                                                                     agaaagaaAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGaa                                                                                            <  3:30957/126‑1 (MQ=255)
                                                                     agaaagaaAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGaa                                                                                            >  4:30957/1‑126 (MQ=255)
                                                                         agaaaATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGa                                                                                             >  2:68779/1‑121 (MQ=255)
                                                                         agaaaATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGa                                                                                             <  1:68779/121‑1 (MQ=255)
                                                                               tGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGcc                                                            >  4:192318/1‑148 (MQ=255)
                                                                               tGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGcc                                                            <  3:192318/148‑1 (MQ=255)
                                                                                gCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGa                                                                                             <  3:241464/114‑1 (MQ=255)
                                                                                gCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGa                                                                                             >  4:241464/1‑114 (MQ=255)
                                                                                   ccATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCtt                                                                                                                                           >  4:715040/1‑65 (MQ=255)
                                                                                   ccATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCtt                                                                                                                                           <  3:715040/65‑1 (MQ=255)
                                                                                     aTTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACt                                                      <  4:679394/148‑1 (MQ=255)
                                                                                       tAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATatcgat                                                                                                                          <  1:242819/78‑1 (MQ=255)
                                                                                       tAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATatcgat                                                                                                                          >  2:242819/1‑78 (MQ=255)
                                                                                       tAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCtt                                                                                    <  2:920347/116‑1 (MQ=255)
                                                                                       tAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATAtt                                                                                    >  1:920347/1‑116 (MQ=255)
                                                                                          ggCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATatcga                                                                                                                           >  3:551435/1‑74 (MQ=255)
                                                                                          ggCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATatcga                                                                                                                           <  4:551435/74‑1 (MQ=255)
                                                                                           gCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGa                                                                                                             >  3:994385/1‑87 (MQ=255)
                                                                                            cacaATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGa                                                                                             >  3:770110/1‑102 (MQ=255)
                                                                                            cacaATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGa                                                                                             <  4:770110/102‑1 (MQ=255)
                                                                                            cacaATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCa                                                                         >  2:659732/1‑122 (MQ=255)
                                                                                            cacaATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCa                                                                         <  1:659732/122‑1 (MQ=255)
                                                                                            cacaATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGt                                               >  1:1098236/1‑148 (MQ=255)
                                                                                             acaATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAg                                                                              >  3:907239/1‑116 (MQ=255)
                                                                                             acaATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAg                                                                              <  4:907239/116‑1 (MQ=255)
                                                                                              caATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGttt                                             >  3:937053/1‑148 (MQ=255)
                                                                                                aTAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTg                                                                                  >  1:105592/1‑109 (MQ=255)
                                                                                                aTAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTg                                                                                  <  2:105592/109‑1 (MQ=255)
                                                                                                aTAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGtt                                              >  3:520694/1‑145 (MQ=255)
                                                                                                 tAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGt                                                                                                                     <  1:430267/73‑1 (MQ=255)
                                                                                                 tAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGt                                                                                                                     >  2:430267/1‑73 (MQ=255)
                                                                                                  aaCTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGcc                                                                                >  3:607786/1‑109 (MQ=255)
                                                                                                  aaCTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGcc                                                                                <  4:607786/109‑1 (MQ=255)
                                                                                                     tttCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGc                                                             <  2:718210/125‑1 (MQ=255)
                                                                                                     tttCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGc                                                             >  1:718210/1‑125 (MQ=255)
                                                                                                      ttCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCAt                                     >  2:331001/1‑148 (MQ=255)
                                                                                                              cTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATg                             <  2:1098236/148‑1 (MQ=255)
                                                                                                                tCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGa                                                                                             <  2:832619/82‑1 (MQ=255)
                                                                                                                tCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGa                                                                                             >  1:832619/1‑82 (MQ=255)
                                                                                                                tCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAg                                                                              <  4:155975/97‑1 (MQ=255)
                                                                                                                    tACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGttt                                             <  2:625212/126‑1 (MQ=255)
                                                                                                                    tACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGttt                                             >  1:625212/1‑126 (MQ=255)
                                                                                                                       cccATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATAtt                                                                                                <  2:1052696/72‑1 (MQ=255)
                                                                                                                       cccATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATAtt                                                                                                >  1:1052696/1‑72 (MQ=255)
                                                                                                                       cccATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGAATCg                    <  1:331001/148‑1 (MQ=255)
                                                                                                                        ccATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGa                                                                                             >  3:1118088/1‑74 (MQ=255)
                                                                                                                        ccATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGa                                                                                             <  4:1118088/74‑1 (MQ=255)
                                                                                                                        ccATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAgg                                                                       <  4:540965/96‑1 (MQ=255)
                                                                                                                        ccATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAgg                                                                       >  3:540965/1‑96 (MQ=255)
                                                                                                                        ccATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATg                             <  2:45477/138‑1 (MQ=255)
                                                                                                                        ccATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGAATCGa                   >  2:130338/1‑148 (MQ=255)
                                                                                                                          aTTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTa                                                                                        >  3:1107796/1‑77 (MQ=255)
                                                                                                                          aTTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTa                                                                                        <  4:1107796/77‑1 (MQ=255)
                                                                                                                           ttGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGaa                                                                                            <  1:311212/72‑1 (MQ=255)
                                                                                                                           ttGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGaa                                                                                            >  2:311212/1‑72 (MQ=255)
                                                                                                                            tGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGAATCGATAAg               >  3:397340/1‑148 (MQ=255)
                                                                                                                            tGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATATTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGAATCGATAAg               <  4:937053/148‑1 (MQ=255)
                                                                                                                                ccTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGa                            <  4:572088/131‑1 (MQ=255)
                                                                                                                                ccTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGa                            >  3:572088/1‑131 (MQ=255)
                                                                                                                                ccTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGAATCGATAAGGTGc           >  1:157159/1‑148 (MQ=255)
                                                                                                                                ccTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGAATCGATAAGGTGc           >  4:1027839/1‑148 (MQ=255)
                                                                                                                                 agtaaaaGACGCGAGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTa                               <  1:472150/123‑1 (MQ=255)
                                                                                                                                   ctAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCa                                                                         <  2:111437/83‑1 (MQ=255)
                                                                                                                                   ctAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCa                                                                         >  1:111437/1‑83 (MQ=255)
                                                                                                                                   ctAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGa                        >  4:534336/1‑132 (MQ=255)
                                                                                                                                   ctAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGa                        <  3:534336/132‑1 (MQ=255)
                                                                                                                                      aaGACGCGAGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTc                                            >  1:747309/1‑109 (MQ=255)
                                                                                                                                      aaGACGCGAGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTc                                            <  2:747309/109‑1 (MQ=255)
                                                                                                                                       aGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGAATCGATAAGGTGCGATCAAt    <  3:1027839/148‑1 (MQ=255)
                                                                                                                                         aCTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTc                                            >  2:217920/1‑106 (MQ=255)
                                                                                                                                         aCTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTc                                            <  1:217920/106‑1 (MQ=255)
                                                                                                                                         aCTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGAATCGATAAGGTGCGATCAATTg  <  1:130338/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
CAACAAGGCGATAACCTGCTCTGGGGCTTCCTCCATGTTTGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCACAATAACTTTCGCTATCTTCGATACCCCATTGATCCTCTAAAGACTCGCGTCTTTTACTTATGATATCGATCGAGTCAAAAGGAAGCACATGATATTGGAAGGTATCTTTGCCAGGTTCAGGCTTTCGCGGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGAATCGATAAGGTGCGATCAATTG  >  NC_000913/1844236‑1844520

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