Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,659,184 C→T 12.9% A12A (GCC→GCT nuoE → NADH‑quinone oxidoreductasesubunit E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,659,1840CT12.9% 154.2 / 12.3 62A12A (GCC→GCTnuoENADH‑quinone oxidoreductasesubunit E
Reads supporting (aligned to +/- strand):  ref base C (32/22);  new base T (4/4);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-01

TACCCGCACCTGCAGCAGATCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGCCATCGGCGAAGTGCTGG  >  NC_002947/4659039‑4659332
                                                                                                                                                 |                                                                                                                                                      
tACCCGCACCTGCAGCAGATCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCt                                                                                                                                                      >  3:170829/1‑148 (MQ=255)
 aCCCGCACCTGCAGCAGATCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTCCCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCAGGCTTATCGAGCCAGACCGTGTCGCCGTgg                                                                                                                                                    >  8:68437/1‑148 (MQ=255)
 aCCCGCACCTGCAGCAGATCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCt                                                                                                                                                      <  7:68437/147‑1 (MQ=255)
 aCCCGCACCTGCAGCAGATCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCt                                                                                                                                                      >  1:86084/1‑147 (MQ=255)
   ccGCACCTGCAGCAGATCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGc                                                                                                                                                  >  4:271139/1‑149 (MQ=255)
        ccTTCAGCAGATCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAAc                                                                                                                                             <  6:342924/149‑1 (MQ=255)
         cTGCAGCAGATCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCTCAGGGCGAAAcg                                                                                                                                            <  7:60735/149‑2 (MQ=255)
         cTGCAGCAGATCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAAcc                                                                                                                                            <  3:271139/149‑1 (MQ=255)
                   tCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGa                                                                                                                                                    >  5:92646/1‑131 (MQ=255)
                   tCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGa                                                                                                                                                    <  6:92646/131‑1 (MQ=255)
                      cTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCa                                                                                                                             >  3:194091/1‑149 (MQ=255)
                              aTCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCTCAGGGCGAAAcggtc                                                                                                                                         <  8:130848/131‑5 (MQ=255)
                              aTCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCTCAGGGCGAAAcggtc                                                                                                                                         >  7:130848/1‑127 (MQ=255)
                                 aaaGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGa                                                                                                                   <  4:194091/148‑1 (MQ=255)
                                  aaGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGaaa                                                                                                                 >  7:148779/1‑149 (MQ=255)
                                     gcagcaTGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATg                                                                                                               >  6:234837/1‑148 (MQ=255)
                                     gcagcaTGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATg                                                                                                               >  2:268696/1‑148 (MQ=255)
                                     gcagcaTGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGc                                                                                                              <  1:157642/149‑1 (MQ=255)
                                         caTGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCa                                                                                                          <  6:82784/149‑1 (MQ=255)
                                                gCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAg                                                                                                   >  7:101007/1‑149 (MQ=255)
                                                  cGACCTCATTGCGTACCTGGGCAGTATCTACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGATTCGCCCTGAGCGAACCCGAGCGC‑‑TCGGCCCTCGCTCAGGAAATTCATCACTACGAGGa                                                                                                 >  4:233346/1‑149 (MQ=255)
                                                        cATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCcgcg                                                                                           >  1:200800/1‑149 (MQ=255)
                                                            gCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCcgcgcgc                                                                                        <  5:174573/148‑1 (MQ=255)
                                                                 ccTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGTCCCGCGCGCGGCGTc                                                                                  <  2:86084/149‑1 (MQ=255)
                                                                        aGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGaa                                                                           <  5:234837/149‑1 (MQ=255)
                                                                        aGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGaa                                                                           <  4:292956/149‑1 (MQ=255)
                                                                          tATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGc                                                                         >  8:48878/1‑149 (MQ=255)
                                                                          tATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCAC‑‑TCGGCCATCGAACACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAga                                                                         >  7:342815/1‑148 (MQ=255)
                                                                                    gTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGaa                                                                >  1:128855/1‑147 (MQ=255)
                                                                                              aCGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGaa                                                       <  8:101007/147‑1 (MQ=255)
                                                                                                gTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACt                                                                                                        <  8:85155/96‑1 (MQ=255)
                                                                                                gTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCAACACt                                                                                                        >  7:85155/1‑96 (MQ=255)
                                                                                                 tGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGc                                                  >  2:235433/1‑149 (MQ=255)
                                                                                                     ccGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTgg                                              <  8:148779/149‑1 (MQ=255)
                                                                                                      cGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGc                                             >  1:220432/1‑149 (MQ=255)
                                                                                                       gCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCt                                            <  3:56923/149‑1 (MQ=255)
                                                                                                        cTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTg                                           >  8:115209/1‑149 (MQ=255)
                                                                                                        cTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTg                                           >  2:213155/1‑149 (MQ=255)
                                                                                                        cTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTg                                           >  3:12763/1‑149 (MQ=255)
                                                                                                         tAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTgg                                          >  3:57476/1‑149 (MQ=255)
                                                                                                         tAAGCATGAACAGCAAGCTGATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTgg                                          >  2:5099/1‑149 (MQ=255)
                                                                                                             cATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGCGGACGGCCAACGAGCACGAAATGAATCACTACGAGGACCCGCCCGCGGCGTCAATCGAAGCCCTGAACATCTTCAAGAAGGAGCGTGGATggtt                                        >  1:98726/1‑147 (MQ=255)
                                                                                                                 aaCAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGa                                  >  3:220481/1‑149 (MQ=255)
                                                                                                                      cACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGgcg                             >  6:333072/1‑149 (MQ=255)
                                                                                                                         gCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCa                          <  7:115209/149‑1 (MQ=255)
                                                                                                                         gCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCa                          <  2:220432/149‑1 (MQ=255)
                                                                                                                             aTCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCa                      >  6:292547/1‑149 (MQ=255)
                                                                                                                             aTCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCa                      >  6:292549/1‑149 (MQ=255)
                                                                                                                             aTCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCAGCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTc                       <  3:233346/148‑1 (MQ=255)
                                                                                                                               ccAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACg                    >  2:351899/1‑149 (MQ=255)
                                                                                                                                 agtcagacCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGc                   >  3:232079/4‑148 (MQ=255)
                                                                                                                                 agacagacCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACgg                               <  4:28864/136‑1 (MQ=255)
                                                                                                                                 agacagacCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACgg                               >  3:28864/1‑136 (MQ=255)
                                                                                                                                     gctcggTTTCGCTCAGGGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGcc                                     <  1:229875/121‑1 (MQ=255)
                                                                                                                                     gctcggTTTCGCTCAGGGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGcc                                     >  2:229875/6‑126 (MQ=255)
                                                                                                                                     agacAGTTTCGCCCTGAGCGAAAACGAGCGC‑‑TCGGCCAACGAGCACGAAATGCATCACTACGAGGACCCCCGCGCGGCGTACATCGAAGCCCAGAAGATAGTCCAGAAGGAGCGATTCTCGGTGCCGGAAGGAGCCATTCACGCCATcc              >  4:247261/1‑148 (MQ=255)
                                                                                                                                      gacCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGCCATCgg             <  1:351899/149‑1 (MQ=255)
                                                                                                                                         ggTTTCGCTCAGGGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGCCATCGGCg           >  7:200628/2‑148 (MQ=255)
                                                                                                                                         cGTTTCGCCATGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCAGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGCCATCGGCGa          <  3:247261/149‑1 (MQ=255)
                                                                                                                                          gTTTCGCTCAGGGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGCCATCGGCg           >  3:100289/1‑147 (MQ=255)
                                                                                                                                          gTTTCGCTCAGGGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGCCATCGGCGaa         <  8:200628/149‑1 (MQ=255)
                                                                                                                                          gTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGCCATCGGCGaa         <  7:48878/149‑1 (MQ=255)
                                                                                                                                                cccTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGTCTGGGTGCCGGACGGCGCCATTCACGCCATCGGCGAAGTGCTg   <  2:200800/149‑1 (MQ=255)
                                                                                                                                                 ccTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGCCATCGGCGAAGTGCttg  <  5:292549/149‑3 (MQ=255)
                                                                                                                                                 ccTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGCCATCGGCGAAGTGCTgg  <  5:292547/149‑1 (MQ=255)
                                                                                                                                                 ccTGAGCGAAACCGAGCAC‑‑TCGGCCAGCGCACACGAAATGCAGCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGCCATCGGCGAAGTGCTgg  <  8:342815/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
TACCCGCACCTGCAGCAGATCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGCCATCGGCGAAGTGCTGG  >  NC_002947/4659039‑4659332

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: