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.2 / ‑5.4 357intergenic (‑126/‑482)dnaA/rpmHDnaA/50S ribosomal subunit protein L34
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (172/184);  minor base A (0/1);  total (172/185)
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 357 total aligned reads displayed.

GGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGAGTCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTA  >  NC_000913/3881395‑3881570
                                                                                        |                                                                                       
ggCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTTACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAg                                                                                    <  1:6211098/93‑1 (MQ=255)
ggCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAg                                                                                    >  2:5072780/1‑93 (MQ=255)
ggCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAg                                                                                    >  1:5806861/1‑93 (MQ=255)
 gCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGa                                                                                   <  1:1521791/93‑1 (MQ=255)
   tGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATc                                                                                 >  1:6583212/1‑93 (MQ=255)
    gAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATcc                                                                                <  1:2927766/93‑1 (MQ=255)
     aaaCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCt                                                                               <  1:6607457/93‑1 (MQ=255)
      aaCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTg                                                                              >  2:2379175/1‑93 (MQ=255)
      aaCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTg                                                                              >  2:2076855/1‑93 (MQ=255)
      aaCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTg                                                                              >  2:3474590/1‑93 (MQ=255)
       aCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGc                                                                             <  1:5660083/93‑1 (MQ=255)
         aTTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaa                                                                           >  1:3499938/1‑93 (MQ=255)
          ttCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaa                                                                          >  2:6261107/1‑93 (MQ=255)
          ttCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaa                                                                          <  1:5220392/93‑1 (MQ=255)
           tCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaaa                                                                         >  1:3795232/1‑93 (MQ=255)
           tCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaaa                                                                         <  2:6471938/93‑1 (MQ=255)
            cATGATTGTTGACGTACGTCGAAAAGTCCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                                        <  1:5338760/93‑1 (MQ=255)
            cATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                                        >  2:7029406/1‑93 (MQ=255)
            cATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                                        <  1:4246061/93‑1 (MQ=255)
             aTGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACg                                                                       <  1:3182412/93‑1 (MQ=255)
              tGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCGCGATCAAGATCCTGCAAAACGa                                                                      >  1:1615167/1‑93 (MQ=255)
              tGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGa                                                                      <  2:542068/93‑1 (MQ=255)
              tGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGa                                                                      >  1:3293610/1‑93 (MQ=255)
              tGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGa                                                                      <  1:2736126/93‑1 (MQ=255)
               gATTGTTGACGTACGTCGAAACGCCCCTGCTTGAGGGTGACGCCCGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGAt                                                                     <  1:965900/93‑1 (MQ=255)
                aTTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATc                                                                    <  1:193674/93‑1 (MQ=255)
                aTTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATc                                                                    >  2:1218484/1‑93 (MQ=255)
                 ttgttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCg                                                                   <  1:2550371/93‑1 (MQ=255)
                 ttgttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCg                                                                   >  1:6711750/1‑93 (MQ=255)
                   gttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCggg                                                                 <  2:3553836/93‑1 (MQ=255)
                   gttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCggg                                                                 >  2:948164/1‑93 (MQ=255)
                     tgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGAc                                                               <  1:6098864/93‑1 (MQ=255)
                     tgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGAc                                                               <  2:4380421/93‑1 (MQ=255)
                      gACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGAcc                                                              <  1:4082348/93‑1 (MQ=255)
                      gACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGAcc                                                              >  1:246181/1‑93 (MQ=255)
                         gtacgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACcgcg                                                           <  2:3256122/93‑1 (MQ=255)
                          tacgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCgg                                                          <  1:2878801/93‑1 (MQ=255)
                           acgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGa                                                         <  2:3127683/93‑1 (MQ=255)
                             gtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATc                                                       >  2:1803068/1‑93 (MQ=255)
                              tCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCa                                                      >  2:6150009/1‑93 (MQ=255)
                              tCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCa                                                      <  2:2187789/93‑1 (MQ=255)
                               cGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCAt                                                     <  2:6808621/93‑1 (MQ=255)
                                gAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATa                                                    >  2:3842589/1‑93 (MQ=255)
                                 aaaaGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAg                                                   <  1:3420049/93‑1 (MQ=255)
                                 aaaaGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAg                                                   >  2:4902205/1‑93 (MQ=255)
                                  aaaGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGc                                                  <  2:6891276/93‑1 (MQ=255)
                                   aaGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGcc                                                 <  2:6576281/93‑1 (MQ=255)
                                    aGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCt                                                <  1:5726654/93‑1 (MQ=255)
                                     gACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTa                                               >  2:4029119/1‑93 (MQ=255)
                                      aCCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTaa                                              >  1:1186861/1‑93 (MQ=255)
                                       cccTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTaaa                                             >  1:5991380/1‑93 (MQ=255)
                                        ccTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCTTGCAAAACGATCGGGACCGCGGATCATAGCCTAAAc                                            <  1:2303867/93‑1 (MQ=255)
                                        ccTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAAc                                            <  2:4754219/93‑1 (MQ=255)
                                        ccTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAAc                                            <  2:7517149/93‑1 (MQ=255)
                                        ccTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAAc                                            >  1:3054075/1‑93 (MQ=255)
                                        ccTACTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAAc                                            >  1:6976717/1‑93 (MQ=255)
                                         cTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACt                                           <  2:7047917/93‑1 (MQ=255)
                                          tGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTg                                          >  1:4910216/1‑93 (MQ=255)
                                          tGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTg                                          >  1:756848/1‑93 (MQ=255)
                                             ttGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTgcgc                                       >  1:2824836/1‑93 (MQ=255)
                                               gAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCaa                                     >  1:1995120/1‑93 (MQ=255)
                                               gAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCaa                                     <  2:5344757/93‑1 (MQ=255)
                                                aGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAag                                    <  1:7105260/93‑1 (MQ=255)
                                                 gggTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAaga                                   <  1:1802615/93‑1 (MQ=255)
                                                  ggTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAagag                                  <  1:1898974/93‑1 (MQ=255)
                                                   gTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAgaga                                 >  1:7194672/1‑93 (MQ=255)
                                                     gACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtc                               <  1:2110374/93‑1 (MQ=255)
                                                     gACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtc                               >  1:5032485/1‑93 (MQ=255)
                                                     gACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtc                               <  2:5120917/93‑1 (MQ=255)
                                                     gACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtc                               <  1:6900840/93‑1 (MQ=255)
                                                        gCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttc                            >  1:1365696/1‑93 (MQ=255)
                                                        gCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttc                            >  1:1730997/1‑93 (MQ=255)
                                                        gCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttc                            <  2:4433085/93‑1 (MQ=255)
                                                         cACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCGTAGCCTAAACTGCGCAAGAGAtcttct                           <  2:5642520/93‑1 (MQ=255)
                                                         cACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttct                           <  1:1116031/93‑1 (MQ=255)
                                                         cACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttct                           >  2:5786747/1‑93 (MQ=255)
                                                          aCGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTg                          >  2:5896058/1‑93 (MQ=255)
                                                          aCGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTg                          >  1:5104924/1‑93 (MQ=255)
                                                           cGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGt                         <  2:6050203/93‑1 (MQ=255)
                                                           cGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGt                         >  2:2661119/1‑93 (MQ=255)
                                                           cGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGt                         <  2:2857298/93‑1 (MQ=255)
                                                            gAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGtt                        >  1:1243910/1‑93 (MQ=255)
                                                             aaCCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGttt                       >  1:6564030/1‑93 (MQ=255)
                                                                cGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtctc                    <  1:7329817/93‑1 (MQ=255)
                                                                 gCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTca                   >  2:1501644/1‑93 (MQ=255)
                                                                   tGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTcaca                 <  2:3639581/93‑1 (MQ=255)
                                                                   tGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTcaca                 <  1:3926792/93‑1 (MQ=255)
                                                                       tGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGAtt             >  2:3022740/1‑93 (MQ=255)
                                                                        gCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGAttt            <  1:5531866/93‑1 (MQ=255)
                                                                         cGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGAtttt           <  2:2708126/93‑1 (MQ=255)
                                                                          ggTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTc          >  1:1040634/1‑93 (MQ=255)
                                                                          ggTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCAGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTc          <  2:5198605/93‑1 (MQ=255)
                                                                           gTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTcc         >  2:3735607/1‑93 (MQ=255)
                                                                           gTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTcc         <  1:4565132/93‑1 (MQ=255)
                                                                            ttATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTctc        >  2:4587455/1‑91 (MQ=255)
                                                                            ttATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTccc        <  2:1054198/93‑1 (MQ=255)
                                                                             tataTGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCg       >  1:6325690/1‑93 (MQ=255)
                                                                                atGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGAtt    >  1:4496899/1‑93 (MQ=255)
                                                                                 tGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGAttt   <  2:7204765/93‑1 (MQ=255)
                                                                                  gCCCGATCAAGATCCTGCAAAACGATCGGGATCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTa  >  1:7284729/1‑93 (MQ=255)
                                                                                        |                                                                                       
GGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGAGTCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTA  >  NC_000913/3881395‑3881570

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