Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,410,230 (A)6→7 intergenic (+279/‑50) prmA → / → dusB methyltransferase for 50S ribosomal subunit protein L11/tRNA‑dihydrouridine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,410,2241.A97.8% 294.1 / ‑1.5 91intergenic (+273/‑56)prmA/dusBmethyltransferase for 50S ribosomal subunit protein L11/tRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (48/41);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

ATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGG  >  NC_000913/3410086‑3410357
                                                                                                                                           |                                                                                                                                     
aTTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGc                                                                                                                                        >  1:568388/1‑139 (MQ=255)
      ttGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCaaaaaa                                                                                                                                 <  2:568388/139‑1 (MQ=255)
           tgtAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCaaaaaaatgcg                                                                                                                            <  1:1229090/139‑6 (MQ=255)
                     aGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATAcgc                                                                                                                   >  1:360178/1‑139 (MQ=255)
                                  tAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTca                                                                                                      <  2:360178/139‑1 (MQ=255)
                                     tatGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTg                                                                                                            >  1:1287114/1‑130 (MQ=255)
                                     tatGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTg                                                                                                            <  2:1287114/130‑1 (MQ=255)
                                     tatGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAg                                                                                                   >  1:393886/1‑139 (MQ=255)
                                     tatGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAg                                                                                                   >  2:201003/1‑139 (MQ=255)
                                     tatGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGt                                                                                                        >  1:462438/1‑133 (MQ=255)
                                     tatGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGt                                                                                                        <  2:462438/133‑1 (MQ=255)
                                         atgatTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATg                                                                                               >  2:685519/1‑139 (MQ=255)
                                         atgatTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATg                                                                                               <  1:685519/139‑1 (MQ=255)
                                         atgatTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCaaaaaa                                                                                                                                 <  2:758526/104‑1 (MQ=255)
                                         atgatTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCaaaaaa                                                                                                                                 >  1:758526/1‑104 (MQ=255)
                                                  aTGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCt                                                                                      >  2:1735738/1‑139 (MQ=255)
                                                               gCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCa                                                                         <  2:1678033/139‑1 (MQ=255)
                                                                caccacGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCAt                                                                        <  1:201003/139‑1 (MQ=255)
                                                                 accacGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATc                                                                       >  2:447129/1‑139 (MQ=255)
                                                                   cacGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCgg                                                                     >  2:735867/1‑139 (MQ=255)
                                                                    acGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCaat                                                                                                        <  1:1127969/103‑3 (MQ=255)
                                                                    acGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGt                                                                                                        >  2:1127969/1‑103 (MQ=255)
                                                                     cGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTg                                                                                                            >  1:1358736/1‑98 (MQ=255)
                                                                     cGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTg                                                                                                            <  2:1358736/98‑1 (MQ=255)
                                                                     cGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGAc                                                                   >  2:1057255/1‑139 (MQ=255)
                                                                         cATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAata                                                               <  2:393886/139‑1 (MQ=255)
                                                                         cATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAata                                                               >  1:512405/1‑139 (MQ=255)
                                                                           ttttAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTg                                                                                                            >  1:278675/1‑92 (MQ=255)
                                                                           ttttAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTg                                                                                                            <  2:278675/92‑1 (MQ=255)
                                                                                    cAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATg                                                                             <  2:647567/114‑1 (MQ=255)
                                                                                    cAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATg                                                                             >  1:647567/1‑114 (MQ=255)
                                                                                      aTTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGt                                                                                             >  2:1574/1‑96 (MQ=255)
                                                                                      aTTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGt                                                                                             <  1:1574/96‑1 (MQ=255)
                                                                                           ttGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCt                                             <  2:512405/139‑1 (MQ=255)
                                                                                                 cATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTg                                                                                                            <  1:1846612/70‑1 (MQ=255)
                                                                                                 cATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTg                                                                                                            >  2:1846612/1‑70 (MQ=255)
                                                                                                 cATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGt                                                                                             <  2:1057846/85‑1 (MQ=255)
                                                                                                 cATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGt                                                                                             >  1:1057846/1‑85 (MQ=255)
                                                                                                  aTCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCa                                      <  1:1735738/139‑1 (MQ=255)
                                                                                                  aTCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCa                                      >  1:1606665/1‑139 (MQ=255)
                                                                                                   tctcAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTg                                                                                                            <  1:663768/68‑1 (MQ=255)
                                                                                                   tctcAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTg                                                                                                            >  2:663768/1‑68 (MQ=255)
                                                                                                   tctcAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGaa                                                    <  2:252881/124‑1 (MQ=255)
                                                                                                   tctcAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGaa                                                    >  1:252881/1‑124 (MQ=255)
                                                                                                   tctcAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAg                                     >  2:656027/1‑139 (MQ=255)
                                                                                                   tctcAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATA‑‑‑CGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAgcgc                                  >  1:1818756/1‑139 (MQ=255)
                                                                                                    ctcAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGt                                                                                             >  2:955988/1‑82 (MQ=255)
                                                                                                    ctcAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGt                                                                                             <  1:955988/82‑1 (MQ=255)
                                                                                                    ctcAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCaaagaa                                                                                                                                 >  2:1226772/1‑45 (MQ=255)
                                                                                                    ctcAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCaaaaaa                                                                                                                                 <  1:1226772/45‑1 (MQ=255)
                                                                                                         agGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCa                               >  1:1366047/1‑139 (MQ=255)
                                                                                                            aTTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGAc                                                                   <  1:1309427/100‑1 (MQ=255)
                                                                                                            aTTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGAc                                                                   >  2:1309427/1‑100 (MQ=255)
                                                                                                             ttGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCgg                                                                     <  2:33623/97‑1 (MQ=255)
                                                                                                             ttGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCgg                                                                     >  1:33623/1‑97 (MQ=255)
                                                                                                             ttGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGc                                               >  1:1757546/1‑119 (MQ=255)
                                                                                                             ttGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGc                                               <  2:1757546/119‑1 (MQ=255)
                                                                                                             ttGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCAtggc                           >  1:1551774/1‑139 (MQ=255)
                                                                                                              tGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATAcgccgc                                                                                                                >  2:834661/1‑53 (MQ=39)
                                                                                                              tGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATAcgccgc                                                                                                                <  1:834661/53‑1 (MQ=39)
                                                                                                              tGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTg                                                                                                            <  2:89647/57‑1 (MQ=39)
                                                                                                              tGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTg                                                                                                            >  1:89647/1‑57 (MQ=39)
                                                                                                              tGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTc                                                                                                       >  2:2375/1‑62 (MQ=255)
                                                                                                              tGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTc                                                                                                       <  1:2375/62‑1 (MQ=255)
                                                                                                              tGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCAtggct                          >  2:1634408/1‑139 (MQ=255)
                                                                                                              tGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCAtggct                          >  2:546942/1‑139 (MQ=255)
                                                                                                               ggTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCAtggctg                         <  2:1606665/139‑1 (MQ=255)
                                                                                                               ggTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCATTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCAtggctg                         <  1:735867/139‑1 (MQ=255)
                                                                                                                gTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCAt                                                                                          >  2:813803/1‑73 (MQ=255)
                                                                                                                gTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCAt                                                                                          <  1:813803/73‑1 (MQ=255)
                                                                                                                gTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCgg                                                                     <  2:485617/94‑1 (MQ=255)
                                                                                                                gTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCgg                                                                     >  1:485617/1‑94 (MQ=255)
                                                                                                                gTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCAtggctgg                        >  1:947906/1‑139 (MQ=255)
                                                                                                                gTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCAtggctgg                        >  1:517275/1‑139 (MQ=255)
                                                                                                                   aaaGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGccc                                >  2:1881576/1‑128 (MQ=255)
                                                                                                                   aaaGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGccc                                <  1:1881576/128‑1 (MQ=255)
                                                                                                                     aGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTa                                                                                                                           <  1:1459816/35‑1 (MQ=25)
                                                                                                                     aGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTa                                                                                                                           >  2:1459816/1‑35 (MQ=25)
                                                                                                                     aGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTaa                                                                                   <  2:643541/75‑1 (MQ=255)
                                                                                                                     aGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTaa                                                                                   >  1:643541/1‑75 (MQ=255)
                                                                                                                     aGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATgcg                                                                           >  1:1073918/1‑83 (MQ=255)
                                                                                                                     aGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATgcg                                                                           <  2:1073918/83‑1 (MQ=255)
                                                                                                                     aGTTTGGCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTa                   >  2:396656/1‑139 (MQ=255)
                                                                                                                          ggCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCAt                                                                              >  1:655085/1‑75 (MQ=255)
                                                                                                                          ggCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCAt                                                                              <  2:655085/75‑1 (MQ=255)
                                                                                                                          ggCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAatat                                                              <  2:599196/91‑1 (MQ=255)
                                                                                                                          ggCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAatat                                                              >  1:599196/1‑91 (MQ=255)
                                                                                                                          ggCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAg                                     <  1:1432666/116‑1 (MQ=255)
                                                                                                                          ggCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAg                                     >  2:1432666/1‑116 (MQ=255)
                                                                                                                          ggCCTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTacagac              >  1:171560/1‑139 (MQ=255)
                                                                                                                             cTTTCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTAcagacaga           <  2:517275/139‑1 (MQ=255)
                                                                                                                              tttCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACagacagac          <  1:656027/139‑1 (MQ=255)
                                                                                                                                tcgtgtcGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATACGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCAtggctg                         >  1:1091769/6‑122 (MQ=255)
                                                                                                                                tCATCTCGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCAtggctg                         <  2:1091769/122‑1 (MQ=255)
                                                                                                                                   tctcGTGCAAAAAAATGCGTAATATA‑‑‑CGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCgg  <  2:1818756/139‑1 (MQ=255)
                                                                                                                                      cGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGt                                                                                                  >  1:358369/1‑43 (MQ=255)
                                                                                                                                      cGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGt                                                                                                  <  2:358369/43‑1 (MQ=255)
                                                                                                                                      cGTGCAAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCgg  <  2:1366047/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                     
ATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGG  >  NC_000913/3410086‑3410357

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: