Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,668,407 G→A 5.7% L93L (CTG→CTA)  nuoM → NADH‑quinone oxidoreductase subunit M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,668,4070GA 5.7% 207.1 / 2.1 70L93L (CTG→CTA) nuoMNADH‑quinone oxidoreductase subunit M
Reads supporting (aligned to +/- strand):  ref base G (29/37);  new base A (2/2);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

GTGGCTGTGGGGTACCGGCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGACCTGTTCCTGTTCTTC  >  NC_002947/4668260‑4668551
                                                                                                                                                   |                                                                                                                                                
ctggctgtggGGTACCTGCGACTACACCGTTGCTCCCGCCCCGGGCTCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTg                                                                                                                                                  <  6:355885/147‑1 (MQ=255)
   gctgtggGGTACCGGCGACTACACACTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGACTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTCGACGGCCTGTCCCTGCTCATGATCCTg                                                                                                                                                  >  5:355885/1‑145 (MQ=255)
     tgtggGGTACCGGCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCAcc                                                                                                                                            >  6:91000/1‑149 (MQ=255)
       tggGGTACCGGCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCg                                                                                                                                           <  5:116718/148‑1 (MQ=255)
        ggggTACCGGCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGc                                                                                                                                         <  4:330158/149‑1 (MQ=255)
         gggTACCGGCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGcc                                                                                                                                        >  7:78654/1‑149 (MQ=255)
           gTACCGGCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTg                                                                                                                                      >  1:351068/1‑149 (MQ=255)
               cGGCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCg                                                                                                                                  <  4:170159/149‑1 (MQ=255)
               cGGCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCg                                                                                                                                  <  8:310474/149‑1 (MQ=255)
                 gCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGgt                                                                                                                                <  6:293723/149‑1 (MQ=255)
                  cGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGgtg                                                                                                                               <  2:351068/149‑1 (MQ=255)
                   gACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGG‑CCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGtgtg                                                                                                                             <  8:30098/149‑1 (MQ=255)
                        cacCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGt                                                                                                                         >  5:218159/1‑149 (MQ=255)
                        cacCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGt                                                                                                                         >  5:245144/1‑149 (MQ=255)
                        cacCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCGCACCGGCCTGCTCGGTGTGCTGt                                                                                                                         <  7:77907/149‑1 (MQ=255)
                            cTTGCTCCCGCCCCGGTCGCCTAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGt                                                                                                                     <  6:209690/149‑1 (MQ=255)
                                cTCCCGCTCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACt                                                                                                                  <  1:82502/148‑1 (MQ=255)
                                cTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGACTACAAAGTCGAGTGGATCAAGCGTTTCGGCATCATCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACt                                                                                                                  <  5:142878/148‑1 (MQ=255)
                                 tCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACtgt                                                                                                                >  7:79612/1‑149 (MQ=255)
                                 tCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACtgt                                                                                                                >  7:260255/1‑149 (MQ=255)
                                   ccGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGGGGATCAAGCGTGTCGGCATCAGCATCCCCAAGGCCCTGGACGGACTGTCGATGATGATGATCCTGCACAACGAGCAGCCCGgtgacaa                                                                                                                           >  3:214458/1‑132 (MQ=255)
                                        ccGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTg                                                                                                         >  5:292366/1‑149 (MQ=255)
                                            gcgcCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGaaa                                                                                                     <  6:292366/149‑1 (MQ=255)
                                            gcgcCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGaaa                                                                                                     <  2:10823/149‑1 (MQ=255)
                                              gcCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAaga                                                                                                   <  5:91000/149‑1 (MQ=255)
                                                            ggCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCAc                                                                                      <  8:260255/148‑1 (MQ=255)
                                                                 tcGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGaaa                                                                                                     >  3:344996/1‑128 (MQ=255)
                                                                 tcGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGaaa                                                                                                     <  4:344996/128‑1 (MQ=255)
                                                                 tcGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACCGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCAGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGc                                                                                >  2:328989/1‑149 (MQ=255)
                                                                   gAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGctt                                                                              <  3:208841/149‑1 (MQ=255)
                                                                     aTACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGcttct                                                                            <  8:79612/149‑1 (MQ=255)
                                                                            gTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCt                                                                     <  8:78654/149‑1 (MQ=255)
                                                                            gTCGAGTGGATCAAGCGTTTCGGCATCACCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCt                                                                     >  1:201028/1‑149 (MQ=255)
                                                                               gAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCaa                                                                  >  4:237590/1‑149 (MQ=255)
                                                                                    gttCAAGCGTTTGGGCATCAGCATCCTCGTGGGCCTGGACGGCGTGTCGCTGCTGATGTTCCTGCTCACCGGCCTGCTCGTTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCttcttc                                                                          <  3:278741/134‑1 (MQ=255)
                                                                                    gATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGa                                                             >  8:77961/1‑149 (MQ=255)
                                                                                    gATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCACGGAGTGCTGACGGAACTGTGTTCCTGGAAAGAGATAAAGCGCCACGTCCGCttcatc                                                                          >  4:278741/1‑136 (MQ=255)
                                                                                      tCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTACTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGAtg                                                           >  2:12522/1‑149 (MQ=255)
                                                                                      tCAAGCGTTTCGGCATCAGCATACACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCATGCTCGGTGTGCTGTCGGTAATGTGTTCCTGGAAAGAGATCCAGCGCCACGTCCGATTCTTCCACCTCAACCTGaag                                                           >  2:169635/1‑147 (MQ=255)
                                                                                         aGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGAtgtg                                                         >  4:93117/1‑148 (MQ=255)
                                                                                         aGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTgg                                                        <  2:201028/149‑1 (MQ=255)
                                                                                          gCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGa                                                       <  7:77961/149‑1 (MQ=255)
                                                                                             tttCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATcc                                                    >  2:160758/1‑149 (MQ=255)
                                                                                                cGGCAACAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTc                                                  <  3:114823/148‑1 (MQ=255)
                                                                                                 ggCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTcg                                                 <  3:237590/148‑1 (MQ=255)
                                                                                                     tCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCggcggc                                            >  6:281442/1‑149 (MQ=255)
                                                                                                        gCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGt                                          <  7:292377/148‑1 (MQ=255)
                                                                                                        gCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGt                                          >  8:292377/1‑148 (MQ=255)
                                                                                                        gCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTg                                         >  8:220149/1‑149 (MQ=255)
                                                                                                            ccACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCg                                     >  5:310498/1‑149 (MQ=255)
                                                                                                              aCCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGgt                                   <  3:93117/149‑1 (MQ=255)
                                                                                                              aCCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGgt                                   <  7:220149/149‑1 (MQ=255)
                                                                                                              aCCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGgt                                   <  5:281442/149‑1 (MQ=255)
                                                                                                                cTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGgtgt                                 <  6:245144/149‑1 (MQ=255)
                                                                                                                   gCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACCGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTCCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGtt                              <  6:218159/149‑1 (MQ=255)
                                                                                                                     ccTCGACGGCCTGACGCTGCTCATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTCTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTCGATCCTCGGCCGCGTCGTCGGTGTGTTcc                            >  8:130600/1‑149 (MQ=255)
                                                                                                                      cTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCt                           >  8:197166/1‑149 (MQ=255)
                                                                                                                          aCGGCCTGTCGCTGCTGATGATCCTACTCACCGGCCTGCTCGGTGTGCGGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGcc                       <  1:12522/149‑1 (MQ=255)
                                                                                                                                  tCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGAcct               >  4:15972/1‑149 (MQ=255)
                                                                                                                                  tCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGAc                 <  1:160758/147‑1 (MQ=255)
                                                                                                                                   cGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGAcctg              <  7:197166/149‑1 (MQ=255)
                                                                                                                                    gctgctGATGATCCTACTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTc                             <  3:105467/133‑1 (MQ=255)
                                                                                                                                    gctgatGATGATCCTACTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGATTCTTCCAACTCAACCTGATGTGGATCCTCGCCCGCTTGGTCGGTGTGTTc                             >  4:105467/1‑133 (MQ=255)
                                                                                                                                      tgctgATGATCCTGCTCACCGGCCTGGTCGGTGTGCTGTCGGTACTGTGTTCCTGGATAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCATCCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGAcctgttc           <  7:130600/149‑1 (MQ=255)
                                                                                                                                        ctgATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGAcctgttcct         >  7:145014/1‑149 (MQ=255)
                                                                                                                                           atgatCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGAcctgttcctgt       <  7:4001/148‑1 (MQ=255)
                                                                                                                                           atgatCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGAcctgttcctgt       >  8:294768/1‑148 (MQ=255)
                                                                                                                                             gatCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTgg                   <  5:83273/134‑1 (MQ=255)
                                                                                                                                             gatCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTgg                   >  6:83273/1‑134 (MQ=255)
                                                                                                                                               tCCTGCTCAGCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGTTCTGTGTGTTCCTGGCCCTGGACCTGTTCCTGttcttc  <  1:169635/149‑1 (MQ=255)
                                                                                                                                               tCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGACCTGTTCCTGttct    >  6:151919/1‑147 (MQ=255)
                                                                                                                                               tCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGACCTGTTCCTGttct    <  5:151919/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GTGGCTGTGGGGTACCGGCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCTTCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGACCTGTTCCTGTTCTTC  >  NC_002947/4668260‑4668551

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