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)
   tGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCgg                                                                                                                      >  2:145396/1‑137 (MQ=255)
      aTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTAGGACCCCGAACTTCCCCCTCCCCATACGAAGACCGGGG‑TTTTTTTTGCCCGGAAt                                                                                                                   >  8:148351/1‑137 (MQ=255)
        ggCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCACCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAAttt                                                                                                                 >  1:162541/1‑137 (MQ=255)
          cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAAttt                                                                                                                 <  5:318969/135‑1 (MQ=255)
          cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAAttt                                                                                                                 >  6:318969/1‑135 (MQ=255)
          cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTgg                                                                                                               >  3:205573/1‑136 (MQ=255)
          cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCACCCCTAGCGCTCCGGACTTGGACCCCGAACTTCCACCTCACCATACGAAGTCCGGGG‑TTTTTATTGACCGGAATTTgg                                                                                                               >  3:260819/1‑136 (MQ=255)
            ggTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCg                                                                                                             >  8:205674/1‑137 (MQ=255)
              tGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACACCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCACTCCCCATACGAAGTCAGGGG‑TTTTTTTTGCCAGGAATTTGCCggg                                                                                                           >  7:222835/1‑137 (MQ=255)
               gAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGa                                                                                                          >  5:55246/1‑137 (MQ=255)
                  gCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCg                                                                                                       >  6:178962/1‑137 (MQ=255)
                               ccTATGTTGTGTTAAGTGTTGGCAGATATCTGGCCCCCGCCCTAGCGGTCCGGACTTGGACCTCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGa                                                                                          <  8:294117/137‑1 (MQ=255)
                                 tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  1:301612/137‑1 (MQ=255)
                                   tgttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCACTCCCCATACGAAGTACGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                      >  4:289966/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCTGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                              >  5:115335/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGATTGCGATGTCGa                                                                              >  1:298557/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                              >  3:116051/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                              >  4:291461/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                              >  3:87490/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                              >  6:75779/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGACCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATTCGATGTCGa                                                                              >  7:207846/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCAGGGGTTTTTTTTT‑CCCGGAATTTGCCGGGAGCGGGTTTTCTCGCTAATGCGATGTCGa                                                                              >  8:85580/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCAGAACTACCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGACCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGAttgcga                                                                              >  2:325960/1‑132 (MQ=255)
                                           tAAGTGTTGGCAGATATCTCGACCCCGCCCTAGCGGTCCGGACTTGGACACCGAAATTCCCCCTCCCCATACGAAGTCCCGGGTTTTTATTT‑AAAGGAATTTTCCGGGAGCGGGTTTTCTCGCGAAATCGATGTAGa                                                                              >  5:34435/1‑137 (MQ=255)
                                                   ggCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccg                                                                      <  3:291461/137‑1 (MQ=255)
                                                     cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCGGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAcctcc                                                                    >  7:225554/1‑137 (MQ=255)
                                                     cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACACCGAACTTCCCCCTCAACCTACGAAGTCCGGGGTTTTTTTTT‑ACCGGAATTTGAAGGGCTCGGGTTTTCTCGCGAAATCTATGTCGAATACAcagcc                                                                    >  4:251992/1‑137 (MQ=255)
                                                      aGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATAAAccgccg                                                                   >  6:145189/1‑137 (MQ=255)
                                                       gATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACTAAGTCCGGGG‑TTTTTTTTGCAAGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAAAAACcgccgc                                                                  >  5:61487/1‑137 (MQ=255)
                                                        atatCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCa                                                                 <  4:190025/137‑1 (MQ=255)
                                                         tatCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                                <  5:269544/137‑1 (MQ=255)
                                                           tCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTc                                                              <  5:218635/137‑1 (MQ=255)
                                                           tCTGGACCCCGCACTAGCGGTCCGGACTTGGACCCCGAACTTACACCTCCACAAACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATATGCCGGGAGCGGGTTTTCTAGCGAATTCGCGGTAGAATACACCGCCGCAtgc                                                              >  2:183753/1‑135 (MQ=255)
                                                                 ccccGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGt                                                        >  5:165979/1‑137 (MQ=255)
                                                                 ccccGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGt                                                        >  3:115512/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                  >  7:189584/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                  >  6:286360/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                  >  5:295748/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                  >  5:63985/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                  >  2:269562/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                  >  1:185405/1‑137 (MQ=255)
                                                                         tAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCgg                                                <  8:133325/137‑1 (MQ=255)
                                                                         tAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCgg                                                <  4:98496/137‑1 (MQ=255)
                                                                                 cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcg                                        <  6:63985/137‑1 (MQ=255)
                                                                                 cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcg                                        >  2:312348/1‑137 (MQ=255)
                                                                                 cGGACTTGGACCCCGAACTTCCCCCTACCAATACGAAGTCCGTGG‑TTTTTTTTGCCCGGAATTTGCCGGGATCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCAGGGTCatgag                                        >  2:307721/1‑133 (MQ=38)
                                                                                  ggACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcgc                                       >  6:267535/1‑137 (MQ=255)
                                                                                   gACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGCATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGCGTCAGGCGCa                                      >  6:270559/1‑137 (MQ=255)
                                                                                    aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCaa                                     >  3:108705/1‑137 (MQ=255)
                                                                                             ccGAACTTCCCCCTCCCATTAGGAAGTCC‑GGGTTTTTTTTTGCCCGTAATTTGCGGGGAGGGGGTTTTCTCGCTAATGCGATGTCTAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTaaa                            <  4:116051/137‑1 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCCCTAAACGAc                        >  7:123455/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                        >  6:163087/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                        >  6:157809/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                        >  1:49006/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                        >  1:2600/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTAGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                        >  6:2252/1‑137 (MQ=255)
                                                                                                   ttCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  2:53902/71‑1 (MQ=255)
                                                                                                   ttCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        >  1:53902/1‑71 (MQ=255)
                                                                                                     cccccTCCCCATACTAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCTATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGc                    <  1:104430/137‑1 (MQ=255)
                                                                                                          tCCCCATACTAAGT‑CGGGTGTTTTTTTTGCCCGTAATTTGCCGGGAGCGGGTTTTCTCGCTAAGGCTATGTCGAATACACCGCCGCATTCGCGGTTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc               <  8:207846/137‑1 (MQ=255)
                                                                                                          tCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc               <  8:92693/137‑1 (MQ=255)
                                                                                                             ccATACGAAGTCCGCGG‑TTTTTTTTGCCCGGAATTTGCCGGCAGCGGGTTTTCTCGCGAATGAGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCTGCTTGTCttg            >  6:230245/1‑136 (MQ=255)
                                                                                                              cATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTcc           >  7:28382/1‑137 (MQ=255)
                                                                                                                    aaGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        >  7:249224/1‑54 (MQ=255)
                                                                                                                    aaGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  8:249224/54‑1 (MQ=255)
                                                                                                                    aaGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtcc     >  1:240175/1‑137 (MQ=255)
                                                                                                                       tCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatc  >  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: ‑