Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 181,943 C→T 26.2% noncoding (52/120 nt) PP_5SB → 5S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947181,9430CT26.2% 186.1 / 80.8 126noncoding (52/120 nt)PP_5SB5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (45/48);  new base T (8/25);  total (53/73)
Fisher's exact test for biased strand distribution p-value = 2.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 135 total aligned reads displayed.

CAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAG  >  NC_002947/181794‑182088
                                                                                                                                                     |                                                                                                                                                  
cAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCACGCACGCCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagc                                                                                                                                                    <  2:1032361/150‑1 (MQ=255)
 aaaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagca                                                                                                                                                   >  2:2292662/1‑150 (MQ=255)
      cccTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagca                                                                                                                                                   <  2:213308/145‑1 (MQ=255)
      cccTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagca                                                                                                                                                   >  1:213308/1‑145 (MQ=255)
                  aTGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTATAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGc                                                                                                                                  >  1:1774042/1‑150 (MQ=255)
                     cAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATGTCTTGACGACCATAGTGCATTGGAACCACCTGACCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGa                                                                                                                               <  2:1977400/150‑1 (MQ=255)
                        gTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATgg                                                                                                                            >  2:232068/1‑150 (MQ=255)
                                 ttGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTggg                                                                                                                   <  1:232068/150‑1 (MQ=255)
                                 ttGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTggg                                                                                                                   <  1:1294763/150‑1 (MQ=255)
                                    tCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTggg                                                                                                                   <  2:2055674/147‑1 (MQ=255)
                                        tAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTccc                                                                                                            <  1:395850/150‑1 (MQ=255)
                                              ggAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCAtgtg                                                                                                      >  2:1763175/1‑150 (MQ=255)
                                               gAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAACGTCTGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTga                                                                                                     >  1:1946287/1‑150 (MQ=21)
                                               gAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTga                                                                                                     <  1:1715504/150‑1 (MQ=255)
                                                    aCAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTa                                                                                                >  1:1096024/1‑150 (MQ=255)
                                                     cAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAg                                                                                               >  1:1978167/1‑150 (MQ=255)
                                                     cAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAg                                                                                               <  2:2361751/150‑1 (MQ=255)
                                                           ttCGCAGGCACGTCCTAGGGCCGGGCCAGCTCTCCGTATTTCTTGACGACCCTAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATc                                                                                         <  2:1946287/150‑1 (MQ=17)
                                                           ttCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGt                                                                                             >  2:361583/1‑146 (MQ=255)
                                                            tCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGACGGTAGTGTGGGGTTTCCCCATGTGAGACTAGTTCATCg                                                                                        >  1:703834/1‑150 (MQ=255)
                                                              gCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTc                                                                                      >  1:1374839/1‑150 (MQ=255)
                                                               cAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCa                                                                                     <  1:615245/150‑1 (MQ=255)
                                                                aGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCaa                                                                                    >  1:2356241/1‑150 (MQ=255)
                                                                           aGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCaa                                                                                    >  1:140687/1‑139 (MQ=255)
                                                                           aGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCaa                                                                                    <  2:140687/139‑1 (MQ=255)
                                                                                     cAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagca                                                                                                                                                   <  2:1504450/66‑1 (MQ=255)
                                                                                     cAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagca                                                                                                                                                   >  1:1504450/1‑66 (MQ=255)
                                                                                         caacaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTga                                                                                                     <  1:1149101/106‑1 (MQ=11)
                                                                                         acacaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTg                                                           >  1:1957614/1‑150 (MQ=255)
                                                                                          aacaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGc                                                          >  2:432885/2‑150 (MQ=255)
                                                                                              gAATTTCTTGACGACCATAGTGCATTGGAACCACCTGCTCGCAGCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATg                                                     <  1:21880/150‑1 (MQ=255)
                                                                                                aTTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCa                                                   <  1:432885/150‑1 (MQ=255)
                                                                                                    cTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACa                                               >  2:2116619/1‑150 (MQ=255)
                                                                                                     ttGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagat                                               <  2:2451923‑M1/150‑32 (MQ=255)
                                                                                                     ttGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagat                                               <  2:1076331‑M1/150‑32 (MQ=255)
                                                                                                     ttGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAg                                              >  1:2045359/1‑150 (MQ=255)
                                                                                                     ttGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAg                                              <  2:529601/150‑1 (MQ=255)
                                                                                                       gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatag                                             <  1:1267372‑M1/150‑34 (MQ=255)
                                                                                                       gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcacaaaaccctacctgcgtgcgcagatagg                                             <  2:1238363‑M1/150‑34 (MQ=255)
                                                                                                       gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAggg                                            >  2:649652/1‑150 (MQ=255)
                                                                                                       gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGGTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAggg                                            <  2:1374839/150‑1 (MQ=255)
                                                                                                       gacgacCATAGAGCATTGGAACCAACTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGCGTAGATCATCATCAAGATTCGtttcgcaaagccgctacctgcatgtgcagatag                                             <  2:1076302‑M1/150‑34 (MQ=255)
                                                                                                        acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatagg                                            <  2:1738139‑M1/150‑35 (MQ=255)
                                                                                                        acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcacaaaacccctacctgcatgtgcagatagg                                            <  2:167515‑M1/150‑35 (MQ=255)
                                                                                                        acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTAAGTAGTAAAACGATGCCTCGCCGATGGTAATATCGCTTGTCCCCATGTGATAGTAGTTCATCGTCAAGATTCAtttctcaaaacccctacctgcatgcgcagatagg                                            <  2:1464343‑M1/150‑35 (MQ=255)
                                                                                                        acgacCATAGAGCATTGGAACCACCTCATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagatagg                                            <  2:1565184‑M1/150‑35 (MQ=255)
                                                                                                         cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                           <  1:2245021‑M1/150‑36 (MQ=255)
                                                                                                         cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                           <  1:26723‑M1/150‑36 (MQ=255)
                                                                                                         cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaaccccgacctgcatgcgcagataggg                                           <  1:1890854‑M1/150‑36 (MQ=255)
                                                                                                          gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagatagggg                                          >  1:2115544‑M1/1‑114 (MQ=255)
                                                                                                          gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatagggg                                          <  1:1016099‑M1/150‑37 (MQ=255)
                                                                                                          gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGa                                                        >  2:1231800/1‑136 (MQ=255)
                                                                                                          gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCt                                         <  2:1825758/150‑1 (MQ=255)
                                                                                                             cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagatagg                                            >  1:269893‑M1/1‑111 (MQ=255)
                                                                                                             cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagatagg                                            <  2:269893‑M1/145‑35 (MQ=255)
                                                                                                             cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAggg                                            <  1:748301/144‑1 (MQ=255)
                                                                                                             cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAggg                                            >  2:748301/1‑144 (MQ=255)
                                                                                                              aTAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTgg                                     <  1:2219569/150‑2 (MQ=255)
                                                                                                               tAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTc                                    >  2:1336184/1‑150 (MQ=255)
                                                                                                                 gagCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTgggtt                                  <  1:104165/150‑5 (MQ=17)
                                                                                                                 gagCATTGGAACCACCTGATCCCATCCCGAACTAAGCAGTGAAATGATGATTCGCCGATGGTAGTGTGGGTTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCtt                                  >  1:1003227/1‑150 (MQ=18)
                                                                                                                     aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatagg                                            >  2:916902‑M1/1‑103 (MQ=255)
                                                                                                                     aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatagg                                            <  1:916902‑M1/137‑35 (MQ=255)
                                                                                                                     aTTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGACATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAg                                              >  1:1334810/1‑134 (MQ=255)
                                                                                                                      ttGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAgg                                             >  2:1643757/1‑134 (MQ=255)
                                                                                                                      ttGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAgg                                             <  1:1643757/134‑1 (MQ=255)
                                                                                                                       tGGAACCGCCTCCTCCCAGACCGAAATCAGCAGCTAAACGATGCATCGACGAGGGTAGCGTCCGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCgg                            <  2:695547/150‑1 (MQ=21)
                                                                                                                        ggAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAg                                              >  1:94807/1‑131 (MQ=255)
                                                                                                                         gAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGcc                                                                   <  1:1732705/110‑1 (MQ=255)
                                                                                                                         gAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGcc                                                                   >  2:1732705/1‑110 (MQ=255)
                                                                                                                         gAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGt                          >  2:1161650/1‑150 (MQ=255)
                                                                                                                         gAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGt                          <  2:2045359/150‑1 (MQ=255)
                                                                                                                         gAACCACCTGATCCCATACCGAACTCAGCAGTGAAACGATGCAt                                                                                                                                     >  1:758152/1‑44 (MQ=255)
                                                                                                                          aaCCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGtt                         <  2:1978167/150‑1 (MQ=255)
                                                                                                                             caccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGc                                                                    >  2:277657/1‑105 (MQ=255)
                                                                                                                             caccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGc                                                                    <  1:277657/105‑1 (MQ=255)
                                                                                                                              accTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCAt                                                                                                         <  1:618720/67‑1 (MQ=37)
                                                                                                                               ccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTg                                                           >  1:1687567/1‑112 (MQ=255)
                                                                                                                               ccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTg                                                           <  2:1687567/112‑1 (MQ=255)
                                                                                                                                cTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAgg                                             >  1:34727/1‑124 (MQ=255)
                                                                                                                                      ccatgccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATc                 <  2:1319160/146‑1 (MQ=255)
                                                                                                                                      ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTGGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTATGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGtt             >  1:2357605/1‑150 (MQ=255)
                                                                                                                                      ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATc                 >  1:1319160/1‑146 (MQ=255)
                                                                                                                                      ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGtt             <  2:1003227/150‑1 (MQ=255)
                                                                                                                                      ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGtt             >  2:1826487/1‑150 (MQ=255)
                                                                                                                                        atcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAg           >  1:2027293/1‑150 (MQ=255)
                                                                                                                                         tcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                           >  2:2213094‑M1/1‑83 (MQ=255)
                                                                                                                                         tcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                           <  1:2213094‑M1/118‑36 (MQ=255)
                                                                                                                                          cccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGaa         <  1:1161650/150‑1 (MQ=255)
                                                                                                                                           ccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagataggg                                           <  2:896231‑M1/116‑36 (MQ=255)
                                                                                                                                           ccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCa                                                   <  2:1320962/107‑1 (MQ=255)
                                                                                                                                           ccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCtt                                        >  2:341857/1‑118 (MQ=255)
                                                                                                                                           ccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAAc        >  2:207089/1‑150 (MQ=255)
                                                                                                                                              aaCTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGa     <  1:2049249/150‑1 (MQ=255)
                                                                                                                                               aCTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGAt    <  1:1445137/150‑1 (MQ=255)
                                                                                                                                                 tCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggggtttt                                      >  2:359363‑M1/1‑75 (MQ=255)
                                                                                                                                                 tCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggggtttt                                      <  1:359363‑M1/115‑41 (MQ=255)
                                                                                                                                                 tCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                           >  1:490909‑M1/1‑75 (MQ=255)
                                                                                                                                                 tCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                           <  2:490909‑M1/110‑36 (MQ=255)
                                                                                                                                                    gcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTgaga                                                                                                   <  1:458168/51‑1 (MQ=33)
                                                                                                                                                     |                                                                                                                                                  
CAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAG  >  NC_002947/181794‑182088

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