Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,120,133 T→C T263A (ACG→GCG)  ftsN ← essential cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,120,1330TC99.3% 461.9 / ‑4.6 146T263A (ACG→GCG) ftsNessential cell division protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (81/63);  total (83/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 159 total aligned reads displayed.

CTTTCACCGGGCCAATGACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAG  >  NC_000913/4120043‑4120222
                                                                                          |                                                                                         
cTTTCACCGGGCCAATGACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGcct                                                                                         >  1:4743029/1‑93 (MQ=255)
cTTTCACCGGGCCAATGACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGcct                                                                                         <  1:1601203/93‑1 (MQ=255)
 tttCACCGGGCCAATGACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGcctc                                                                                        >  1:3060822/1‑93 (MQ=255)
  ttCACCGGGCCAATGACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCctct                                                                                       <  1:1447489/93‑1 (MQ=255)
   tCACCGGGCCAATGACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTg                                                                                      >  2:3691991/1‑93 (MQ=255)
     aCCGGGCCAATGACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGcc                                                                                    <  1:2605992/93‑1 (MQ=255)
      ccGGGCCAATGACCACACGATTCCAGCCATTGTTGGTGGTGGTTTTCGAGTCAAAGCCTTCTAACGCCAGCTGATCACGTACCGCctcttgct                                                                                   >  2:1229929/1‑89 (MQ=255)
      ccGGGCCAATGACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCt                                                                                   >  2:2444058/1‑93 (MQ=255)
      ccGGGCCAATGACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCt                                                                                   <  1:3242354/93‑1 (MQ=255)
         ggCCAATGACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGtt                                                                                <  1:2365691/93‑1 (MQ=255)
          gCCAATGACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTc                                                                               <  1:2160847/93‑1 (MQ=255)
           ccAATGACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTcc                                                                              >  2:3895160/1‑93 (MQ=255)
               tGACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCACAGCCTTCGAACGCCAGCTGAGCACGTACCGACTCTGCCTGTTCCgcgc                                                                          >  1:4285958/1‑93 (MQ=255)
                gACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACTGCCTCTGCCCGTTCCGCGcc                                                                         >  2:3073429/1‑93 (MQ=255)
                gACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGcc                                                                         <  1:118343/93‑1 (MQ=255)
                gACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGcc                                                                         <  1:791817/93‑1 (MQ=255)
                gACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGcc                                                                         >  2:3073428/1‑93 (MQ=255)
                gACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGcc                                                                         <  1:1384007/93‑1 (MQ=255)
                  ccACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCctc                                                                       >  1:339141/1‑93 (MQ=255)
                    ccacGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTg                                                                     <  1:114714/92‑1 (MQ=255)
                    acacGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTg                                                                     >  2:3901697/1‑93 (MQ=255)
                     cacGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTga                                                                    <  2:1468013/93‑1 (MQ=255)
                     cacGATGCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTga                                                                    <  1:3377580/93‑1 (MQ=255)
                      acGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTgaa                                                                   >  2:2831198/1‑93 (MQ=255)
                      acGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTgaa                                                                   >  2:512465/1‑93 (MQ=255)
                       cGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTgaac                                                                  <  2:499456/93‑1 (MQ=255)
                        gATTCCAGCCATTGTTGGTGGTGATTTTAGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTgaacg                                                                 <  1:782414/93‑1 (MQ=255)
                        gATACCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTgaacg                                                                 <  2:2858743/93‑1 (MQ=255)
                         aTTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTgaacga                                                                <  1:2444058/93‑1 (MQ=255)
                          ttCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTgaacgaa                                                               >  2:4256098/1‑93 (MQ=255)
                          ttCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTgaacgaa                                                               >  1:3173093/1‑93 (MQ=255)
                          ttCCAGCCATTGTTGGTGGTGAGTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTgaacgaa                                                               >  2:2989713/1‑93 (MQ=255)
                           tCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTTCCTGTTCCGCGCCTCTgaacgaac                                                              <  1:2438343/93‑1 (MQ=255)
                           tCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTgaacgaac                                                              <  1:2874731/93‑1 (MQ=255)
                           tCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTgaacgaac                                                              <  1:3337452/93‑1 (MQ=255)
                              aGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGc                                                           <  2:4519704/93‑1 (MQ=255)
                              aGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGc                                                           >  1:4268388/1‑93 (MQ=255)
                               gCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCa                                                          <  1:2758361/93‑1 (MQ=255)
                                ccATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCAc                                                         <  1:717708/93‑1 (MQ=255)
                                ccATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCAc                                                         <  1:1873192/93‑1 (MQ=255)
                                 cATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACt                                                        >  1:3756371/1‑93 (MQ=255)
                                 cATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACt                                                        >  2:3550797/1‑93 (MQ=255)
                                 cATTGTTGGGGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTTGACGAACCGCAcg                                                        >  2:3371182/1‑92 (MQ=255)
                                    tgttgGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCa                                                     <  1:529964/93‑1 (MQ=255)
                                       tggtggTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCa                                                  <  2:3109947/93‑1 (MQ=255)
                                        ggtggtGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCAt                                                 <  2:2561388/93‑1 (MQ=255)
                                        ggtggtGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCAt                                                 <  1:561194/93‑1 (MQ=255)
                                          tggtgATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATcc                                               <  2:2737329/93‑1 (MQ=255)
                                          tggtgATTGTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATcc                                               <  1:614328/93‑1 (MQ=255)
                                           ggtgATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCa                                              <  1:2259022/93‑1 (MQ=255)
                                           ggtgATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCa                                              >  2:1721473/1‑93 (MQ=255)
                                            gtgATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAg                                             >  1:2018558/1‑93 (MQ=255)
                                             tgATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAgc                                            >  1:2829302/1‑93 (MQ=255)
                                             tgATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAgc                                            >  2:1937799/1‑93 (MQ=255)
                                             tgATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAgc                                            <  2:1374634/93‑1 (MQ=255)
                                               aTTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAgcgg                                          >  2:958544/1‑93 (MQ=255)
                                               aTTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAgcgg                                          >  1:2682864/1‑93 (MQ=255)
                                               aTTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAgcgg                                          <  1:3597737/93‑1 (MQ=255)
                                                  ttCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGt                                       <  1:3442268/93‑1 (MQ=255)
                                                   tCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGtt                                      <  1:1702865/93‑1 (MQ=255)
                                                   tCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGtt                                      <  1:4161934/93‑1 (MQ=255)
                                                    cGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTc                                     <  1:2173820/93‑1 (MQ=255)
                                                    cGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTc                                     >  1:662041/1‑93 (MQ=255)
                                                    cGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCAGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTc                                     >  1:662040/1‑93 (MQ=255)
                                                       gTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTc                                  <  2:1035041/93‑1 (MQ=255)
                                                        tCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCt                                 >  2:378035/1‑93 (MQ=255)
                                                        tCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCt                                 >  2:2033800/1‑93 (MQ=255)
                                                         cacaGCATTCGACCGCCAGCTGAGCACGTACCGCCTATGCCTGATCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGAGGCGTTCGTCtt                                <  2:4218017/90‑1 (MQ=255)
                                                          aaaGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCttt                               >  1:4016657/1‑93 (MQ=255)
                                                          aaaGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCttt                               >  1:3042390/1‑93 (MQ=255)
                                                          aaaGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCttt                               <  1:235441/93‑1 (MQ=255)
                                                           aggCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCCGCGGCGTTCGTCtttt                              <  2:3130473/91‑1 (MQ=255)
                                                           aaGCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCtttt                              <  1:2628330/93‑1 (MQ=255)
                                                             gCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCtttttt                            <  2:3201361/93‑1 (MQ=255)
                                                             gCCTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCtttttt                            >  2:2456813/1‑93 (MQ=255)
                                                               cTTCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTtct                          >  2:2585531/1‑93 (MQ=255)
                                                                ttCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTtctc                         >  1:3018953/1‑93 (MQ=255)
                                                                 tCGAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTcc                        >  1:489784/1‑93 (MQ=255)
                                                                   gAACGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTccgc                      <  2:2328618/93‑1 (MQ=255)
                                                                    aaCGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTccgcc                     >  1:2694999/1‑93 (MQ=255)
                                                                    aaCGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTccgcc                     >  1:458796/1‑93 (MQ=255)
                                                                    aaCGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTccgcc                     <  1:895953/93‑1 (MQ=255)
                                                                     aCGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTccgccg                    >  2:2118441/1‑93 (MQ=255)
                                                                      cGCCAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCACACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGt                   >  2:4538602/1‑93 (MQ=255)
                                                                        ccAGCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCg                 <  1:3779027/93‑1 (MQ=255)
                                                                           gCTGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGtt              >  1:1589140/1‑93 (MQ=255)
                                                                             tGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGtttt            >  1:1682336/1‑93 (MQ=255)
                                                                             tGAGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGtttt            >  1:199953/1‑93 (MQ=255)
                                                                               aGCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTgg          >  1:1210443/1‑93 (MQ=255)
                                                                                gCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGgc         <  1:895206/93‑1 (MQ=255)
                                                                                gCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGgc         <  2:3483315/93‑1 (MQ=255)
                                                                                gCACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTCgc         >  2:1316019/1‑93 (MQ=255)
                                                                                 cACGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGgcg        <  1:4572031/93‑1 (MQ=255)
                                                                                  aCGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGgcgc       >  1:1608847/1‑93 (MQ=255)
                                                                                   cGTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGcgcg      <  1:1452699/93‑1 (MQ=255)
                                                                                    gTACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCGCCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGt     <  2:2210335/93‑1 (MQ=255)
                                                                                     tACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTc    >  1:3946199/1‑93 (MQ=255)
                                                                                     tACCGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTc    >  1:3797784/1‑93 (MQ=255)
                                                                                      cccGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCa   <  2:1164535/92‑1 (MQ=255)
                                                                                       ccGCCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAg  >  2:4386140/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTTTCACCGGGCCAATGACCACACGATTCCAGCCATTGTTGGTGGTGATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAG  >  NC_000913/4120043‑4120222

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