Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,383,070 T→G 5.7% E77D (GAA→GAC)  PP_1204 ← lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,383,0700TG 5.7% 219.9 / 4.6 70E77D (GAA→GAC) PP_1204lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (37/29);  new base G (2/2);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

TTTGCGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGG  >  NC_002947/1382925‑1383212
                                                                                                                                                 |                                                                                                                                              
tttGCGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCg                                                                                                                                             >  3:231800/1‑149 (MQ=255)
              cccAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAgcg                                                                                                                               >  8:159378/1‑149 (MQ=255)
               ccAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAgcgc                                                                                                                              >  2:108106/1‑149 (MQ=255)
                cAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCt                                                                                                                             >  1:80900/1‑149 (MQ=255)
                cAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAAGTCACCGAATTCGAAgtgac                                                                                                                             >  3:278219/1‑145 (MQ=255)
                  aGGCGAAGACGTTGACCACACTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGATGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCttt                                                                                                                           >  2:180595/1‑149 (MQ=255)
                   ggCGACGACGTTGACCACATTCGCGCCGACACCCCCCACCACGGCCACTGCACCCAGCGGGATCTCCTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCtttt                                                                                                                          <  8:243742/149‑1 (MQ=255)
                   ggCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAAGTCACCGAATTCGAAgtgacccg                                                                                                                          <  4:278219/149‑8 (MQ=255)
                                  cacaTTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGca                                                                                                           >  7:88751/1‑149 (MQ=255)
                                    caTTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGt                                                                                                                                           <  1:324467/115‑1 (MQ=255)
                                    caTTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGt                                                                                                                                           >  2:324467/1‑115 (MQ=255)
                                    caTTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGcaca                                                                                                         <  5:85605/149‑1 (MQ=255)
                                         gcgcCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCAcgcg                                                                                                             <  1:58234/140‑1 (MQ=255)
                                         gcgcCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCAcgcg                                                                                                             >  2:58234/1‑140 (MQ=255)
                                                   ccGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCa                                                                                                                     >  5:122751/1‑122 (MQ=255)
                                                   ccGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCa                                                                                                                     <  6:122751/122‑1 (MQ=255)
                                                   ccGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGCGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAg                                                                                           >  6:277755/1‑148 (MQ=255)
                                                     gCCACCACGGCCACTGCACCCAGCTGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTc                                                                                        <  8:88751/149‑1 (MQ=255)
                                                             ggCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCt                                                                                <  2:236247/149‑1 (MQ=255)
                                                             ggCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCt                                                                                <  1:108106/149‑1 (MQ=255)
                                                               ccACTGCACCCAGCGGTATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTc                                                                              <  3:30967/149‑1 (MQ=255)
                                                               ccACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTc                                                                              >  5:243411/1‑149 (MQ=255)
                                                               ccACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTc                                                                              >  4:187485/1‑149 (MQ=255)
                                                               ccACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTc                                                                              >  1:19975/1‑149 (MQ=255)
                                                                cACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCg                                                                             >  1:292700/1‑149 (MQ=255)
                                                                 aCTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTg                                                                                                                                          >  7:238943/1‑87 (MQ=255)
                                                                 aCTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTg                                                                                                                                          <  8:238943/87‑1 (MQ=255)
                                                                         cAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGa                                                                    <  2:80900/149‑1 (MQ=255)
                                                                                aTCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTc                                                                              <  6:216653/132‑1 (MQ=255)
                                                                                aTCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTc                                                                              >  5:216653/1‑132 (MQ=255)
                                                                                   tcGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTg                                                          >  5:332289/1‑149 (MQ=255)
                                                                                     gTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTg                                                          >  5:172996/1‑147 (MQ=255)
                                                                                     gTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTg                                                          >  3:209436/1‑147 (MQ=255)
                                                                                      tAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAg                                                       <  3:187485/149‑1 (MQ=255)
                                                                                        gATCTCCCGCCAGGGTTGGTAAGGCTTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAgg                                                      <  5:277755/148‑1 (MQ=255)
                                                                                        gATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAgg                                                      <  6:332289/148‑1 (MQ=255)
                                                                                         aTCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTc                                                                              >  4:89983/1‑123 (MQ=255)
                                                                                         aTCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTc                                                                              <  3:89983/123‑1 (MQ=255)
                                                                                            tcCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATg                                                 >  2:233328/1‑149 (MQ=255)
                                                                                            tcCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATg                                                 >  7:289830/1‑149 (MQ=255)
                                                                                             cccGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGc                                                >  5:23013/1‑149 (MQ=255)
                                                                                               ccgctgggTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTc                                              >  7:256113/6‑149 (MQ=255)
                                                                                                 ccAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTc                                              >  7:39179/1‑147 (MQ=255)
                                                                                                    gggTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCTGCTCGACGACATCGATCTGCAGGGTATGCTCGAg                                           >  6:91350/1‑147 (MQ=255)
                                                                                                    gggTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCa                                         <  5:91350/149‑1 (MQ=255)
                                                                                                          gTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTgc                                   >  7:214108/1‑149 (MQ=255)
                                                                                                          gTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTgc                                   >  7:244949/1‑149 (MQ=255)
                                                                                                           tAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTgcg                                  <  6:243411/149‑1 (MQ=255)
                                                                                                            aaGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCa                                         <  3:218974/141‑1 (MQ=255)
                                                                                                            aaGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCa                                         >  4:218974/1‑141 (MQ=255)
                                                                                                                 gtgtAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGAc                            <  8:289830/149‑1 (MQ=255)
                                                                                                                   gtAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGAcgc                          >  7:47150/1‑149 (MQ=255)
                                                                                                                    tAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGAcgcg                         >  1:119819/1‑149 (MQ=255)
                                                                                                                         gTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGctcct                    >  4:230501/1‑149 (MQ=255)
                                                                                                                            aTAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAg                                                       <  5:73467/111‑1 (MQ=255)
                                                                                                                            aTAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAg                                                       >  6:73467/1‑111 (MQ=255)
                                                                                                                            aTAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCAt                 <  6:23013/149‑1 (MQ=255)
                                                                                                                             tAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATg                >  8:276639/1‑149 (MQ=255)
                                                                                                                               acgacgGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGCCGCGCTCCTCATgg               <  7:276639/148‑1 (MQ=255)
                                                                                                                               acgacgGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGt              <  8:256113/149‑1 (MQ=255)
                                                                                                                                   aaaagcgcTTCGAAGTCACCGAATTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTcgct          >  6:268599/8‑149 (MQ=255)
                                                                                                                                     aagcgcTTCGAAGTCACCGAATTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTcgctgc        <  5:268599/144‑1 (MQ=255)
                                                                                                                                          cTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTgg   <  8:214108/149‑1 (MQ=255)
                                                                                                                                          cTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTgg   <  2:119819/149‑1 (MQ=255)
                                                                                                                                          cTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTgg   <  6:172996/149‑1 (MQ=255)
                                                                                                                                          cTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTgg   >  7:129091/1‑149 (MQ=255)
                                                                                                                                           ttCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGCCGCGCTCCTCATGGTCGCTGCGCTggg  >  6:90426/1‑149 (MQ=255)
                                                                                                                                             cGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTggg  <  1:193997/147‑1 (MQ=255)
                                                                                                                                                aTTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTgg   >  8:202779/1‑143 (MQ=255)
                                                                                                                                                aTTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTgg   <  7:202779/143‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TTTGCGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGG  >  NC_002947/1382925‑1383212

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