Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB143 120,983 (T)9→8 intergenic (+280/+92) polA → / ← engB DNA polymerase I/putative GTP‑binding protein EngB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143120,9750T.93.9% 270.0 / 13.1 67intergenic (+272/+100)polA/engBDNA polymerase I/putative GTP‑binding protein EngB
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base . (48/14);  total (51/16)
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

AACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATC  >  ACB143/120849‑121105
                                                                                                                              |                                                                                                                                  
aaCTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGttttttttgcc                                                                                                                          >  1:21769/1‑134 (MQ=255)
   tGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCgg                                                                                                                      >  2:145396/1‑137 (MQ=255)
      aTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTAGGACCCCGAACTTCCCCCTCCCCATACGAAGACCGGGGTTTTTTTTGCCCGGAAt                                                                                                                   >  8:148351/1‑137 (MQ=255)
        ggCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCACCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAAttt                                                                                                                 >  1:162541/1‑137 (MQ=255)
          cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAAttt                                                                                                                 <  5:318969/135‑1 (MQ=255)
          cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAAttt                                                                                                                 >  6:318969/1‑135 (MQ=255)
          cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTgg                                                                                                               >  3:205573/1‑136 (MQ=255)
          cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCACCCCTAGCGCTCCGGACTTGGACCCCGAACTTCCACCTCACCATACGAAGTCCGGGGTTTTTATTGACCGGAATTTgg                                                                                                               >  3:260819/1‑136 (MQ=255)
            ggTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCg                                                                                                             >  8:205674/1‑137 (MQ=255)
              tGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACACCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCACTCCCCATACGAAGTCAGGGGTTTTTTTTGCCAGGAATTTGCCggg                                                                                                           >  7:222835/1‑137 (MQ=255)
               gAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGa                                                                                                          >  5:55246/1‑137 (MQ=255)
                  gCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCg                                                                                                       >  6:178962/1‑137 (MQ=255)
                               ccTATGTTGTGTTAAGTGTTGGCAGATATCTGGCCCCCGCCCTAGCGGTCCGGACTTGGACCTCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGa                                                                                          <  8:294117/137‑1 (MQ=255)
                                 tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  1:301612/137‑1 (MQ=255)
                                   tgttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCACTCCCCATACGAAGTACGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                      >  4:289966/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCTGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                              >  5:115335/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGATTGCGATGTCGa                                                                              >  1:298557/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                              >  3:116051/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                              >  4:291461/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                              >  3:87490/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                              >  6:75779/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGACCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATTCGATGTCGa                                                                              >  7:207846/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCAGGGGTTTTTTTTTCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCTAATGCGATGTCGa                                                                              >  8:85580/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCAGAACTACCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGACCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGAttgcga                                                                              >  2:325960/1‑132 (MQ=255)
                                           tAAGTGTTGGCAGATATCTCGACCCCGCCCTAGCGGTCCGGACTTGGACACCGAAATTCCCCCTCCCCATACGAAGTCCCGGGTTTTTATTTAAAGGAATTTTCCGGGAGCGGGTTTTCTCGCGAAATCGATGTAGa                                                                              >  5:34435/1‑137 (MQ=255)
                                                   ggCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccg                                                                      <  3:291461/137‑1 (MQ=255)
                                                     cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCGGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAcctcc                                                                    >  7:225554/1‑137 (MQ=255)
                                                     cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACACCGAACTTCCCCCTCAACCTACGAAGTCCGGGGTTTTTTTTTACCGGAATTTGAAGGGCTCGGGTTTTCTCGCGAAATCTATGTCGAATACAcagcc                                                                    >  4:251992/1‑137 (MQ=255)
                                                      aGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATAAAccgccg                                                                   >  6:145189/1‑137 (MQ=255)
                                                       gATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACTAAGTCCGGGGTTTTTTTTGCAAGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAAAAACcgccgc                                                                  >  5:61487/1‑137 (MQ=255)
                                                        atatCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCa                                                                 <  4:190025/137‑1 (MQ=255)
                                                         tatCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                                <  5:269544/137‑1 (MQ=255)
                                                           tCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTc                                                              <  5:218635/137‑1 (MQ=255)
                                                           tCTGGACCCCGCACTAGCGGTCCGGACTTGGACCCCGAACTTACACCTCCACAAACGAAGTCCGGGGTTTTTTTTGCCCGGAATATGCCGGGAGCGGGTTTTCTAGCGAATTCGCGGTAGAATACACCGCCGCAtgc                                                              >  2:183753/1‑135 (MQ=255)
                                                                 ccccGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGt                                                        >  5:165979/1‑137 (MQ=255)
                                                                 ccccGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGt                                                        >  3:115512/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                  >  7:189584/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                  >  6:286360/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                  >  5:295748/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                  >  5:63985/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                  >  2:269562/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                  >  1:185405/1‑137 (MQ=255)
                                                                         tAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCgg                                                <  8:133325/137‑1 (MQ=255)
                                                                         tAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCgg                                                <  4:98496/137‑1 (MQ=255)
                                                                                 cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcg                                        <  6:63985/137‑1 (MQ=255)
                                                                                 cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcg                                        >  2:312348/1‑137 (MQ=255)
                                                                                 cGGACTTGGACCCCGAACTTCCCCCTACCAATACGAAGTCCGTGGTTTTTTTTGCCCGGAATTTGCCGGGATCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCAGGGTCatgag                                        >  2:307721/1‑133 (MQ=38)
                                                                                  ggACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcgc                                       >  6:267535/1‑137 (MQ=255)
                                                                                   gACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGCATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGCGTCAGGCGCa                                      >  6:270559/1‑137 (MQ=255)
                                                                                    aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCaa                                     >  3:108705/1‑137 (MQ=255)
                                                                                             ccGAACTTCCCCCTCCCATTAGGAAGTCCGGGTTTTTTTTTGCCCGTAATTTGCGGGGAGGGGGTTTTCTCGCTAATGCGATGTCTAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTaaa                            <  4:116051/137‑1 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCCCTAAACGAc                        >  7:123455/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                        >  6:163087/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                        >  6:157809/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                        >  1:49006/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                        >  1:2600/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTAGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                        >  6:2252/1‑137 (MQ=255)
                                                                                                   ttCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  2:53902/71‑1 (MQ=255)
                                                                                                   ttCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        >  1:53902/1‑71 (MQ=255)
                                                                                                     cccccTCCCCATACTAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCTATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGc                    <  1:104430/137‑1 (MQ=255)
                                                                                                          tCCCCATACTAAGTCGGGTGTTTTTTTTGCCCGTAATTTGCCGGGAGCGGGTTTTCTCGCTAAGGCTATGTCGAATACACCGCCGCATTCGCGGTTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc               <  8:207846/137‑1 (MQ=255)
                                                                                                          tCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc               <  8:92693/137‑1 (MQ=255)
                                                                                                             ccATACGAAGTCCGCGGTTTTTTTTGCCCGGAATTTGCCGGCAGCGGGTTTTCTCGCGAATGAGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCTGCTTGTCttg            >  6:230245/1‑136 (MQ=255)
                                                                                                              cATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTcc           >  7:28382/1‑137 (MQ=255)
                                                                                                                    aaGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        >  7:249224/1‑54 (MQ=255)
                                                                                                                    aaGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  8:249224/54‑1 (MQ=255)
                                                                                                                    aaGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtcc     >  1:240175/1‑137 (MQ=255)
                                                                                                                       tCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatc  >  7:261659/1‑137 (MQ=255)
                                                                                                                              |                                                                                                                                  
AACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATC  >  ACB143/120849‑121105

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: