Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,555,389 T→A 15.9% intergenic (+27/+47) sbcB → / ← PP_1366 exodeoxyribonuclease I/transcriptional regulator MvaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,555,3890TA15.9% 159.7 / 15.7 67intergenic (+27/+47)sbcB/PP_1366exodeoxyribonuclease I/transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base T (32/24);  new base A (7/4);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTC  >  NC_002947/1555242‑1555537
                                                                                                                                                   |                                                                                                                                                      
gggCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGTAGGCCTGGCAACTGCATGCTCTACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAAt                                                                                                                                                        <  3:1714/148‑1 (MQ=255)
    gACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGcc                                                                                                                                                   >  1:51344/1‑149 (MQ=255)
      cTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCgg                                                                                                                                                 <  8:314169/149‑1 (MQ=255)
       ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGc                                                                                                                                                >  1:405690/1‑149 (MQ=255)
         cGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGc                                                                                                                                                <  1:374699/147‑1 (MQ=255)
          gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGc                                                                                                                                                >  2:305418/1‑146 (MQ=255)
          gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGc                                                                                                                                                <  1:305418/146‑1 (MQ=255)
            aaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAAt                                                                                                                                                        >  4:238010/1‑136 (MQ=255)
            aaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAAt                                                                                                                                                        <  3:238010/136‑1 (MQ=255)
             aaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑Gttttt                                                                                                                                        <  6:196050/149‑1 (MQ=255)
             aaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑Gttttt                                                                                                                                        <  2:51344/149‑1 (MQ=255)
             aaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑Gttttt                                                                                                                                        <  3:124464/149‑1 (MQ=255)
               ggCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTttgt                                                                                                                                      >  1:118692/1‑149 (MQ=255)
                    gcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGc                                                                                                                                 >  2:325383/1‑149 (MQ=255)
                                 ggAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGc                                                                                                                    <  3:371075/149‑1 (MQ=255)
                                   aGGGCGCCGACGAGGCAGGCGGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGt                                                                                                                  >  1:202199/1‑149 (MQ=255)
                                       cgcCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGc                                                                                                              <  7:158164/149‑1 (MQ=255)
                                       cgcCGACGAGGCAGGCAGGCGCGTGTTGGAGGACTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTTATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGc                                                                                                              <  1:727/149‑1 (MQ=255)
                                             cgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATa                                                                                                         <  6:30261/148‑1 (MQ=255)
                                                ggcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATTCGGATAAt                                                                                                       >  7:57048/1‑147 (MQ=255)
                                                 gcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAg                                                                                                                             <  7:341868/124‑1 (MQ=255)
                                                 gcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAg                                                                                                                             >  8:341868/1‑124 (MQ=255)
                                                   aggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTc                                                                                                  <  4:215990/149‑1 (MQ=255)
                                                    gtcaggcGCGTGTTGCAGGCCTGGCAACTGCATTCTCGACAATTTCAGAAACGCTATTCAATCGGCTGATGGCAGATAAAAAAAAAAGCCGGCAATGTCCAAGC‑‑TTTATTATTGTTAGAACCAAGAGCACGAAAGCCGAAAATCAGTc                                                                                                  >  4:110016/3‑148 (MQ=255)
                                                      caggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcag                                                                                               >  6:189314/1‑149 (MQ=255)
                                                          cgcgTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAgg                                                                                           <  3:177828/149‑1 (MQ=255)
                                                             gtgtTGGAGTCATGGGAACTGCATGCGCCACAATTGCAGGAACGCTATTCATTCGGATGATGGCAGATATACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAATCAAGAGGGCGTATGCGGATAATCAGACCATCAggg                                                                                          <  6:356447/147‑1 (MQ=255)
                                                             gtgtTGGAGGCCTGGGAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGtt                                                                                        >  5:101564/1‑149 (MQ=255)
                                                             gtgtTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAggg                                                                                          <  3:110016/147‑1 (MQ=255)
                                                             gtgtTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGtt                                                                                        >  6:327403/1‑149 (MQ=255)
                                                               gtTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGc                                                                                      >  1:377804/1‑149 (MQ=255)
                                                                   gAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAg                                                                                  >  2:388237/1‑149 (MQ=255)
                                                                       ccTGGCAACTGCATGCTCGTCAATTGCAGGAACGCTATTGAATCGGGTGATGGCATATAAACAAAAACGCCGGCAATGTGCCGGCTAGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCt                                                                                <  1:325383/149‑1 (MQ=255)
                                                                        cTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCAAGCttt                                                                              >  3:224993/1‑148 (MQ=255)
                                                                         tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTc                                                                             >  8:80431/1‑148 (MQ=255)
                                                                            cAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAAcc                                                                         <  2:118692/149‑1 (MQ=255)
                                                                             aaCTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCa                                                                        <  4:224993/149‑1 (MQ=255)
                                                                                tGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGt                                                                     >  6:204421/1‑149 (MQ=255)
                                                                                tGCATGCTCGACAATTACAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGt                                                                     >  4:203903/1‑149 (MQ=255)
                                                                                   aTGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTcg                                                                   <  5:300220/148‑1 (MQ=255)
                                                                                   aTGCTCGACAATTGCAGGAACCCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTcg                                                                   >  6:300220/1‑148 (MQ=255)
                                                                                         gACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGccc                                                             >  7:102422/1‑148 (MQ=255)
                                                                                         gACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGccc                                                             >  2:20772/1‑148 (MQ=255)
                                                                                                               aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCATGGTTGCCCAGCTTTCAACCACGTCGACGCCCCAAGTGGCTTTCCACTATTTCAg                                      >  3:52931/1‑149 (MQ=255)
                                                                                                                ctCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGACGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                     <  8:349679/148‑1 (MQ=255)
                                                                                                                aTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                     >  5:290984/1‑149 (MQ=255)
                                                                                                                    gCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCt                                 >  6:243284/1‑149 (MQ=255)
                                                                                                                        aTGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtgtg                             <  5:327403/149‑1 (MQ=255)
                                                                                                                             aGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCttt                                                  >  1:193115/1‑123 (MQ=255)
                                                                                                                             aGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCttt                                                  <  2:193115/123‑1 (MQ=255)
                                                                                                                              gATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgcc                       <  5:189314/149‑1 (MQ=255)
                                                                                                                               cacaaCAAAAACGCCGGCACATTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCACAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTATCAGGGTCTTGTGGTTgccg                      >  6:215413/4‑149 (MQ=39)
                                                                                                                               aTAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCAAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgccg                      >  5:162024/1‑149 (MQ=255)
                                                                                                                                  aaCAAAAACGCCGGCACATTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtgtg                             <  2:128908/139‑1 (MQ=255)
                                                                                                                                  aaCAAAAACGCCGGCACATTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGCTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtgtg                             >  1:128908/1‑139 (MQ=255)
                                                                                                                                   aCAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCAGGTCGCCGCCCCAAGTGGCTTTCGACTCTTTCAGGGTCTTGTGGTTgccgccgt                  >  7:170288/1‑147 (MQ=255)
                                                                                                                                   aCAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCtt                  >  8:281326/1‑149 (MQ=255)
                                                                                                                                   aCAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCtt                  >  4:354182/1‑149 (MQ=255)
                                                                                                                                     aaaaaCGCCGGCACATTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTg                >  2:234636/1‑149 (MQ=255)
                                                                                                                                     aaaaaCGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                     <  1:122531/128‑1 (MQ=255)
                                                                                                                                     aaaaaCGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                     >  2:122531/1‑128 (MQ=255)
                                                                                                                                      aaaaCGCCGGCACATTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg               >  7:293207/1‑149 (MQ=255)
                                                                                                                                      aaaaCGCCGGCACATTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg               >  5:261781/1‑149 (MQ=255)
                                                                                                                                         aCGCCGGCACATTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg               >  7:263023/1‑146 (MQ=255)
                                                                                                                                         aCGCCGGCACATTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg               <  8:263023/146‑1 (MQ=255)
                                                                                                                                          cGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCgg                                                                                                            >  6:280011/1‑52 (MQ=255)
                                                                                                                                          cGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCgg                                                                                                            <  5:280011/52‑1 (MQ=255)
                                                                                                                                           gCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTc           >  3:211/1‑148 (MQ=255)
                                                                                                                                            ccGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGa         <  1:20772/149‑1 (MQ=255)
                                                                                                                                             cGGCAATGTGACGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTACAGCAGGGTTGCCCAGCTTTAAAACACGTCGCCGCCCCAACTGGGTTTCAAATCTTTCATGGTCTTGTGGTTGCCGCCTTTGGTTTCGa         >  6:281703/1‑148 (MQ=255)
                                                                                                                                                   tgtgCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTc  <  2:202199/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGC‑‑GTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTC  >  NC_002947/1555242‑1555537

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: