Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 132,467 G→C 5.6% intergenic (+41/‑111) cc → / → dsbA cytochrome c4/thiol:disulfide interchange protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947132,4670GC 5.6% 225.5 / 3.6 72intergenic (+41/‑111)cc/dsbAcytochrome c4/thiol:disulfide interchange protein
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base C (2/2);  total (36/36)
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.15e-01

ACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTCGC  >  NC_002947/132332‑132612
                                                                                                                                            |                                                                                                                                                 
aCGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGTCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCtttca                                                                                                                                      <  3:60430/149‑3 (MQ=255)
    tGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtcc                                                                                                                                   <  4:154681/148‑1 (MQ=255)
      gCGATGCCATGACCATGCGCACCATCTCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCATCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgt                                                                                                                                 <  5:30774/148‑1 (MQ=255)
      gCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgtc                                                                                                                                <  7:268222/149‑1 (MQ=255)
       cGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgtcc                                                                                                                               >  1:355248/1‑149 (MQ=255)
            ccATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAg                                                                                                                           <  2:71209/148‑1 (MQ=255)
             cATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGcc                                                                                                                         >  4:150431/1‑149 (MQ=255)
              aTGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCt                                                                                                                                               >  3:28648/1‑126 (MQ=255)
              aTGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCt                                                                                                                                               <  4:28648/126‑1 (MQ=255)
                 aCCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTa                                                                                                                     <  3:269213/149‑1 (MQ=255)
                  ccATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCtttca                                                                                                                                      <  3:30559/131‑3 (MQ=255)
                  ccATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCtttca                                                                                                                                      >  4:30559/1‑129 (MQ=255)
                     tGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTAGGGATTAAATTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTac                                                                                                                    >  1:116012/1‑148 (MQ=255)
                          aCCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTagatg                                                                                                            >  5:366035/1‑145 (MQ=255)
                            cATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAAc                                                                                                          >  6:214764/1‑149 (MQ=255)
                              tCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTc                                                                                                        >  3:363263/1‑149 (MQ=255)
                                 cTGGCAAGCTGCGCAACCATGACATCGAAGCGTGGGCCAGCTACATCCAGGCCCT‑‑GCATTAACTATCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCTGTCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAg                                                                                                     <  8:330403/149‑1 (MQ=255)
                                 cTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAg                                                                                                     <  5:327335/149‑1 (MQ=255)
                                    gCAAGCTGAGCAACCATGACATCCAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCCGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGGCCGCAGCCGTTACACTACAGAACTCGAGc                                                                                                    <  6:294630/147‑1 (MQ=255)
                                     cAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGTGAGAATTTTTTTGGCAAAATTATTAAAGCGCGGCCGGGGGGCCTTTTTTGCAGTGCAAACCGGTTACAATACACAACTCGAg                                                                                                     >  1:111459/1‑148 (MQ=255)
                                      aaGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCtt                                                                                                >  6:197296/1‑149 (MQ=255)
                                            aGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGc                                                                                          <  2:121512/149‑1 (MQ=255)
                                             gCAACCATGCCATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGcc                                                                                         <  2:305838/149‑1 (MQ=255)
                                             gCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGcc                                                                                         <  6:186124/149‑1 (MQ=255)
                                              cAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCt                                                                                        >  5:18861/1‑149 (MQ=255)
                                                 ccATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGtc                                                                                     >  2:323574/1‑149 (MQ=255)
                                                  cATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGGCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTTCAGAACTCGAGCCCTTCCCGGCCTGtct                                                                                    <  1:213096/149‑1 (MQ=255)
                                                     gACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAg                                                                                 >  5:130551/1‑149 (MQ=255)
                                                        aTCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTg                                                                               >  6:219719/1‑148 (MQ=255)
                                                        aTCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTg                                                                               >  3:68713/1‑148 (MQ=255)
                                                        aTCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTg                                                                               >  2:235372/1‑148 (MQ=255)
                                                          cGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGAGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAg                                                                            >  6:376062/1‑149 (MQ=255)
                                                           gAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCTGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAg                                                                            <  2:355248/148‑1 (MQ=255)
                                                                gTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCTGGCCTGTCTCAGTGCAGGTcgcg                                                                      <  1:283636/149‑1 (MQ=255)
                                                                    gCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGTCCTTCCCGGCCTGTCTCAGTGCAGGTCGcgtcgt                                                                  >  1:325562/1‑149 (MQ=255)
                                                                      cAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCAGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAAATCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTgg                                                                >  1:322073/1‑149 (MQ=255)
                                                                       aGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGc                                                               >  1:364360/1‑149 (MQ=255)
                                                                           aCATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGAcc                                                           >  4:324510/1‑149 (MQ=255)
                                                                            cATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGAc                                                            <  3:132670/147‑1 (MQ=255)
                                                                            cATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGAc                                                            <  6:76878/147‑1 (MQ=255)
                                                                             aTCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCa                                                         >  3:143300/1‑149 (MQ=255)
                                                                                cAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGa                                                      <  1:210915/149‑1 (MQ=255)
                                                                                   ggCCT‑‑GCAATAACTTTCGGT‑‑TAATGATGTAGGTCAATGATGAAAGGGCGGCCGGGCTCCCCTTTTTTCCGTCCGCAGCCGTTACAGTACAGAACTCGAGCCCATCCCGGCCTGACTCTCTCCACCTCGCGACGTGGCGACCCATAACTg                                                   >  2:98555/1‑149 (MQ=255)
                                                                                    gCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCAGTCCGCAGACGTTACACTACAGAACTCGAGCCCTTACCGGCCTGTCTCAGTGCAGGTCGAGTCGTGGCGCCCCATGACTGc                                                  >  4:34541/1‑149 (MQ=255)
                                                                                           cATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCa                                               <  6:130551/147‑1 (MQ=255)
                                                                                           cATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCa                                               <  4:68713/147‑1 (MQ=255)
                                                                                           cATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAg                                              <  5:336586/148‑1 (MQ=255)
                                                                                            aTTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAg                                              >  1:282051/1‑147 (MQ=255)
                                                                                            aTTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAg                                              >  2:63205/1‑147 (MQ=255)
                                                                                            aTTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGa                                            <  4:363263/149‑1 (MQ=255)
                                                                                            aTTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGa                                            <  2:364360/149‑1 (MQ=255)
                                                                                                  tttCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAg                                      >  7:104696/1‑149 (MQ=255)
                                                                                                    tCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGc                                     <  5:197296/148‑1 (MQ=255)
                                                                                                     cGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCa                                    <  6:168848/148‑1 (MQ=255)
                                                                                                       gT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTTCACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGc                                 >  2:361474/1‑149 (MQ=255)
                                                                                                       gT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGc                                 <  5:219719/149‑1 (MQ=255)
                                                                                                       gT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGc                                 <  3:150431/149‑1 (MQ=255)
                                                                                                              tgttgtAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTa                                                                                                                >  4:82735/1‑65 (MQ=255)
                                                                                                              tgttgtAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTa                                                                                                                <  3:82735/65‑1 (MQ=255)
                                                                                                              tgttgtAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTaaa                            >  5:210040/1‑149 (MQ=255)
                                                                                                              tgttgtAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTaa                             >  3:66796/1‑148 (MQ=255)
                                                                                                               gttgtAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCTTCCGCAGCCGTTACACTACAGAACTCGAGCCCATCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTaaa                            <  2:302480/148‑1 (MQ=255)
                                                                                                                     ggTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATtc                     <  2:116012/149‑1 (MQ=255)
                                                                                                                         cAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATtct                    >  4:202044/1‑147 (MQ=255)
                                                                                                                            ggaaaaaaGGGCAGCCCGGCCGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGGGTAAACTGATTCTCAgcg               >  7:309127/6‑149 (MQ=255)
                                                                                                                             gatgaAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgc              >  7:185393/1‑149 (MQ=255)
                                                                                                                              atgaAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgct             <  6:144769/149‑1 (MQ=255)
                                                                                                                                gaAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGcgtc                                                                    <  7:59253/92‑1 (MQ=255)
                                                                                                                                gaAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGcgtc                                                                    >  8:59253/1‑92 (MQ=255)
                                                                                                                                  aaGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTTCACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgctgcgc         <  1:361474/149‑1 (MQ=255)
                                                                                                                                    gggCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCACGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGcgctgcgctg       <  8:104696/149‑1 (MQ=255)
                                                                                                                                         gcCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTcgc  >  7:274831/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
ACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCT‑‑GCATTAACTTTCGGT‑‑TAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTCGC  >  NC_002947/132332‑132612

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