Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,602,700 T→C 16.6% Q16R (CAG→CGG)  ftsE ← cell division ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,602,7000TC16.6% 151.0 / 8.2 72Q16R (CAG→CGG) ftsEcell division ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (35/24);  new base C (2/11);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 5.34e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTCATGGCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGA  >  NC_000913/3602552‑3602840
                                                                                                                                                    |                                                                                                                                            
gTCATGTCCGCTAAACCAGATTTTCCCGGCGCTGGACCGCTCACTCCCACAGCTGAGCTTCTGGAGGGTGCTTTTCCCTGCGCCGTAATGAGCGGTCAGAAACGCCGTCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCt                                                                                                                                              <  2:951551/149‑1 (MQ=255)
 tCATGGCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGGACTTTTCCCTGCGCCGGAATGACCGGACAGAAACGCCATCTCACCCGGTGGCATATGGAACGTAACGCCCTGCAGCGCCTg                                                                                                                                             >  2:716825/1‑149 (MQ=255)
  cATGGCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGt                                                                                                                                            <  2:319869/149‑1 (MQ=255)
     ggCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTACAGCGCCTGtcta                                                                                                                                         >  2:649957/1‑148 (MQ=255)
      gCCGCTAAACCAGCGTTTCCCGGCGCTGGGCCGCTCATCCCAACGGGTCAGCTTCAGGAGCGGACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGGCTcc                                                                                                                                        <  2:392364/149‑1 (MQ=255)
        cgcgAAACCAGGTTGTCCCGGCGCTGGGCCGCTCAATCCGACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAGCGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTcccc                                                                                                                                      <  2:18904/149‑2 (MQ=255)
          cTAAACCGGATTTTCCCGGCGCTGGGCCGCTCAATCCCACGGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCAcc                                                                                                                                    <  1:1050716/149‑1 (MQ=255)
            aaaCCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCga                                                                                                                                  <  1:656063/149‑1 (MQ=255)
              aCCAGATTTTCCCGGCGCTGGGGCGCTCAATCCCTCAGATCAGCTTCAGGAGGGTACTTTTCCTTGCGCCGGGATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCgaga                                                                                                                                <  1:261684/149‑1 (MQ=255)
                   aTTTTCCCGGCGCTGGGCCGCTCAATCCCGCAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTGACCCGGCTGCATATGGAACGTAACGCCCTGCAGCTCCTGTCTCCCACCGAGATAAGc                                                                                                                           >  1:822896/1‑149 (MQ=255)
                   aTTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGc                                                                                                                           >  2:733215/1‑149 (MQ=255)
                       tCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCGGCGCCTGCCTCCCACCGAGATAAGCCTTg                                                                                                                       >  2:541706/1‑149 (MQ=255)
                        cccGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGCGGGTACTTTCCCCTGCGCCGCAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGc                                                                                                                      <  1:1056356/149‑1 (MQ=255)
                           gtcccTGGCCCGTTCAATCCCATAGCTCCGCCTCAGGAGCGCACTCTCCCCCTCGCCGTAATGACCGGTCATTACCCCCATCTCACCCGGCTGCATATGGAACGTTCCTCCCTTAATCATCCGGCTACCACCGAGATAAGCCTTGCTGa                                                                                                                   <  2:312128/145‑1 (MQ=255)
                                 gggCCGGTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACTCGGCTGCATATTGAACGTAACGCCCTGTAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGtt                                                                                                             <  1:457434/149‑1 (MQ=255)
                                   gCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCGCACCGAGATAAGCCTTGATGACATGTTCa                                                                                                           >  1:422885/1‑149 (MQ=255)
                                     cGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCaaa                                                                                                         <  1:68021/149‑1 (MQ=255)
                                         catgcccACAGATCAGGTTCAGGGGGGTATGTTTCCCTGCGCCGGAATGACCGGTAAGAAACGCCATCTCACCCGGCGGCAGATGGAACGTATCGTCTTGCAGGATCCGGCTACCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGa                                                                                                     <  2:564776/145‑1 (MQ=255)
                                         cAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTATTTTTTTTCATGATCCGGCGACCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGa                                                                                                     <  2:484088/149‑1 (MQ=255)
                                          aaTCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATCAGCCTTGCTGACATGTTCAAAGCGaa                                                                                                    >  1:1141052/1‑149 (MQ=255)
                                           ctCCCACAGCTCCGCTTCACGAGGCTATTCCTCCCTGCGCCTGACTCACCGTTCAGTAACCTCTTCTCCCCCGGCTGCGTATGGACCGTATTTCTTTTAATGATACGGCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAAt                                                                                                   <  2:618429/148‑1 (MQ=255)
                                              ccACAGATCAGTTCCAGGAGGGTACTTTTCCCTGCGCCGGAAGGACAGGTCAGACACTGCATCTCGTCCGGCTGCATATGGAACGTTTTTTTCTTAATGATCCGGCGCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCAt                                                                                                <  2:637159/149‑1 (MQ=255)
                                              ccACAGATCAGCTTCAGGGGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTTTCTCTTTGAAGCACCCGGCGACCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCAt                                                                                                <  2:278441/149‑1 (MQ=255)
                                              ccACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGCCGTGTTCAAAGCGAATCAt                                                                                                >  1:621684/1‑149 (MQ=255)
                                               cacaGATCAGTGTCAGGGTGGTACTTTTCCCGGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTTACTTTTTTAATGATCCGGCGACCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCAtt                                                                                               <  2:898354/149‑1 (MQ=255)
                                               cacaGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGGGATAAGCCTTGCTGACATGTTCAAAGCGAATCAtt                                                                                               >  1:1258293/1‑149 (MQ=255)
                                                    aTCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAGACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGGGATAAGCCTGGCTGACAGGTTCAAAGCGAATCGTTGTTaa                                                                                          >  1:113512/1‑149 (MQ=255)
                                                    aTCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCGCCCGGCTGCATATGGAACGTAACGCCCTGGCGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACAGGTTCAAAGCCAATCATTGTTaa                                                                                          >  1:156622/1‑149 (MQ=255)
                                                     tCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCACGAAGCGGCCGGCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAAt                                                                                         >  1:1271243/1‑149 (MQ=255)
                                                      cAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATc                                                                                        <  1:550019/149‑1 (MQ=255)
                                                      cAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATc                                                                                        <  1:926322/149‑1 (MQ=255)
                                                      cAGCGTCAGCAGCGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATc                                                                                        <  1:354314/149‑1 (MQ=255)
                                                           tCAGGAGGGTACTGTTCCGTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTTATTCTCTTCATGGCCCGGCGACCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCtctc                                                                                   <  2:623669/149‑1 (MQ=255)
                                                               gaggGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGc                                                                               >  2:1125474/1‑149 (MQ=255)
                                                                    tACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTATCTCCCTGCAGGGTCCGGCGCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCaaaaa                                                                          <  2:1289320/149‑1 (MQ=255)
                                                                     aCTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCCTTGTTAATCCTCTCGGGCAAACAg                                                                         >  1:365686/1‑149 (MQ=255)
                                                                     aCTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCCGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCCCCGAGGTAAGCCTTGCCGACAGGTTCGAAGCGACTCATTGTTAATGATGTCGGGCCAAAAg                                                                         >  2:231965/1‑149 (MQ=255)
                                                                       ttttCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAATGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTCACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTg                                                                       >  1:564363/1‑149 (MQ=255)
                                                                       ttttCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGAAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTGTCGGGCAAAAAgag                                                                       >  2:990290/1‑147 (MQ=255)
                                                                         ttCCCTGCGCCGGAATGACCCGTCAGAAACGCCATCTCACCCGGCTGCCTATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACGCGTTCAAGGCGAATCATTGTTAATCCTCTCGGGCCAACAGTGcc                                                                     >  2:1275764/1‑149 (MQ=255)
                                                                            ccTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCctct                                                                  >  2:1038387/1‑149 (MQ=255)
                                                                            ccTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGGTAATCCTCTCGGGCAAAAAGTGCctca                                                                  >  2:118830/1‑148 (MQ=255)
                                                                             cTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTCACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTCAATCCTCTCGGGCAACAAGTGCctctt                                                                 >  1:440298/1‑148 (MQ=255)
                                                                                 gcCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGACGCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATaaa                                                             >  1:982125/1‑149 (MQ=255)
                                                                                 gcCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATaaa                                                             <  1:978273/149‑1 (MQ=255)
                                                                                  ccGGAATGACCGGTCGGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAg                                                            <  1:610134/149‑1 (MQ=255)
                                                                                  ccGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAACAGTGCCGCTGTAAAg                                                            >  1:732226/1‑149 (MQ=255)
                                                                                   cGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAgt                                                           <  1:733215/149‑1 (MQ=255)
                                                                                      aaTGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAgtcgt                                                        >  2:1131245/1‑149 (MQ=255)
                                                                                      aaTGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAatcgt                                                        >  1:101995/1‑149 (MQ=255)
                                                                                       aTGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGTCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAgtcgtc                                                       <  2:1055799/149‑1 (MQ=255)
                                                                                       aTGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCCCTCGGGCAAAAAGTGCCTCTATAAAgtcgtc                                                       >  1:346938/1‑149 (MQ=255)
                                                                                           ccGGTCAGAAACGCCATCTCACGCGGCTGCCTATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGcc                                                   <  1:6029/149‑1 (MQ=255)
                                                                                            gggTCAGAAAAGCCATCTCTCCCGGCTGCATATGGACCGTAATGCCTTGCAGGGCCCGGCGCCCACCGAGATAAGCCTTGCTGACATGCTCAAAGCGAATCATTGTTAATACTCGGAAGCAAAAAGTGCCTCTATAAAGTCGTCCGCCt                                                  <  2:922483/148‑1 (MQ=255)
                                                                                             ggTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGGAAAGCCCTGCAGCGCCTCTCTCCCACCGAGATAAACCTGGCTGACATGTCCAAACCGAACCATCCTCAATCCTCTCGGGCACAAACTGCCACAGTAAAGCCATCGGCCta                                                 >  2:485308/1‑148 (MQ=255)
                                                                                                cAGAAACGCCATCTCACCCGGCTGCATATCGAACGTATCGCCCTGCAGCGCCTGGCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCACAAGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTaaa                                              <  2:746708/149‑1 (MQ=255)
                                                                                                 aGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGACCTACCGCGCCCGTATCCGTCCTACAAAATCCATGCTAACATGTTCAAAGCGAATAATCGATAATACCCTCGGGCAAAGAGGGCCTCGCTAAAGTCCGCCGCCTTaaaa                                             >  2:148909/1‑148 (MQ=255)
                                                                                                  gAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGACCCAGTTTCACATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACg                                            <  1:945075/149‑1 (MQ=255)
                                                                                                   aaaCGCCATCTCACCCGGCTGCATATGGGACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGAAAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATACTCTCGGGCAAAAAGAGCCTCAAAAAAGTCGTCCGCCTTAAACga                                           >  1:113496/1‑148 (MQ=255)
                                                                                                     aCGCCATCTCACCCGGCTGCATATGGAACGTACCGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCATTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGAc                                         <  1:256881/149‑1 (MQ=255)
                                                                                                          aTCTCACCCGGCTGCATATGGAACGTAACGCCCTGCATCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCaaa                                    >  1:827925/1‑149 (MQ=255)
                                                                                                                cccGGCTGCATATGGATCGTATCGCCCTGCAGCGTCCGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCa                              <  2:259748/149‑1 (MQ=255)
                                                                                                                 ccGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCaa                             >  1:9454/1‑149 (MQ=255)
                                                                                                                        cATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTc                      <  1:235917/149‑1 (MQ=255)
                                                                                                                          tatGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGc                    <  1:1125474/149‑1 (MQ=255)
                                                                                                                            tGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATACAGTCGTCCGCCTTAAATGGACGCAAATCCTCCATCCGTCCACCg                  >  1:253118/1‑149 (MQ=255)
                                                                                                                             ggAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGa                 >  2:877385/1‑149 (MQ=255)
                                                                                                                              gAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGac                >  2:381529/1‑149 (MQ=255)
                                                                                                                                 cGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACACGCTCAAAGCGAATCCTCGTTAATCCGCCCGAGCAAAACCTGCCTCTATAAGGCACTCCGCCTGAACAGGACGCAAACGCTCAATACGGTCGCCGGCAcc             >  2:1211677/1‑149 (MQ=255)
                                                                                                                                       gCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAAGCCTCAATCCGTGCGCCGACACCACTGTa       >  1:21883/1‑149 (MQ=255)
                                                                                                                                        cccTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAg      <  1:1112051/149‑1 (MQ=255)
                                                                                                                                         ccTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTAAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGc     <  1:425814/149‑1 (MQ=255)
                                                                                                                                            gCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCAAGTCAGTTTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGCCGCAAATCCTCAATACGTTCGCCGACCCCAATGTAGCGGa  >  1:1013782/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GTCATGGCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGA  >  NC_000913/3602552‑3602840

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