Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,327,127 C→T 56.6% intergenic (+97/‑196) PP_16SE → / → PP_23SE 16S ribosomal RNA/23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,327,1270CT56.6% ‑5.4 / 87.0 76intergenic (+97/‑196)PP_16SE/PP_23SE16S ribosomal RNA/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (17/15);  new base T (22/20);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

AGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAA  >  NC_002947/1326981‑1327235
                                                                                                                                                  |                                                                                                            
aGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGACTACTCTTTGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATc                                                                                                           <  1:968806/150‑1 (MQ=11)
     tAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATg                                                                                                      >  1:112289/1‑150 (MQ=17)
     tAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATg                                                                                                      >  2:29994/1‑150 (MQ=17)
           ggTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGAtt                                                                                                <  2:112289/150‑1 (MQ=17)
            gTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGAttt                                                                                               >  1:416195/1‑150 (MQ=32)
                  gTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCt                                                                                         >  1:265628/1‑150 (MQ=32)
                   tAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTg                                                                                            >  2:20822/1‑146 (MQ=17)
                   tAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTg                                                                                            <  1:20822/146‑1 (MQ=17)
                      gggAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTg                                                                                     >  1:1241801/1‑150 (MQ=32)
                        gAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTc                                                                                   >  2:1335305/1‑150 (MQ=17)
                        gAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTc                                                                                   <  2:416195/150‑1 (MQ=32)
                           ccTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGCCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCACTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGa                                                                                >  2:824477/1‑150 (MQ=17)
                            cTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGAt                                                                               >  2:264971/1‑150 (MQ=17)
                                ggCTGGATCACCTCCTTAATCGACGACATCCGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGtt                                                                           <  2:985139/150‑1 (MQ=17)
                                 gCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTc                                                                          <  2:678480/150‑1 (MQ=17)
                                 gCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTc                                                                          >  1:1201959/1‑150 (MQ=17)
                                     gATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGGCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCCTGGTATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTa                                                                      <  1:824477/150‑1 (MQ=255)
                                           ctcctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAGAAAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAAtt                                                                >  1:728404/1‑150 (MQ=255)
                                           ctcctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAAtt                                                                >  1:466682/1‑150 (MQ=17)
                                            tcctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACCCGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTc                                                               <  1:944549/150‑1 (MQ=17)
                                            tcctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTc                                                               <  2:265628/150‑1 (MQ=32)
                                            tcctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAGAAGACGTTTTTCCTGGTCAGAACCTGGCAATGAGCATACGACACGATGGATGTGTACTGGCTTTTGTGAGATCGTTCTTTAAAAAttt                                                               >  2:944549/1‑149 (MQ=11)
                                             cctTAATCGACGCCATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCg                                                              <  1:648820/150‑1 (MQ=17)
                                             cctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCg                                                              <  2:950024/150‑1 (MQ=17)
                                              ctTAATCGACGACATCAGCCTGCTTAGGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCgg                                                             >  2:1194917/1‑150 (MQ=12)
                                              ctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTGAGATCGTTCTTTAAAAATTCgg                                                             >  2:937002/1‑150 (MQ=17)
                                              ctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCCGATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCgg                                                             >  1:312541/1‑150 (MQ=17)
                                               ttAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGa                                                            <  2:466682/150‑1 (MQ=17)
                                                 aaTCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGata                                                          <  1:607297/150‑1 (MQ=32)
                                                  aTCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGatat                                                         >  1:1381499/1‑150 (MQ=32)
                                                     gacgacATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAtgt                                                       >  1:42225/1‑149 (MQ=17)
                                                     gacgacATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAtgt                                                       <  2:42225/149‑1 (MQ=17)
                                                      gcgacATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAg                                                                                 <  2:1161605/121‑1 (MQ=18)
                                                      gcgacATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAg                                                                                 >  1:1161605/2‑122 (MQ=17)
                                                      acgacATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTga                                                     >  1:1256924/1‑150 (MQ=17)
                                                      acgacATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTga                                                     <  1:264971/150‑1 (MQ=17)
                                                      acgacATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTga                                                     <  2:1201959/150‑1 (MQ=17)
                                                      acgacATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTga                                                     >  2:688574/1‑150 (MQ=32)
                                                        gacATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTg                                                                                            >  2:1007423/1‑109 (MQ=18)
                                                        gacATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTg                                                                                            <  1:1007423/109‑1 (MQ=17)
                                                           aTCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgatagat                                                >  1:517264/1‑150 (MQ=17)
                                                           aTCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGa                                                            >  2:1167716/1‑138 (MQ=33)
                                                           aTCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGa                                                            <  1:1167716/138‑1 (MQ=33)
                                                           aTCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgatagat                                                >  1:292124/1‑150 (MQ=32)
                                                           aTCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgatagat                                                >  2:1401831/1‑150 (MQ=32)
                                                            tCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgatagata                                               <  1:999161/150‑1 (MQ=17)
                                                             cAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTg                                                                                            <  1:561498/104‑1 (MQ=33)
                                                             cAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTg                                                                                            >  2:561498/1‑104 (MQ=33)
                                                             cAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGatat                                              <  1:1401831/150‑1 (MQ=32)
                                                                  tgctgATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACt                                         >  2:775012/1‑150 (MQ=17)
                                                                   gctgATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTg                                        <  2:292124/150‑1 (MQ=32)
                                                                      gatgaGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAac                                     >  2:1409940/1‑150 (MQ=32)
                                                                         gaGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACAcc                                  >  2:96321/1‑150 (MQ=32)
                                                                          aGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCa                                 <  1:171225/150‑1 (MQ=17)
                                                                           gCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAg                                <  2:1241801/150‑1 (MQ=32)
                                                                               ccACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGatat                                              <  1:1087317/132‑1 (MQ=18)
                                                                               ccACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGatat                                              >  2:1087317/1‑132 (MQ=17)
                                                                               ccACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTc                            >  1:1256728/1‑150 (MQ=17)
                                                                               ccACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTc                            >  1:1351411/1‑150 (MQ=17)
                                                                                cacaCGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCa                           <  2:1256924/150‑1 (MQ=17)
                                                                                 acacTAATTTCTTGATTCATTCTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGATCATTCTCATCGAATGTTGATTTCTGGCTTTTGTCATATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCAc                          <  1:542305/150‑1 (MQ=11)
                                                                                   acGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCActg                        >  1:1388964/1‑150 (MQ=17)
                                                                                      aaTTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCActgctg                     <  2:1003179/150‑1 (MQ=17)
                                                                                           cTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCtttt                                                                                      >  2:419613/1‑80 (MQ=255)
                                                                                           cTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCtttt                                                                                      <  1:419613/80‑1 (MQ=255)
                                                                                           cTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGa                >  1:246758/1‑150 (MQ=32)
                                                                                            ttGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGAt               >  1:860027/1‑150 (MQ=32)
                                                                                              gATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCa             <  1:1289680/150‑1 (MQ=32)
                                                                                               aTTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAg            <  1:1335305/150‑1 (MQ=17)
                                                                                               aTTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAg            <  2:1381499/150‑1 (MQ=32)
                                                                                                 tCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGc          <  1:1409940/150‑1 (MQ=32)
                                                                                                  cATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCt         >  1:1132126/1‑150 (MQ=32)
                                                                                                    ttGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTaa       <  2:576568/150‑1 (MQ=17)
                                                                                                       tCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGt    >  1:496455/1‑150 (MQ=17)
                                                                                                       tCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGt    >  1:1134003/1‑150 (MQ=17)
                                                                                                         gAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTaa  >  1:522880/1‑150 (MQ=11)
                                                                                                                                                  |                                                                                                            
AGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAA  >  NC_002947/1326981‑1327235

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: