Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,881,483 Δ1 bp intergenic (‑126/‑482) dnaA ← / → rpmH DnaA/50S ribosomal subunit protein L34

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,881,4830G.99.7% 1652.4 / ‑3.9 350intergenic (‑126/‑482)dnaA/rpmHDnaA/50S ribosomal subunit protein L34
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (178/171);  total (179/171)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 356 total aligned reads displayed.

GGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGAGTCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTA  >  NC_000913/3881395‑3881570
                                                                                        |                                                                                       
ggCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAg                                                                                    >  2:4267715/1‑93 (MQ=255)
ggCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAg                                                                                    <  2:6492524/93‑1 (MQ=255)
ggCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAg                                                                                    >  2:6613668/1‑93 (MQ=255)
 gCTGAAACATTCATGATTGTTGACGTCCGTCGAAAAGACCCTGCTTGAGGGTGACGCTCGAACCGCTGTCTGCGGTTATATGCCCGATCAAGa                                                                                   <  2:5807245/93‑1 (MQ=255)
   tGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATc                                                                                 >  1:3973804/1‑93 (MQ=255)
   tGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATc                                                                                 <  2:4020370/93‑1 (MQ=255)
    gAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATcc                                                                                <  1:3233563/93‑1 (MQ=255)
     aaaCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCt                                                                               >  2:1768596/1‑93 (MQ=255)
     aaaCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCt                                                                               >  2:1071277/1‑93 (MQ=255)
      aaaaTTCATGATTGTTGACGTACGTCGAAAAGAACCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTg                                                                              >  1:3193212/4‑93 (MQ=255)
      aaCATTCATGATTGTTGACGTACGTCGAAAATGCCCTTCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTg                                                                              <  2:2238509/93‑1 (MQ=255)
      aaCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTg                                                                              >  2:1010590/1‑93 (MQ=255)
      aaCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTg                                                                              <  1:4126313/93‑1 (MQ=255)
      aaCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTg                                                                              >  2:2668633/1‑93 (MQ=255)
       aCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGc                                                                             >  1:3523774/1‑93 (MQ=255)
        cATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCa                                                                            >  1:3848420/1‑93 (MQ=255)
         aTTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaa                                                                           >  2:4781083/1‑93 (MQ=255)
         aTTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaa                                                                           >  1:3036276/1‑93 (MQ=255)
          ttCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaa                                                                          >  2:3435232/1‑93 (MQ=255)
          ttCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaa                                                                          >  2:3654733/1‑93 (MQ=255)
          ttCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaa                                                                          >  1:2607194/1‑93 (MQ=255)
           tCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCTGTTATATGCCCGATCAAGATCCTGCaaaa                                                                         <  1:6186041/93‑1 (MQ=255)
            cATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                                        <  2:2178705/93‑1 (MQ=255)
            cATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                                        <  1:2142045/93‑1 (MQ=255)
            cATGATGGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                                        <  2:5380859/93‑1 (MQ=255)
              tGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGa                                                                      >  2:390054/1‑93 (MQ=255)
              tGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGa                                                                      >  1:6269361/1‑93 (MQ=255)
              tGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGa                                                                      >  1:6137213/1‑93 (MQ=255)
                aTTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATc                                                                    <  2:5040373/93‑1 (MQ=255)
                aTTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATc                                                                    <  2:4708492/93‑1 (MQ=255)
                aTTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATc                                                                    <  2:1458515/93‑1 (MQ=255)
                aTTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATc                                                                    <  1:5452392/93‑1 (MQ=255)
                 ttgttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCg                                                                   <  2:5437001/93‑1 (MQ=255)
                   gttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCggg                                                                 >  1:5829633/1‑93 (MQ=255)
                    ttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATGGGGa                                                                <  2:3491080/93‑1 (MQ=255)
                    ttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGa                                                                >  2:5626396/1‑93 (MQ=255)
                    ttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGa                                                                >  1:799703/1‑93 (MQ=255)
                     tgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGAc                                                               >  1:3387848/1‑93 (MQ=255)
                     tgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGAc                                                               >  1:435311/1‑93 (MQ=255)
                      gACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGAcc                                                              >  1:149013/1‑93 (MQ=255)
                       acgtacgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACcg                                                             >  2:2356898/1‑93 (MQ=255)
                       acgtacgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACcg                                                             <  1:836634/93‑1 (MQ=255)
                         gtacgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACcgcg                                                           >  2:5315207/1‑93 (MQ=255)
                           acgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGa                                                         >  2:4525657/1‑93 (MQ=255)
                           acgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGa                                                         >  1:6588170/1‑93 (MQ=255)
                            cgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGAt                                                        <  1:6351218/93‑1 (MQ=255)
                            cgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGAt                                                        >  1:4597640/1‑93 (MQ=255)
                            cgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGAt                                                        <  1:1832137/93‑1 (MQ=255)
                              tCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCa                                                      <  2:2467962/93‑1 (MQ=255)
                              tCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCa                                                      <  1:1630453/93‑1 (MQ=255)
                               cGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCAt                                                     <  1:236168/93‑1 (MQ=255)
                                gAAAATACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATa                                                    <  1:3118379/93‑1 (MQ=255)
                                gAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATa                                                    <  1:1953613/93‑1 (MQ=255)
                                gAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATa                                                    <  2:1244352/93‑1 (MQ=255)
                                gAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATa                                                    <  1:5533301/93‑1 (MQ=255)
                                 aaaaGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAg                                                   >  1:3311985/1‑93 (MQ=255)
                                    aGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCt                                                >  2:4128177/1‑93 (MQ=255)
                                      aCCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTaa                                              <  1:1036382/93‑1 (MQ=255)
                                      aCCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTaa                                              <  2:744423/93‑1 (MQ=255)
                                         cTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACt                                           >  2:2536436/1‑93 (MQ=255)
                                         cTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACt                                           >  1:3640678/1‑93 (MQ=255)
                                          tGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTg                                          >  2:6387063/1‑93 (MQ=255)
                                           gCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTgc                                         >  2:3221971/1‑93 (MQ=255)
                                            cTTGAGGGTGACGCCCGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTgcg                                        <  1:3724003/93‑1 (MQ=255)
                                             ttGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATGGGGACCGCGGATCATAGCCTAAACTgcgc                                       >  1:4230109/1‑93 (MQ=255)
                                             ttGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTgcgc                                       >  2:2014628/1‑93 (MQ=255)
                                              tGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCa                                      >  1:2830723/1‑93 (MQ=255)
                                               gAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCaa                                     <  2:4863141/93‑1 (MQ=255)
                                                   gTGACGCACCAGCCGCTGTCTGCGGTTATATGCCTGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAgaga                                 <  2:620091/93‑1 (MQ=255)
                                                    tGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAt                                <  1:1264054/93‑1 (MQ=255)
                                                     gACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtc                               <  1:5054147/93‑1 (MQ=255)
                                                     gACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtc                               >  2:2731977/1‑93 (MQ=255)
                                                     gACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtc                               <  2:5057105/93‑1 (MQ=255)
                                                     gACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtc                               >  2:703762/1‑93 (MQ=255)
                                                      aCGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtct                              <  1:5172073/93‑1 (MQ=255)
                                                      aCGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtct                              <  2:3770332/93‑1 (MQ=255)
                                                        gCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttc                            >  1:1360612/1‑93 (MQ=255)
                                                         cACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttct                           <  2:5225417/93‑1 (MQ=255)
                                                          aCGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTg                          >  2:4159107/1‑93 (MQ=255)
                                                          aCGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTg                          >  1:5896688/1‑93 (MQ=255)
                                                           cGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGt                         >  1:3794805/1‑93 (MQ=255)
                                                            gAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGtt                        >  1:2035865/1‑93 (MQ=255)
                                                             aaCCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGttt                       >  1:4729613/1‑93 (MQ=255)
                                                             aaCCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGttt                       <  2:2063435/93‑1 (MQ=255)
                                                             aaCCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGttt                       >  1:2766760/1‑93 (MQ=255)
                                                              aCCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtc                      <  1:4488986/93‑1 (MQ=255)
                                                               ccGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtct                     <  2:1868964/93‑1 (MQ=255)
                                                                 gCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTca                   <  2:2907210/93‑1 (MQ=255)
                                                                    gTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGAGCAAGAGATCTTCTGTTTCTCACAg                >  1:1739469/1‑93 (MQ=255)
                                                                     tCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGa               >  2:3199079/1‑93 (MQ=255)
                                                                       tGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGAtt             >  2:2396795/1‑93 (MQ=255)
                                                                        gCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGAttt            <  2:6048562/93‑1 (MQ=255)
                                                                           gTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTcc         <  2:1565635/93‑1 (MQ=255)
                                                                           gTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTcc         <  2:2855477/93‑1 (MQ=255)
                                                                            ttATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTccc        >  2:4444447/1‑93 (MQ=255)
                                                                             tataTGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCg       <  1:4960997/93‑1 (MQ=255)
                                                                              atatGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGa      >  2:822103/1‑93 (MQ=255)
                                                                               tatGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGAt     >  2:2119018/1‑93 (MQ=255)
                                                                                 tGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTTCCGCttt   >  2:3087388/1‑93 (MQ=255)
                                                                                 tGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGAttt   <  1:6053847/93‑1 (MQ=255)
                                                                                        |                                                                                       
GGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGAGTCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTA  >  NC_000913/3881395‑3881570

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