Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,466,210 A→C 6.0% G13G (GGT→GGG)  PP_5559 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,466,2100AC 6.0% 190.8 / 2.2 62G13G (GGT→GGG) PP_5559hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/35);  new base C (2/2);  total (25/37)
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.56e-01

GGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCATTTGTTTTTGTACTGAC  >  NC_002947/3466066‑3466344
                                                                                                                                                |                                                                                                                                      
ggggTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGTCGGCAACGGCCAGCACAGGGCGCAGGCTGCGa                                                                                                                                        >  8:127256/1‑145 (MQ=255)
ggggTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGTCGGCAACGGCCAGCACAGGGCGCAGGCTGCGa                                                                                                                                        <  7:127256/145‑1 (MQ=255)
ggggTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCGGGCTGCGACCgg                                                                                                                                    >  8:24320/1‑149 (MQ=255)
ggggTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCgg                                                                                                                                    >  2:314182/1‑149 (MQ=255)
      tGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACCGGTCGCAGCCTGCGCCCtgtgctg                                                                                                                               <  8:73221/148‑4 (MQ=255)
       gAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGTACAGGGCGCAGGCTGCGACCGGTGTGGt                                                                                                                              <  6:247856/148‑1 (MQ=255)
       gAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGtt                                                                                                                             >  1:170962/1‑149 (MQ=255)
          tATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGc                                                                                                                          >  4:172569/1‑149 (MQ=255)
                 aTGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTgg                                                                                                                      <  3:172569/146‑1 (MQ=255)
                    ttCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCgg                                                                                                                <  7:85270/149‑1 (MQ=255)
                       ggcggcCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGgt                                                                                                              <  4:195164/148‑1 (MQ=255)
                       ggcggcCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGgt                                                                                                              <  2:331948/148‑1 (MQ=255)
                        gcggcCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGgtgt                                                                                                            >  2:55845/1‑149 (MQ=255)
                         cggcCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGgtgtg                                                                                                           <  4:63966/149‑1 (MQ=255)
                             cccTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCATGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTc                                                                                                       >  7:224034/1‑149 (MQ=255)
                             cccTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTc                                                                                                       <  2:48603/149‑1 (MQ=255)
                              ccTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCg                                                                                                      <  8:103982/149‑1 (MQ=255)
                               cTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCgg                                                                                                     >  6:331542/1‑149 (MQ=255)
                               cTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCgg                                                                                                     >  4:199743/1‑149 (MQ=255)
                                  gAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGc                                                                                                    <  5:72748/147‑1 (MQ=255)
                                    aGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGt                                                                                                <  1:11900/149‑1 (MQ=255)
                                         gACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGctc                                                                                           <  3:55206/149‑1 (MQ=255)
                                            gCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTcc                                                                                        <  2:296037/149‑1 (MQ=255)
                                                 gtgggtATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTg                                                                                    >  3:127230/1‑148 (MQ=255)
                                                  tgggtATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGaa                                                                                  >  3:6588/1‑149 (MQ=255)
                                                     gtATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGa                                                                                                                                        >  3:192192/1‑92 (MQ=255)
                                                     gtATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGa                                                                                                                                        <  4:192192/92‑1 (MQ=255)
                                                     gtATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGGGTTTCGGCATGTTGCTCTCCAGTGAAt                                                                                 >  4:327352/1‑147 (MQ=255)
                                                         cGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGccc                                                                           >  3:311408/1‑149 (MQ=255)
                                                             gcCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGc                                                                       <  4:311408/149‑1 (MQ=255)
                                                              ccTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCt                                                                      <  4:6588/149‑1 (MQ=255)
                                                                   tGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACtt                                                                 <  8:107047/149‑1 (MQ=255)
                                                                   tGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACtt                                                                 <  5:416240/149‑1 (MQ=255)
                                                                   tGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACtt                                                                 <  4:66465/149‑1 (MQ=255)
                                                                       tAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCCCGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCTGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGt                                                             >  2:312086/1‑149 (MQ=255)
                                                                       tAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCTGt                                                             >  6:203601/1‑149 (MQ=255)
                                                                         aaCCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGCCCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTc                                                           >  7:175063/1‑149 (MQ=255)
                                                                         aaCCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTc                                                           <  1:55845/149‑1 (MQ=255)
                                                                             aaTACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCGGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGc                                                       <  7:24320/149‑1 (MQ=255)
                                                                                 cccAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGc                                                   >  5:116671/1‑149 (MQ=255)
                                                                                  ccAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGt                                                             <  3:214103/138‑1 (MQ=255)
                                                                                  ccAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGt                                                             >  4:214103/1‑138 (MQ=255)
                                                                                         aaTAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCaa                                            <  1:268847/148‑1 (MQ=255)
                                                                                          aTAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAAt                                           <  5:94852/148‑1 (MQ=255)
                                                                                          aTAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATg                                          <  2:86820/149‑1 (MQ=255)
                                                                                                  cAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGt                                                                                                >  3:83996/1‑87 (MQ=255)
                                                                                                  cAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGt                                                                                                <  4:83996/87‑1 (MQ=255)
                                                                                                          aCGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTACGGGCCATTCAATGGATAAATGCTATGGAg                          >  6:157274/1‑149 (MQ=255)
                                                                                                                ccGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATc                    <  4:154105/149‑1 (MQ=255)
                                                                                                                    cAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTAGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCAttt                >  5:49424/1‑149 (MQ=255)
                                                                                                                       cGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATg                             <  3:328842/133‑1 (MQ=255)
                                                                                                                       cGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATg                             >  4:328842/1‑133 (MQ=255)
                                                                                                                       cGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCATTTGtt             <  5:331542/149‑1 (MQ=255)
                                                                                                                        ggCCAGCACAGGGCGCGGGCTGCGACAGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCATTTGttt            >  1:227502/1‑149 (MQ=255)
                                                                                                                        ggCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCATTTGttt            <  1:42570/149‑1 (MQ=255)
                                                                                                                         gCCAGCACAGGGCGCATGCTGCGACCGGTGTGGTTGGCGTGGTTGCTGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCATTTGtttt           <  1:5328/149‑1 (MQ=255)
                                                                                                                         gCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCa                                                                                       >  1:101202/1‑73 (MQ=255)
                                                                                                                         gCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCa                                                                                       <  2:101202/73‑1 (MQ=255)
                                                                                                                         gCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCATTTGtttt           <  3:199743/149‑1 (MQ=255)
                                                                                                                          ccAGCACAGGGCGCAGGCTGCGACGGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCATTTGttttt          <  1:312086/149‑1 (MQ=255)
                                                                                                                          ccAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTCCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCATTTGttttt          <  1:314182/149‑1 (MQ=255)
                                                                                                                             gCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCATTTGTTTTTGTa       <  4:127230/149‑1 (MQ=255)
                                                                                                                             acaccGGTCGCAGCCTGCGCCCTGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCATTTGtttagaat       >  5:109244/2‑144 (MQ=255)
                                                                                                                             acaccGGTCGCAGCCTGCGCCCTGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATCAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCATTTGtt             <  6:109244/142‑1 (MQ=255)
                                                                                                                                 aGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCCGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCATTTGTTTTTGTactga   >  8:50525/1‑149 (MQ=255)
                                                                                                                                  gggCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCATTTGTTTTTGTactgac  <  5:412765/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                      
GGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTGCGACCGGTGTGGTTGGCGTGGTTGCGGGTGTGTTTCGGCATGTTGCTCTCCAGTGAATGAGCCCGAGCTTACTTCAGTTCGTGCGGGCCATTCAATGGATAAATGCTATGGAGCCGATCATTTGTTTTTGTACTGAC  >  NC_002947/3466066‑3466344

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