Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,874,875 G→T 5.8% S451* (TCG→TAG) ‡ ilvA‑II ← threonine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,874,8750GT 5.8% 275.0 / 3.4 85S451* (TCG→TAG) ‡ilvA‑IIthreonine deaminase
Reads supporting (aligned to +/- strand):  ref base G (34/46);  new base T (2/3);  total (36/49)
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.48e-01

TTTGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGTT  >  NC_002947/5874727‑5875023
                                                                                                                                                    |                                                                                                                                                    
tttGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCga                                                                                                                                                     <  2:1456096/150‑1 (MQ=255)
tttGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCga                                                                                                                                                     <  1:1228287/150‑1 (MQ=255)
tttGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCga                                                                                                                                                     <  2:1456572/150‑1 (MQ=255)
tttGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCga                                                                                                                                                     >  2:1114514/1‑150 (MQ=255)
  tGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCCCTAGGGTTGCGGTAGTGGAACATCgaga                                                                                                                                                   <  1:1399994/150‑1 (MQ=255)
   gTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGTGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGAt                                                                                                                                                  >  2:401280/1‑150 (MQ=255)
   gTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGAt                                                                                                                                                  <  2:210380/150‑1 (MQ=255)
       gCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAGCCGATCTTGGCCAGCGCTGCCGGCACCAGGGGACGTGCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTc                                                                                                                                              >  2:1770269/1‑150 (MQ=255)
       gCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTGGGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTc                                                                                                                                              >  1:959005/1‑150 (MQ=255)
        ccGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTcc                                                                                                                                             >  2:1570161/1‑150 (MQ=255)
         cGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCa                                                                                                                                            <  1:1133292/150‑1 (MQ=255)
         cGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCa                                                                                                                                            <  1:1289701/150‑1 (MQ=255)
          gggTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAg                                                                                                                                           >  1:1759383/1‑150 (MQ=255)
            gTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCTCGATGTTCCacta                                                                                                                                         <  1:279063/150‑4 (MQ=255)
            gTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGACCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCTCGATGTTCCacta                                                                                                                                         <  2:1696727/150‑4 (MQ=255)
                       cGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCt                                                                                                                              <  1:1114514/150‑1 (MQ=255)
                        gTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCtt                                                                                                                             >  2:1793504/1‑150 (MQ=255)
                        gTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCtt                                                                                                                             >  1:341824/1‑150 (MQ=255)
                        gTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCtt                                                                                                                             >  2:1760314/1‑150 (MQ=255)
                           ccAGTACGTGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCttgtt                                                                                                                          <  1:401280/150‑1 (MQ=255)
                           ccAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCttgtt                                                                                                                          <  2:1796514/150‑1 (MQ=255)
                            cAGTACGGGTAACCGATCTTGGCCAGCGCTGACAACACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCttgttg                                                                                                                         >  1:1400832/1‑150 (MQ=255)
                             aGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGa                                                                                                                        <  1:1512122/150‑1 (MQ=255)
                              tttcGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAg                                                                                                                       <  1:967222/147‑1 (MQ=255)
                                  gggTAACCGATCTTGGCCAGCGCTGCCGGCAAACCGCTACTGGCTTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAg                                                                                                                   <  2:1400832/150‑1 (MQ=255)
                                   ggTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGt                                                                                                                  <  1:868033/150‑1 (MQ=255)
                                       aCCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGaa                                                                                                              <  1:174193/150‑1 (MQ=255)
                                       aCCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGaa                                                                                                              <  1:1563463/150‑1 (MQ=255)
                                       aCCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGaa                                                                                                              <  2:1668717/150‑1 (MQ=255)
                                       aCCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGaa                                                                                                              <  1:1133488/150‑1 (MQ=255)
                                        ccGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTACAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCga                                                                                                                                                     <  2:1157662/110‑1 (MQ=255)
                                        ccGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTACAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCga                                                                                                                                                     >  1:1157662/1‑110 (MQ=255)
                                         cGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACa                                                                                                            <  1:5229/150‑1 (MQ=255)
                                         cGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACa                                                                                                            <  2:1759383/150‑1 (MQ=255)
                                         cGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACa                                                                                                            >  1:136684/1‑150 (MQ=255)
                                             cTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGc                                                                                                        <  1:178543/150‑1 (MQ=255)
                                                   cAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTgg                                                                                                  <  1:1229781/150‑1 (MQ=255)
                                                           ccGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCATGACCGCTCgg                                                                                          >  2:433432/1‑150 (MQ=255)
                                                               cACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGaa                                                                                      <  2:277358/150‑1 (MQ=255)
                                                                 ccAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAAtt                                                                                    <  2:1871986/150‑1 (MQ=255)
                                                                  cAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCggg                                                                                         >  2:1688076/1‑144 (MQ=255)
                                                                  cAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCggg                                                                                         <  1:1688076/144‑1 (MQ=255)
                                                                     gTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGaa                                                                                <  2:1707737/150‑1 (MQ=255)
                                                                      tGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAg                                                                               >  1:1102536/1‑150 (MQ=255)
                                                                      tGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAg                                                                               <  1:1760314/150‑1 (MQ=255)
                                                                      tGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAg                                                                               >  1:239191/1‑150 (MQ=255)
                                                                          gTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCgca                                                                           >  1:1830907/1‑150 (MQ=255)
                                                                           ttCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCgcag                                                                          >  1:363476/1‑150 (MQ=255)
                                                                               tcctccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCAcc                                                                      <  2:239191/150‑1 (MQ=255)
                                                                                 ctccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCAt                                                                    >  1:788998/1‑150 (MQ=255)
                                                                                  tccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGCGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACTAtt                                                                   >  1:644858/1‑150 (MQ=255)
                                                                                  tccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCAtt                                                                   >  2:1292535/1‑150 (MQ=255)
                                                                                        aCCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTtcgtcg                                                             >  1:1211947/1‑150 (MQ=255)
                                                                                           tGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTg                                                          <  1:899738/150‑1 (MQ=255)
                                                                                           tGCAGCCCCGACACCACGCGCCCGTCAGTAAAGATATATACGCGGTAGGGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATCCGAAGCGCAGCACCATTTCGTCGCTg                                                          <  2:959005/150‑1 (MQ=255)
                                                                                               gCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCAc                                                      >  2:1819877/1‑150 (MQ=255)
                                                                                                 cccGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGAGCAAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCg                                                    <  1:495806/150‑1 (MQ=255)
                                                                                                  ccGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCgg                                                   >  1:357557/1‑150 (MQ=255)
                                                                                                   cGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGTAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGc                                                  >  2:309087/1‑150 (MQ=255)
                                                                                                      caccacGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGTTGGCACCGGCCGa                                               <  2:594575/150‑1 (MQ=255)
                                                                                                      caccacGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGa                                               <  2:341824/150‑1 (MQ=255)
                                                                                                         cacGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTg                                            <  1:653073/150‑1 (MQ=255)
                                                                                                           cgcgCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGgc                                          <  2:102146/150‑1 (MQ=255)
                                                                                                            gcgcCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGcgcccgc                                                                                                                                    <  2:102902/59‑1 (MQ=255)
                                                                                                            gcgcCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGcgcccgc                                                                                                                                    >  1:102902/1‑59 (MQ=255)
                                                                                                            gcgcCCGTCAGCCGCGCCATGGTTGCGGTAGTAGAACATAGAGATGTTCCAGCGCCCGCCCACCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGgcc                                         >  2:690123/1‑150 (MQ=255)
                                                                                                                 cGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGc                                                                                                 <  2:1348197/89‑1 (MQ=255)
                                                                                                                 cGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGc                                                                                                 <  2:1348118/89‑1 (MQ=255)
                                                                                                                 cGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGc                                                                                                 >  1:1348118/1‑89 (MQ=255)
                                                                                                                 cGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGc                                                                                                 >  1:1348197/1‑89 (MQ=255)
                                                                                                                 cGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGa                                    <  2:357557/150‑1 (MQ=255)
                                                                                                                  gTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACa                                                                                                            <  2:1227615/77‑1 (MQ=255)
                                                                                                                  gTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACa                                                                                                            >  1:1229048/1‑77 (MQ=255)
                                                                                                                  gTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACa                                                                                                            >  1:1227615/1‑77 (MQ=255)
                                                                                                                  gTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACa                                                                                                            <  2:1229048/77‑1 (MQ=255)
                                                                                                                             cATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGAc                           <  2:582038/147‑1 (MQ=255)
                                                                                                                             cATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGAc                           >  1:582038/1‑147 (MQ=255)
                                                                                                                             cATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGaa                        <  1:1570161/150‑1 (MQ=255)
                                                                                                                              aTGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACNAat                       <  2:1190688/150‑1 (MQ=255)
                                                                                                                              aTGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAat                       <  2:513400/150‑1 (MQ=255)
                                                                                                                                 gTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATg                    >  2:76271/1‑150 (MQ=255)
                                                                                                                                       gTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCtt              >  1:643131/1‑150 (MQ=255)
                                                                                                                                         gctGGAACATCTCGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGa                            >  2:563217/3‑134 (MQ=255)
                                                                                                                                         gctGGAACATCTCGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGa                            <  1:563217/132‑1 (MQ=255)
                                                                                                                                         aGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTgg            <  1:1292535/150‑1 (MQ=255)
                                                                                                                                           tGGAACATCTCGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGcc          >  1:617510/1‑150 (MQ=255)
                                                                                                                                            ggAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTg                                                          <  2:1819113/101‑1 (MQ=255)
                                                                                                                                            ggAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTg                                                          >  1:1819113/1‑101 (MQ=255)
                                                                                                                                             gAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTtc                                                                 >  2:1781570/1‑93 (MQ=255)
                                                                                                                                             gAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTtc                                                                 <  1:1781570/93‑1 (MQ=255)
                                                                                                                                                   cGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGtt  <  2:1830907/150‑1 (MQ=255)
                                                                                                                                                   cGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGtt  <  2:136684/150‑1 (MQ=255)
                                                                                                                                                    gagaTGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTg                                                          <  1:1405256/93‑1 (MQ=255)
                                                                                                                                                    gagaTGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTg                                                          <  1:1405619/93‑1 (MQ=255)
                                                                                                                                                    gagaTGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTg                                                          >  2:1405256/1‑93 (MQ=255)
                                                                                                                                                    gagaTGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTg                                                          >  2:1405619/1‑93 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGTT  >  NC_002947/5874727‑5875023

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: