Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,799,302 G→C 8.1% A330G (GCT→GGT) ‡ ansB ← glutaminase‑asparaginase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,799,3020GC 8.1% 193.4 / 6.1 62A330G (GCT→GGT) ‡ansBglutaminase‑asparaginase
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base C (2/3);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

GGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATACCGAACCGT  >  NC_002947/2799162‑2799449
                                                                                                                                               |                                                                                                                                                   
ggCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGAcc                                                                                                                                             <  3:347424/149‑1 (MQ=255)
ggCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGAcc                                                                                                                                             <  2:307609/149‑1 (MQ=255)
 gCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGGCTCACGAcct                                                                                                                                            <  8:100153/149‑2 (MQ=255)
  cccTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGAcc                                                                                                                                             >  3:350471/1‑147 (MQ=255)
  cccTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGAc                                                                                                                                              >  7:165619/1‑146 (MQ=255)
  cccTCAGGGCCGGTGCGCCGGCGCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGAc                                                                                                                                              <  8:165619/146‑1 (MQ=255)
     tCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGCGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTc                                                                                                                                        >  1:127270/1‑149 (MQ=255)
     tCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTCCCATCGCCAGGATCCGCGCTTGTTGCGGGTTCCTGTCGTGAGCCCCGCCCCGGTc                                                                                                                                        >  4:310220/1‑149 (MQ=255)
     tCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTc                                                                                                                                        >  1:218856/1‑149 (MQ=255)
       agagCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGGCTCACGAcctgaacc                                                                                                                                       <  3:184753/148‑7 (MQ=255)
       agagCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGGCTCACGAcctgaacc                                                                                                                                       <  1:147810/148‑7 (MQ=255)
       agagCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTc                                                                                                                                        <  6:45834/147‑1 (MQ=255)
         agCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTc                                                                                                                                    <  4:251285/149‑1 (MQ=255)
              gcgcACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATACTCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCtt                                                                                                                                  <  1:126033/149‑1 (MQ=255)
              gcgcACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTc                                                                                                                               <  8:6304/149‑1 (MQ=255)
                   ccGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCg                                                                                                                          >  6:393763/1‑149 (MQ=255)
                         cTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGtt                                                                                                                    <  3:334203/149‑1 (MQ=255)
                                   aGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACg                                                                                                          <  2:287445/149‑1 (MQ=255)
                                     cACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCa                                                                                                        <  2:59035/149‑1 (MQ=255)
                                        tCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTa                                                                                                     >  4:274304/1‑149 (MQ=255)
                                        tCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTa                                                                                                     >  1:65935/1‑149 (MQ=255)
                                        tCAGTACTCCCAGAAGATCCGCTGCAGCTCTTAGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCAACGCCAGGATACTCCCGTATTGCGGGATCACGACGTGAGCCACGCCCCAGTCGTTCTTGTAATCCGGATCTACGCCAATACGCAGTa                                                                                                     >  5:13162/1‑149 (MQ=255)
                                          aGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTTAGCCACGACCCAGTCGTTCTTGTCATACGGCTGTTCGGCATTACGCa                                                                                                        <  6:173918/144‑1 (MQ=255)
                                          aGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTAc                                                                                                    >  2:251122/1‑148 (MQ=255)
                                          aGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACg                                                                                                   <  2:218856/149‑1 (MQ=255)
                                          aGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTAATCCGGCTGTTCGGCATTACGCa                                                                                                        >  5:173918/1‑144 (MQ=255)
                                              cTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAac                                                                                               >  6:352143/1‑149 (MQ=255)
                                                  cAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCtt                                                                                                                                  <  3:225884/110‑1 (MQ=255)
                                                  cAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCtt                                                                                                                                  >  4:225884/1‑110 (MQ=255)
                                                     aagaTCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGACACGACCCAGTCGTTATTGTCGTCCGGCTGTTCGGCATTACGCAGTACGAAACCAcccgg                                                                                        >  8:149133/1‑147 (MQ=255)
                                                             cTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACAt                                                                                <  3:171228/149‑1 (MQ=255)
                                                                 aGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGcga                                                                            <  5:86852/149‑1 (MQ=255)
                                                                        tGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCg                                                                      >  3:176315/1‑148 (MQ=255)
                                                                          cTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGga                                                                    <  5:352143/148‑1 (MQ=255)
                                                                          cTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGACATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGgat                                                                   >  4:202949/1‑149 (MQ=255)
                                                                             tCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGgatgat                                                                <  4:125905/149‑1 (MQ=255)
                                                                               cTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATAAGGCACTACGAACCCACCCGGGGTGACATGCGCCGACCGGATGctct                                                              >  4:98628/1‑149 (MQ=255)
                                                                                  ggTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGc                                                                                                                        <  5:176665/88‑1 (MQ=255)
                                                                                  ggTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGc                                                                                                                        >  6:176665/1‑88 (MQ=255)
                                                                                    gcTTGGGCA‑‑‑TCGCTACCATCGCCAGTATCAGAGCTTTTTGAGGTTTCAGGTCGTGAGCCACGACCCAGTAGGTCTTGTAATCCTGCTTTTCGGCATTACTAAGTACGAAACCACACTTTTTGACATGCGACGAGCTGATGATCTGAACg                                                         <  6:77500/148‑1 (MQ=255)
                                                                                            a‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCtt                                                 <  4:176936/149‑1 (MQ=255)
                                                                                            a‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCtt                                                 <  2:127270/149‑1 (MQ=255)
                                                                                                 cGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTgc                                               >  2:153853/1‑149 (MQ=255)
                                                                                                  gCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTgc                                               <  4:350471/148‑1 (MQ=255)
                                                                                                      ccATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATc                                                               >  8:54455/1‑128 (MQ=255)
                                                                                                      ccATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATc                                                               <  7:54455/128‑1 (MQ=255)
                                                                                                      ccATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCa                                            <  6:161755/147‑1 (MQ=255)
                                                                                                       cATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGtt                                         <  2:65935/149‑1 (MQ=255)
                                                                                                        aTCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTc                                        >  5:370269/1‑149 (MQ=255)
                                                                                                                  tCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGcc                              <  3:310220/149‑1 (MQ=255)
                                                                                                                    cgcgCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCgg                                                                     <  8:82996/108‑1 (MQ=255)
                                                                                                                    cgcgCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCgg                                                                     >  7:82996/1‑108 (MQ=255)
                                                                                                                     gcgcTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCgg                            >  5:34645/1‑148 (MQ=255)
                                                                                                                     gcgcTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCgg                            >  3:216295/1‑148 (MQ=255)
                                                                                                                                    tCAGGTCGTGAGCCACGACCCAGTCGTTATTGTCATCCGGCTGTTCGGCATTACTCAGTACGACAACACCCTGGTTGACATGCGACGCGCGGATGATCTGAACGCCGTTCTTGGGCAGTTCTTGCAGCGCCGGGACAAGCCGcgccga             >  4:10237/1‑148 (MQ=255)
                                                                                                                                     tgggTCGTGGCTCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATa           >  1:305971/3‑149 (MQ=255)
                                                                                                                                      aGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATAc          >  1:374120/1‑149 (MQ=255)
                                                                                                                                      aGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATAc          <  5:393763/149‑1 (MQ=255)
                                                                                                                                       ggTCGTGGCTCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGGATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATAcc         >  1:286200/1‑149 (MQ=255)
                                                                                                                                           gTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCgg                            <  6:280702/126‑1 (MQ=255)
                                                                                                                                           gTGAGCCACGACCCAGTCGTTCTTGTCATACGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCgg                            >  5:280702/1‑126 (MQ=255)
                                                                                                                                              aGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACACTGGTTGACATGCGAAGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACCATACCGAACCGt  >  5:38445/1‑149 (MQ=255)
                                                                                                                                               gCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGtt                                                                                     >  8:144837/1‑65 (MQ=255)
                                                                                                                                               gCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGtt                                                                                     <  7:144837/65‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATACCGAACCGT  >  NC_002947/2799162‑2799449

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