Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 182,034 Δ1 bp 9.3% intergenic (+23/+42) PP_5SB → / ← PP_0160 5S ribosomal RNA/ferrioxamine receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947182,0340G. 9.3% 310.2 / 30.0 107intergenic (+23/+42)PP_5SB/PP_01605S ribosomal RNA/ferrioxamine receptor
Reads supporting (aligned to +/- strand):  ref base G (50/44);  new base . (6/4);  total (57/50)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

Only 100 of 107 total aligned reads displayed.

AATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTG  >  NC_002947/181889‑182181
                                                                                                                                                 |                                                                                                                                                    
aaTTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGc                                                                                                                                                 <  2:464129/150‑1 (MQ=255)
        gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAggg                                                                                                                                         <  2:1847464/150‑1 (MQ=255)
         acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAgggg                                                                                                                                        >  1:743823/1‑150 (MQ=255)
           gacCATAGAGAATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACACAAAGCCCCTGTCTGCATCGCAGACAggggtc                                                                                                                                      >  2:145145/1‑148 (MQ=17)
              cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAg                                                                                                                                           >  1:2725159/1‑142 (MQ=255)
              cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAg                                                                                                                                           <  2:2725159/142‑1 (MQ=255)
               aTAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAggg                                                                                                                                         <  1:2339648/143‑1 (MQ=255)
               aTAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAggg                                                                                                                                         >  2:2339648/1‑143 (MQ=255)
                tAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTc                                                                                                                                 >  1:2284270/1‑150 (MQ=255)
                   agCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAACCGATGCATCGCCGATGGTAGTGTGGCGTTCCCCCATGTGACAGTACGTCATCGTCAAGATTCAAACCCCAACCCCCTGTCCGCATCCCAGACAGTGGCTTTgggtgt                                                                                                                              <  2:2491953/150‑6 (MQ=14)
                    gCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCtttt                                                                                                                             <  1:198295/150‑1 (MQ=255)
                     cATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTc                                                                                                                            >  2:2167967/1‑150 (MQ=255)
                      aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGa                                                                                                                                              <  1:1345990/131‑1 (MQ=255)
                        tGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCgg                                                                                                                         <  1:963037/150‑1 (MQ=255)
                               accTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCaa                                                                                                                  <  1:1127958/150‑1 (MQ=255)
                               accTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCaa                                                                                                                  >  1:875935/1‑150 (MQ=255)
                                    attccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCATTTCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCg                                                                                                             <  2:1776223/147‑1 (MQ=255)
                                    atcccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCg                                                                                                             <  2:2284270/150‑1 (MQ=255)
                                      cccatcccGAACTCAGAAGTGAAACGACGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGt                                                                                                           <  1:1352233/150‑1 (MQ=255)
                                        catcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTa                                                                                                         <  1:2167967/150‑1 (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)
                                                aCTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGAt                                                                                                 >  1:1175057/1‑150 (MQ=255)
                                                    agcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGt                                                                                                           >  2:1061105/1‑136 (MQ=255)
                                                    agcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGt                                                                                                           <  1:1061105/136‑1 (MQ=255)
                                                     gcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTg                                                                                            >  1:2206264/1‑150 (MQ=255)
                                                       agTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGAc                                                                                          >  1:2748533/1‑150 (MQ=255)
                                                             aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATTTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTctgg                                                                                                                         >  2:2266537/1‑110 (MQ=255)
                                                             aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATTTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCTgg                                                                                                                         <  1:2266537/113‑1 (MQ=255)
                                                             aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAg                                                                                                                                           >  2:2980230/1‑95 (MQ=255)
                                                             aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAg                                                                                                                                           <  1:2980230/95‑1 (MQ=255)
                                                                 tGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTc                                                                                >  1:1114415/1‑150 (MQ=255)
                                                                    aTCGCCGAGGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTcacc                                                                             >  2:2149370/1‑149 (MQ=255)
                                                                        ccGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTgc                                                                         <  1:2419403/150‑1 (MQ=255)
                                                                          gATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTgcgc                                                                       >  2:92835/1‑150 (MQ=255)
                                                                              gTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCtt                                                                   <  2:1175057/150‑1 (MQ=255)
                                                                                  tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCAcc                                                               >  2:1900990/1‑150 (MQ=255)
                                                                                  tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCAcc                                                               >  2:1994204/1‑150 (MQ=255)
                                                                                    tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTgc                                                                         >  1:1514748/1‑138 (MQ=255)
                                                                                    tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTgc                                                                         <  2:1514748/138‑1 (MQ=255)
                                                                                    tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCa                                                             >  1:32507/1‑150 (MQ=255)
                                                                                    tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCa                                                             >  1:812594/1‑150 (MQ=255)
                                                                                    tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCa                                                             >  2:1891722/1‑150 (MQ=255)
                                                                                    tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCGGTCTGCGATGCAGCCAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTTGCGCTCACATTGCGCTCTTCACCCa                                                             >  1:781445/1‑150 (MQ=255)
                                                                                     ggggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCt                                                                                             <  1:3049447/117‑1 (MQ=255)
                                                                                     ggggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTATCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCt                                                                                             >  2:3049447/1‑117 (MQ=255)
                                                                                             cccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCAc                                                                >  1:1497713/1‑138 (MQ=255)
                                                                                             cccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCAc                                                                <  2:1497713/138‑1 (MQ=255)
                                                                                             cccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCagt                                                    >  2:1803151/1‑150 (MQ=255)
                                                                                              ccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAGGGGCttt                                                                                                                                    >  2:887500/1‑69 (MQ=21)
                                                                                              ccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAGGGGCttt                                                                                                                                    <  1:887500/69‑1 (MQ=21)
                                                                                              ccATGTGAGAGTAGGTCATCGTCAAGATTCAAACACAAAGCCCCTGTCTGCATCGCAGACAg                                                                                                                                           >  2:2674365/1‑62 (MQ=18)
                                                                                              ccATGTGAGAGTAGGTCATCGTCAAGATTCAAACACAAAGCCCCTGTCTGCATCGCAGACAg                                                                                                                                           <  1:2674365/62‑1 (MQ=18)
                                                                                                aTGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCagtagt                                                 <  2:2748533/150‑1 (MQ=255)
                                                                                                  gtgAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTca                                                                               <  1:1187096/118‑1 (MQ=255)
                                                                                                  gtgAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTca                                                                               >  2:1187096/1‑118 (MQ=255)
                                                                                                  gtgAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTAACCCATCTGGCCGAAGTTc                                               >  2:660197/1‑150 (MQ=255)
                                                                                                   tgAGACTATTTCATCGTCAAGATTCAAACCCTAACCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCa                                              <  2:1114415/150‑1 (MQ=255)
                                                                                                    gagaGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAg                                             >  1:444307/1‑150 (MQ=255)
                                                                                                      gagTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGCTCAGGc                                           <  2:2206264/150‑1 (MQ=255)
                                                                                                          aGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGc                                       >  1:1532738/1‑150 (MQ=255)
                                                                                                            gTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCa                                              <  1:1127385/141‑1 (MQ=255)
                                                                                                            gTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCa                                              >  2:1127385/1‑141 (MQ=255)
                                                                                                                   tCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCa                                                                 <  1:2752418/115‑1 (MQ=255)
                                                                                                                   tCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGt                              <  2:743823/150‑1 (MQ=255)
                                                                                                                     aggATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTcaca                                                                             <  2:2629847/99‑1 (MQ=255)
                                                                                                                      aGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTTACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGAt                           >  2:550975/1‑150 (MQ=255)
                                                                                                                       gATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGAt                                                                                                 >  2:582177/1‑79 (MQ=255)
                                                                                                                       gATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGAt                                                                                                 <  1:582177/79‑1 (MQ=255)
                                                                                                                        aTTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGt                         >  2:38061/1‑150 (MQ=255)
                                                                                                                         ttCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTa                        <  1:2708127/150‑1 (MQ=255)
                                                                                                                         ttCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTa                        <  1:2708368/150‑1 (MQ=255)
                                                                                                                          tCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAg                       >  1:1592849/1‑150 (MQ=255)
                                                                                                                             aaCCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTCTGTCTTTTCCGGGTTTCAATATCGGTTAGACCTGACAGCTGACCGTGGCGCTCACATCGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACACCAGGGGATGTAGCtt                    >  1:2698417/1‑150 (MQ=255)
                                                                                                                                ccAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGa                                                                                                  <  2:1875948/69‑1 (MQ=255)
                                                                                                                                ccAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGa                                                                                                  >  1:1875948/1‑69 (MQ=255)
                                                                                                                                  aaaGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCAtt               <  1:1891722/150‑1 (MQ=255)
                                                                                                                                  aaaGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCAGAGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCAtt               <  1:660197/150‑1 (MQ=255)
                                                                                                                                  aaaGCCCCTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAgg                                            <  1:1324805/121‑1 (MQ=255)
                                                                                                                                  aaaGCCCCTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAgg                                            >  2:1324805/1‑121 (MQ=255)
                                                                                                                                  aaaGCCCCTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCATGgcgc                                                                                  <  2:3030936/83‑1 (MQ=25)
                                                                                                                                  aaaGCCCCTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCATGgcgc                                                                                  >  1:3030936/1‑83 (MQ=25)
                                                                                                                                   aaGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTg              <  2:1592849/150‑1 (MQ=255)
                                                                                                                                   aaGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTg              <  1:92835/150‑1 (MQ=255)
                                                                                                                                       cccTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTaa          >  1:2447891/1‑150 (MQ=255)
                                                                                                                                        ccTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTGGCTTTCATTGGTAAg         >  1:2353619/1‑150 (MQ=255)
                                                                                                                                        ccTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAg         <  2:875935/150‑1 (MQ=255)
                                                                                                                                        ccTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTCCAGGCTGGCACATGAgg                               <  1:1615938/128‑1 (MQ=255)
                                                                                                                                        ccTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAg         >  2:2465722/1‑150 (MQ=255)
                                                                                                                                        ccTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTCCAGGCTGGCACATGAgg                               >  2:1615938/1‑128 (MQ=255)
                                                                                                                                         cTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAgg        >  2:928959/1‑150 (MQ=255)
                                                                                                                                         cTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAgg        >  1:1099864/1‑150 (MQ=255)
                                                                                                                                            tCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTCTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTgc                                                                         <  2:2168660/82‑1 (MQ=255)
                                                                                                                                            tCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTCTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTgc                                                                         >  1:2168660/1‑82 (MQ=255)
                                                                                                                                              tGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGgttgtt   >  2:795993/1‑150 (MQ=255)
                                                                                                                                               gCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCTCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGttgttg  <  1:2149370/150‑1 (MQ=255)
                                                                                                                                               gCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGttgttg  >  2:148833/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTG  >  NC_002947/181889‑182181

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: