Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,659,186 T→A 12.7% L13Q (CTG→CAG)  nuoE → NADH‑quinone oxidoreductasesubunit E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,659,1860TA12.7% 153.5 / 14.2 63L13Q (CTG→CAG) nuoENADH‑quinone oxidoreductasesubunit E
Reads supporting (aligned to +/- strand):  ref base T (30/25);  new base A (4/4);  total (34/29)
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.72e-01

TACCCGCACCTGCAGCAGATCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGCCATCGGCGAAGTGCTGGG  >  NC_002947/4659039‑4659333
                                                                                                                                                   |                                                                                                                                                     
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                                                                                                                                                     <  6:92646/131‑1 (MQ=255)
                   tCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGa                                                                                                                                                     >  5:92646/1‑131 (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                                                                            <  4:292956/149‑1 (MQ=255)
                                                                        aGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGaa                                                                            <  5:234837/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                                            >  3:12763/1‑149 (MQ=255)
                                                                                                        cTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTg                                            >  8:115209/1‑149 (MQ=255)
                                                                                                        cTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTg                                            >  2:213155/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                           <  2:220432/149‑1 (MQ=255)
                                                                                                                         gCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCa                           <  7:115209/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                                      >  2:229875/6‑126 (MQ=255)
                                                                                                                                     gctcggTTTCGCTCAGGGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGcc                                      <  1:229875/121‑1 (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)
                                                                                                                                                   tGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGCCATCGGCGAAGTGCTggg  >  2:55584/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TACCCGCACCTGCAGCAGATCCCTTCGGTGATCAAAGGCAGCATGGTTGCCGACCTCATTGCGTACCTGGGCAGTATCGACTTCGTTATGGCTGACGTGGACCGCTAAGCATGAACAGCACGCTTATCCAGACAGACCGTTTCGCCCTGAGCGAAACCGAGCGC‑‑TCGGCCATCGAGCACGAAATGCATCACTACGAGGACCCGCGCGCGGCGTCCATCGAAGCCCTGAAGATCGTCCAGAAGGAGCGTGGCTGGGTGCCGGACGGCGCCATTCACGCCATCGGCGAAGTGCTGGG  >  NC_002947/4659039‑4659333

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: