Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,504,068 T→A 6.3% intergenic (+35/‑228) rpsI → / → petA 30S ribosomal protein S9/ubiquinol‑cytochrome c reductase iron‑sulfur subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,504,0680TA 6.3% 189.1 / 2.7 63intergenic (+35/‑228)rpsI/petA30S ribosomal protein S9/ubiquinol‑cytochrome c reductase iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base A (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

CTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGACG  >  NC_002947/1503920‑1504216
                                                                                                                                                    |                                                                                                                                                      
cTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCtt                                                                                                                                                        <  7:74576/149‑1 (MQ=255)
   cGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTgg                                                                                                                                                      <  1:183849/148‑1 (MQ=255)
         gcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCggg                                                                                                                                               <  7:341041/149‑1 (MQ=255)
               cgtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑tttt                                                                                                                                        <  1:347562/148‑1 (MQ=255)
               cgtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑tttt                                                                                                                                        <  1:353045/148‑1 (MQ=255)
               cgtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑tttt                                                                                                                                        <  1:7220/148‑1 (MQ=255)
                   gtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑ttttttttt                                                                                                                                   <  3:288314/149‑1 (MQ=255)
                    tGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑tttttttt                                                                                                                                    >  5:372217/1‑147 (MQ=255)
                           tACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTg                                                                                                                            >  3:55732/1‑148 (MQ=255)
                           tACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTg                                                                                                                            >  8:246952/1‑148 (MQ=255)
                             cGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAAc                                                                                                                         >  6:340884/1‑149 (MQ=255)
                              gTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAAc                                                                                                                         >  3:70410/1‑148 (MQ=255)
                               tCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCACTACTCCAAGCGTTAATACCGCTACGGCAGTCGAAAAAAGCCCGGTGACTTGGAACCGGGC‑‑TTTTTTTATAGCCGTGAACt                                                                                                                        >  3:129338/1‑148 (MQ=255)
                                cACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTaa                                                                                                                      <  4:19460/149‑1 (MQ=255)
                                 aCCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAAc                                                                                                                     >  3:272936/1‑149 (MQ=255)
                                        aCGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTcag                                                                                                              <  4:55732/149‑1 (MQ=255)
                                          gCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTcagca                                                                                                            <  6:199798/149‑1 (MQ=255)
                                          gCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTcag                                                                                                              >  5:105965/1‑147 (MQ=255)
                                            tCGTGAAGTTGAGCGTAAGAAAGAGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGcagt                                                                                                          >  2:269018/1‑149 (MQ=255)
                                              gTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGcagtca                                                                                                        <  5:340884/149‑1 (MQ=255)
                                               tGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAAc                                                                                                                     <  8:326972/135‑1 (MQ=255)
                                               tGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAAc                                                                                                                     >  7:326972/1‑135 (MQ=255)
                                                gAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGcagtcagt                                                                                                      >  1:92966/1‑149 (MQ=255)
                                                       aGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTcagca                                                                                                            <  8:220896/136‑1 (MQ=255)
                                                       aGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTcagca                                                                                                            >  7:220896/1‑136 (MQ=255)
                                                        gCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTATCCTTTAACTAACCTGTAAGCATTCAGGTGCCGGtt                                                                                              >  2:77536/1‑149 (MQ=255)
                                                        gCGTAAGAAAGGGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCtt                                                                                              <  1:269018/149‑1 (MQ=255)
                                                           tAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTc                                                                                           >  6:136343/1‑149 (MQ=255)
                                                           tAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTc                                                                                           <  7:169655/149‑1 (MQ=255)
                                                            aaGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGACCCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCg                                                                                          >  7:110/1‑149 (MQ=255)
                                                              gAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGaagtcagt                                                                                                      >  6:329056/1‑135 (MQ=255)
                                                              gAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCa                                                                                        >  4:149890/1‑149 (MQ=255)
                                                                 aGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGCACCGGGC‑‑TTTGTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCa                                                                                        >  4:220662/1‑146 (MQ=255)
                                                                   tGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTAGAACTAACATGTCAGCAGTAAGTGGCCGCTTGTCGCATAGAc                                                                                   >  2:142382/1‑149 (MQ=255)
                                                                      gTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTAAGTGGCCGCTTGTCGCATAGACg                                                                                  >  2:118888/1‑147 (MQ=255)
                                                                       tCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAg                                                                               <  4:132152/149‑1 (MQ=255)
                                                                               aaGCGCGTAAGCGTCCTCATTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGTAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGc                                                                        <  1:77536/148‑1 (MQ=255)
                                                                                aGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCa                                                                       >  7:93034/1‑148 (MQ=255)
                                                                                aGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAAc                                                                                                                                                   <  1:67407/74‑1 (MQ=255)
                                                                                aGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAAc                                                                                                                                                   >  2:67407/1‑74 (MQ=255)
                                                                                   gcgTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGc                                                                   >  5:1710/1‑149 (MQ=255)
                                                                                     gTAAGCGTCCTCAGTACTCCAAGCGTTAATTCCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAg                                                                                     >  1:34614/1‑129 (MQ=255)
                                                                                     gTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAg                                                                                     <  2:34614/129‑1 (MQ=255)
                                                                                      gaaGCGTCCTCAGTACTCCAAGCGTTAATACCCCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCt                                                                                                                   <  6:40485/97‑1 (MQ=255)
                                                                                              cTCAGTACTCCTAGAGTAACTACCGCTTCGGGAGGCGAAAACAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTc                                                         <  8:203121/148‑1 (MQ=255)
                                                                                                 aGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGccc                                                     >  6:225225/1‑149 (MQ=255)
                                                                                                        cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCAgg                                              >  7:235251/1‑149 (MQ=255)
                                                                                                         aaGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                             <  1:142382/149‑1 (MQ=255)
                                                                                                          acccttAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                            <  8:56449/145‑1 (MQ=255)
                                                                                                          aGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                            <  5:136343/149‑1 (MQ=255)
                                                                                                              ttAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGt                                           >  8:134161/1‑146 (MQ=255)
                                                                                                                 aTACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTTAACTAACCTGTCAGCAGTAAGTGGCCGCTTGTAGCATAGAAGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCAc                                     >  3:88668/1‑149 (MQ=255)
                                                                                                                    cggCTTCGGCAGTCGAAAAACGCCCGCTTCCTTGGAACCGGGCTTTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCAcc                                    <  3:191992/147‑1 (MQ=255)
                                                                                                                           tgCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgtg                           <  8:110/148‑1 (MQ=255)
                                                                                                                           ggCAGTCGAAAAAAGCTCGGTTCCTTGGAACAGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgt                            <  2:110419/148‑1 (MQ=255)
                                                                                                                           ggCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgtg                           >  5:153140/1‑149 (MQ=255)
                                                                                                                           ggCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgt                            <  4:70410/148‑1 (MQ=255)
                                                                                                                              aGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgtgtgt                        <  3:149890/149‑1 (MQ=255)
                                                                                                                                tCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTgg                      >  5:5901/1‑149 (MQ=255)
                                                                                                                                  gAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGAt                    >  8:244939/1‑149 (MQ=255)
                                                                                                                                     aaaaGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTg                 >  6:154519/1‑149 (MQ=255)
                                                                                                                                      aatgCCCGGTATCTTTGATCCGGGC‑‑TTTTTTTTTAGCCTTTATCTAACCTGTCAGCAGTCAGTGGCCGCTTGTCTCATAGACGTAGTTAAAGCAAAGCGTTGGTTGTCTCCCTTTCCGTGGTAATCACCTTTTCAGTGTGTGGATTTGt                <  8:235251/146‑1 (MQ=255)
                                                                                                                                           ccGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCTGCTTGTTGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCtt           <  1:118888/149‑1 (MQ=255)
                                                                                                                                           ccGGGTCCAAGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAACCGTGGGTTGTCGCCGTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCt            >  3:217948/1‑148 (MQ=255)
                                                                                                                                                   ttGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGAAATGACGTAGATCAAGAAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACATTGTCAGTTTGTGGATTTGTTTCTTACCATGAc   >  3:90613/1‑149 (MQ=255)
                                                                                                                                                   ttGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGACGCTTGTCGAATAGACCTAGATCAAGCAAAGAGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTACCATGACg  >  4:263762/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGACG  >  NC_002947/1503920‑1504216

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: