Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 181,943 C→T 23.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,9430CT23.2% 239.3 / 74.1 138noncoding (52/120 nt)PP_5SB5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (49/56);  new base T (21/11);  total (70/68)
Fisher's exact test for biased strand distribution p-value = 7.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 143 total aligned reads displayed.

AAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTG  >  NC_002947/181795‑182091
                                                                                                                                                     |                                                                                                                                                     
aaaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagca                                                                                                                                                      >  2:1127958/1‑150 (MQ=255)
  aaCCCCTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGt                                                                                                                                                    >  2:364025/1‑150 (MQ=255)
    ccccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGa                                                                                                                                                  >  2:1397264/1‑150 (MQ=37)
     cccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaa                                                                                                                                                 >  1:1847464/1‑150 (MQ=37)
       cTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAAc                                                                                                                                               >  1:235727/1‑150 (MQ=255)
       cTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGAAGGGCATCAGGATAGGAGTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATATCCCACTCAGCAGTTAAAc                                                                                                                                               >  1:2491953/1‑150 (MQ=25)
         aTCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGa                                                                                                                                             >  1:1237406/1‑150 (MQ=255)
                           gtgTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTa                                                                                                                             <  1:2124477/146‑1 (MQ=255)
                           gtgTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTa                                                                                                                             >  2:2124477/3‑148 (MQ=255)
                                ttGTCTTTAAGTAGGAACTACAGAGATTGGCAGACACGTCCGATAGACGGGCCATTAGACAGAATTTCTTGAAGAAGAACTACCCGAGCATGCACGCGATCCCATCCCGAACTCAGCAGTGAAACGATGCATAGACGATTGGAGTGTggg                                                                                                                      >  1:1466527/1‑150 (MQ=255)
                                    cTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGa                                                                                                                                             >  2:1761878/1‑123 (MQ=255)
                                    cTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGa                                                                                                                                             <  1:1761878/123‑1 (MQ=255)
                                     tttAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTc                                                                                                                 >  1:1276141/1‑150 (MQ=255)
                                     tttAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTc                                                                                                                 <  2:1405559/150‑1 (MQ=255)
                                     tttAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTc                                                                                                                 >  2:2770229/1‑150 (MQ=255)
                                       tAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTccc                                                                                                               <  1:1700849/150‑1 (MQ=255)
                                              gAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTga                                                                                                        >  2:2708368/1‑150 (MQ=255)
                                              gAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTga                                                                                                        >  2:2708127/1‑150 (MQ=255)
                                              gAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTga                                                                                                        >  2:3010034/1‑150 (MQ=255)
                                                     agagATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTa                                                                                                                             >  2:666343/1‑122 (MQ=255)
                                                     agagATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAgg                                                                                                 >  2:1195553/1‑150 (MQ=255)
                                                      gagaTTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTc                                                                                                                 <  2:937695/133‑1 (MQ=255)
                                                      gagaTTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTGTc                                                                                                                 >  1:937695/1‑133 (MQ=255)
                                                        gaTTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCa                                                                                              <  1:364025/150‑1 (MQ=255)
                                                          ttCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATc                                                                                            <  1:2770229/150‑1 (MQ=255)
                                                          ttCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACCATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATc                                                                                            <  1:60312/150‑1 (MQ=255)
                                                             gCAGGCACGTCCTAGGGACGGGCAAGCAAACAGAATTTCTTGACGACCATAGAGCATTGGAACCACATGATCCCATCCCGACCTCAGCAGTGAAACGATGCATCGCCGATGGTAATGTGGGGTTTCCCCATGTGAGATAAGGTCATCGTc                                                                                         >  2:1057545/1‑150 (MQ=17)
                                                               aGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCaa                                                                                       <  1:851008/150‑1 (MQ=255)
                                                                  cACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGAt                                                                                    <  2:2212372/150‑1 (MQ=255)
                                                                  cACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGAGCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCCAAGGTAGTGTGGAGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGAt                                                                                    <  1:1057545/150‑1 (MQ=255)
                                                                   aCGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGAtt                                                                                   <  1:1195553/150‑1 (MQ=255)
                                                                      tCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCaa                                                                                >  2:1971132/1‑150 (MQ=255)
                                                                       ccTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCaaa                                                                               <  2:235727/150‑1 (MQ=255)
                                                                        cctaggACGTGCCTGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAAc                                                                              <  1:697181/146‑1 (MQ=255)
                                                                         tAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCa                                                                                                                                         >  1:1409781/1‑90 (MQ=255)
                                                                         aaGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGTGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAAcc                                                                             >  2:2099875/2‑150 (MQ=255)
                                                                              aCGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAg                                                                        >  2:1840783/1‑150 (MQ=255)
                                                                              aCGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAg                                                                        <  1:1971132/150‑1 (MQ=255)
                                                                                gggCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAg                                                                                                  >  2:214264/1‑122 (MQ=255)
                                                                                gggCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAg                                                                                                  <  1:214264/122‑1 (MQ=255)
                                                                                        acacaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTggg                                                                                                                      >  1:2359607/1‑94 (MQ=37)
                                                                                        acacaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTggg                                                                                                                      <  2:2359607/94‑1 (MQ=37)
                                                                                          cacaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGc                                                             >  2:198295/1‑150 (MQ=255)
                                                                                               aaTTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGc                                                       <  2:464129/150‑1 (MQ=255)
                                                                                                 tttCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTcc                                                                                                                >  1:174612/1‑92 (MQ=34)
                                                                                                 tttCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGcatc                                                          >  1:1004500/1‑143 (MQ=255)
                                                                                                   tCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcag                                                    <  2:1757258‑M2/150‑30 (MQ=255)
                                                                                                    cTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcaga                                                   >  2:941591‑M2/1‑120 (MQ=255)
                                                                                                       gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctatctgcatgcgcagatag                                                >  2:1850588‑M2/1‑117 (MQ=255)
                                                                                                       gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagatag                                                >  2:2961504‑M2/1‑117 (MQ=255)
                                                                                                       gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAggg                                               <  2:1847464/150‑1 (MQ=255)
                                                                                                        acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatagg                                               >  1:476172‑M2/1‑116 (MQ=255)
                                                                                                        acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt                                                                                                                                            <  1:1807437/57‑1 (MQ=38)
                                                                                                        acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAgggg                                              >  1:743823/1‑150 (MQ=255)
                                                                                                        acgacCAAAGGGAATTGAAACCACCTGATCCCATCCCGAACTCAGTATTGAAACGATGCATCGCCGTTGGTAGTGTGGGGTTCCCCCATGTGAGATTGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                              >  2:146024‑M2/1‑115 (MQ=255)
                                                                                                         cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagataggg                                              >  1:1156050‑M2/1‑115 (MQ=255)
                                                                                                         cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                              >  2:2016738‑M2/1‑115 (MQ=255)
                                                                                                          gacCATAGAGAATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACACAAAGCCCCTGTCTGCATCGCAGACAggggtc                                            >  2:145145/1‑148 (MQ=17)
                                                                                                            ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcacaaaacccctacctgcatgcgcagataggggtt                                           <  2:585917‑M2/150‑39 (MQ=255)
                                                                                                             cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggggttt                                          >  2:183779‑M2/1‑111 (MQ=255)
                                                                                                              aTAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAggg                                               <  1:2339648/143‑1 (MQ=255)
                                                                                                               tAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggggttttg                                        >  2:1727320‑M2/1‑109 (MQ=255)
                                                                                                               tAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTggg                                                                                                                      >  2:323766/1‑72 (MQ=35)
                                                                                                               tAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTc                                       >  1:2284270/1‑150 (MQ=255)
                                                                                                                 gagCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTgc                                                                 >  1:332225/1‑122 (MQ=255)
                                                                                                                  agCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAACCGATGCATCGCCGATGGTAGTGTGGCGTTCCCCCATGTGACAGTACGTCATCGTCAAGATTCAAACCCCAACCCCCTGTCCGCATCCCAGACAGTGGCTTTgggtgt                                    <  2:2491953/150‑6 (MQ=14)
                                                                                                                   gCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCtttt                                   <  1:198295/150‑1 (MQ=255)
                                                                                                                    cATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGccc                                                                     <  2:1733205/116‑1 (MQ=255)
                                                                                                                    cATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGccc                                                                     >  1:1733205/1‑116 (MQ=255)
                                                                                                                    cATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTc                                  >  2:2167967/1‑150 (MQ=255)
                                                                                                                     aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagataggg                                              <  1:1157635‑M2/138‑36 (MQ=255)
                                                                                                                     aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagataggg                                              >  2:1157635‑M2/1‑103 (MQ=255)
                                                                                                                     aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGa                                                    <  1:1345990/131‑1 (MQ=255)
                                                                                                                     aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGa                                                    >  2:1345990/1‑131 (MQ=255)
                                                                                                                      ttGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGAt                                                                                                                                 >  2:2465627/1‑54 (MQ=38)
                                                                                                                      ttGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGAt                                                                                                                                 <  1:2465627/54‑1 (MQ=38)
                                                                                                                       tGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCgg                               <  1:963037/150‑1 (MQ=255)
                                                                                                                         gAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAgtgt                                                                                                                         <  2:1036809/59‑1 (MQ=35)
                                                                                                                            ccaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcacaaaacccctacctgcatgcgcagatag                                                <  2:1540255‑M2/129‑34 (MQ=255)
                                                                                                                            ccaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcacaaaacccctacctgcatgcgcagatag                                                >  1:1540255‑M2/1‑96 (MQ=255)
                                                                                                                             caccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCt                                                                   <  1:2949057/109‑1 (MQ=255)
                                                                                                                             caccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCt                                                                   >  2:2949057/1‑109 (MQ=255)
                                                                                                                              accTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCaa                                                                                       >  1:2736291/1‑88 (MQ=35)
                                                                                                                              accTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCaa                                                                                       <  2:2736291/88‑1 (MQ=35)
                                                                                                                                   attccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCATTTCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCg                   <  2:1776223/147‑1 (MQ=255)
                                                                                                                                   atcccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCg                   <  2:2284270/150‑1 (MQ=255)
                                                                                                                                     cccatcccGAACTCAGAAGTGAAACGACGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGt                 <  1:1352233/150‑1 (MQ=255)
                                                                                                                                       catcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTga                                                                                                        >  1:1583935/1‑62 (MQ=37)
                                                                                                                                       catcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTga                                                                                                        <  2:1583935/62‑1 (MQ=37)
                                                                                                                                       catcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTa               <  1:2167967/150‑1 (MQ=255)
                                                                                                                                        atcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                              >  2:2294777‑M2/1‑84 (MQ=255)
                                                                                                                                        atcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTc                                                                <  2:2931963/101‑1 (MQ=255)
                                                                                                                                        atcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTc                                                                >  1:2931963/1‑101 (MQ=255)
                                                                                                                                         tcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGa             >  2:2419403/1‑150 (MQ=255)
                                                                                                                                          cccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCGGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGaa            <  1:3010034/150‑1 (MQ=255)
                                                                                                                                           ccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTCCCGGGTTTCAATATCGGTTAGAAc           <  2:1276141/150‑1 (MQ=255)
                                                                                                                                              aaCTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGa        <  1:1841511/150‑1 (MQ=255)
                                                                                                                                              aaCTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGa        <  2:1466527/150‑1 (MQ=255)
                                                                                                                                                   agcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGt                 <  1:1061105/136‑1 (MQ=255)
                                                                                                                                                    gcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTg  >  1:2206264/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTG  >  NC_002947/181795‑182091

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