Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB122 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
*ACB122120,9750T.95.2% 261.8 / 8.8 65intergenic (+272/+100)polA/engBDNA polymerase I/putative GTP‑binding protein EngB
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base . (35/25);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 8.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCAGAGTAACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAG  >  ACB122/120841‑121108
                                                                                                                                      |                                                                                                                                     
cTCAGAGTAACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGttt                                                                                                                                     <  6:118206/137‑1 (MQ=255)
      gTAACTGAATGGCTGGTGAAGCCCTTAGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCAGCCCTAGCGGTCCGGACTTGGACCCCGACCTTCCACCTCCCCATACGAAGTCCGGGGgtttgtttg                                                                                                                               >  6:78250/1‑136 (MQ=255)
                ggCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAAttt                                                                                                                    >  7:281396/1‑137 (MQ=255)
                ggCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAAttt                                                                                                                    >  7:184633/1‑137 (MQ=255)
                ggCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCAATACGAACTCCGTGG‑TTTTTTTTGCCCGGAAttt                                                                                                                    >  8:65556/1‑137 (MQ=255)
                ggCTGGTGAAGCCCATCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGAGTTGGACCCCGAACTTCCACCTCCCCATACGAAGACCGGGG‑TTTTTGTTGCGCGGAAttt                                                                                                                    >  1:343145/1‑137 (MQ=255)
                    ggTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCg                                                                                                                >  5:30805/1‑137 (MQ=255)
                     gTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCgg                                                                                                               >  2:230823/1‑137 (MQ=255)
                        aaGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCACCATACGAAGTCCGGGG‑TTTTTTTTGGCCGGAATTTGCCGGGAg                                                                                                            >  1:314449/1‑137 (MQ=255)
                           cccTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCgg                                                                                                         >  8:81106/1‑137 (MQ=255)
                            ccTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCggg                                                                                                        >  7:145428/1‑137 (MQ=255)
                              ttCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGtt                                                                                                      >  5:329740/1‑137 (MQ=255)
                                         tATGTTGTGTTAAGTGTTGGCAGATATCTTGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                           <  2:327261/137‑1 (MQ=255)
                                         tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                           <  4:128658/137‑1 (MQ=255)
                                         tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                           <  8:184633/137‑1 (MQ=255)
                                         tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGCTTTCTCGCGAAt                                                                                           <  2:122541/137‑1 (MQ=255)
                                           tgttgtGTTAAGTGTTCGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                         <  6:166362/137‑1 (MQ=255)
                                             ttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCAGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCgg                                                                                       >  8:311302/1‑136 (MQ=255)
                                              tgtgTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGAt                                                                                      <  6:30805/137‑1 (MQ=255)
                                                 gtTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTc                                                                                   <  2:133013/137‑1 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCg                                                                                                                          >  3:230065/1‑96 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCg                                                                                                                          <  4:230065/96‑1 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATTCGATGACGa                                                                                 >  4:256914/1‑137 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  3:89452/1‑137 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  5:142725/1‑137 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  5:108346/1‑137 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  4:63692/1‑137 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  3:83076/1‑137 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  2:166578/1‑137 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCACGCCCTAGCGGTCCTGACTTGGACCCCGAACTTCCCCCTCCCCATACTAAGTCAGGGG‑TTTTTTTTGCACGGAATTTTACGGGATATGGTTTTATCGCTAATGCGATGGCGa                                                                                 >  2:103761/1‑137 (MQ=255)
                                                             cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgcc                                                                       >  5:56308/1‑137 (MQ=255)
                                                               gATATCTGGCCCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGAGAATGCGATGTCGAATACACcgccgc                                                                     >  1:343186/1‑137 (MQ=255)
                                                                 tatCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                                   <  1:230823/137‑1 (MQ=255)
                                                                 tagcTGGACCCCGTCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                                   <  7:100530/134‑1 (MQ=255)
                                                                    cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCggg                                                                                                              <  5:251127/91‑1 (MQ=255)
                                                                    cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCggg                                                                                                              >  6:251127/1‑91 (MQ=255)
                                                                    cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTcg                                                                >  5:240088/1‑137 (MQ=255)
                                                                      ggacggcgcgCTAGCGGTCGGGACTTGGACCGCGAAGTAGCCCGTCCCCATAGGAAGTCC‑GGGTTTTTTTTTGCCCGGAATTTGGCGGGAGCGGGTTTTCTGGCGAATGCGATGTCGAATACACCGCCGCATTggcg                                                              <  8:261826/131‑1 (MQ=255)
                                                                             gCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAc                                                       >  1:243421/1‑137 (MQ=255)
                                                                               ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCggg                                                                                                              <  8:194043/80‑1 (MQ=255)
                                                                               ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCggg                                                                                                              >  7:194043/1‑80 (MQ=255)
                                                                               ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                     >  3:51482/1‑137 (MQ=255)
                                                                               ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                     >  5:260978/1‑137 (MQ=255)
                                                                               ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCAGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                     >  6:246056/1‑137 (MQ=255)
                                                                               ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGACCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                     >  3:83616/1‑137 (MQ=255)
                                                                                cTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCg                                                    <  3:63692/137‑1 (MQ=255)
                                                                                    cGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTc                                                <  4:272630/137‑1 (MQ=255)
                                                                                        ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                            <  2:243421/137‑1 (MQ=255)
                                                                                        ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                            <  1:250344/137‑1 (MQ=255)
                                                                                        ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                            <  6:329740/137‑1 (MQ=255)
                                                                                        ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                            <  1:279102/137‑1 (MQ=255)
                                                                                        ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                            <  1:156054/137‑1 (MQ=255)
                                                                                        ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                            <  8:281396/137‑1 (MQ=255)
                                                                                         cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcg                                           >  7:314836/1‑137 (MQ=255)
                                                                                                          attCCCCCTCCACATACGAAGTCCGTGG‑TTTTTTTTGCCCGTAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAcc                          <  8:314836/136‑1 (MQ=39)
                                                                                                             ccccgTCCCCATACGATGT‑CGGTGGTTTTTTTTGCCGTTATTTTGCCGGTTGGGGGTTTTCTCGCTAATGGTTTGTGGATAACCCCGCCGCATTCGGGGGTAAACCCGGGTCTGGCGCAAACCACTAAACGACCGGc                       <  4:83616/137‑1 (MQ=255)
                                                                                                             cccccTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGc                       <  6:240088/137‑1 (MQ=255)
                                                                                                             ccccaTCCCAATACGAAGT‑CGGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGTTTTTGTCGCTAATGGGATGTCGATTACACCTCCGCATTCGCGGGTAAACCGGGGTCAGGCGCAAACCACTAAACGACCGGc                       <  7:65556/137‑1 (MQ=255)
                                                                                                                 cTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGt                   >  1:42688/1‑137 (MQ=255)
                                                                                                                  tCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCGGGGTGCGGGTTTTCTCGCGATTGCGATGTCGAATACCCCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc                  <  6:130193/137‑1 (MQ=255)
                                                                                                                  tCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGTTGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc                  <  4:315194/137‑1 (MQ=255)
                                                                                                                     ccATACGAAGTCC‑GGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCTATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTc               <  7:255802/137‑1 (MQ=255)
                                                                                                                      cATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTaaa                               >  7:128710/1‑120 (MQ=255)
                                                                                                                      cATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTaaa                               <  8:128710/120‑1 (MQ=255)
                                                                                                                      cATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTcc              >  4:336209/1‑137 (MQ=255)
                                                                                                                      cATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTcc              >  6:13186/1‑137 (MQ=255)
                                                                                                                       aTACCAAGTCC‑GGGTTTTTTTTTGCCCGGATTTTGCCGGGAGGGGGTTTTCTCGCTAATGGGATGTCCAAAACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCa             <  6:270878/137‑1 (MQ=255)
                                                                                                                                  cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatcca   >  5:133193/2‑137 (MQ=255)
                                                                                                                                   ggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAg  >  2:69155/1‑137 (MQ=255)
                                                                                                                                      |                                                                                                                                     
CTCAGAGTAACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAG  >  ACB122/120841‑121108

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: ‑