Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,410,240 Δ1 bp intergenic (+289/‑40) 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,2400T.93.3% 240.3 / 14.5 60intergenic (+289/‑40)prmA/dusBmethyltransferase for 50S ribosomal subunit protein L11/tRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base . (34/22);  total (37/23)
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.90e-01

GTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGAT  >  NC_000913/3410096‑3410385
                                                                                                                                                |                                                                                                                                                 
gtgtAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAAcgccgc                                                                                                                                            <  2:8261/151‑1 (MQ=255)
  gtAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCt                                                                                                                                          >  1:1071610/1‑151 (MQ=255)
  gtAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCt                                                                                                                                          <  1:295586/151‑1 (MQ=255)
               aCTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGt                                                                                                                              >  1:691425/1‑150 (MQ=255)
                  aTTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATg                                                                                                                           <  1:985905/150‑1 (MQ=255)
                   ttttCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTAAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGt                                                                                                                         >  2:1010813/1‑151 (MQ=255)
                     ttCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTc                                                                                                                        <  2:1055306/150‑1 (MQ=255)
                      tCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCa                                                                                                                       <  2:860026/150‑1 (MQ=255)
                         aaTATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCAtt                                                                                                                     <  2:768930/149‑1 (MQ=255)
                         aaTATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCAt                                                                                                                      <  2:561646/148‑1 (MQ=255)
                         aaTATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTc                                                                                                                   >  2:720389/1‑151 (MQ=255)
                            atgatgATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTa                                                                                                                <  1:518469/151‑1 (MQ=255)
                              gatgatTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAAc                                                                                                              >  2:174946/1‑151 (MQ=255)
                              gatgatTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAAc                                                                                                              >  2:630622/1‑151 (MQ=255)
                               atgatTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAAc                                                                                                              <  1:85295/150‑1 (MQ=255)
                               atgatTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACt                                                                                                             >  1:486487/1‑151 (MQ=255)
                                  atTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCAt                                                                                                          >  1:500756/1‑151 (MQ=255)
                                      ttATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATgcgc                                                                                                      >  1:602448/1‑151 (MQ=255)
                                      ttATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATgcg                                                                                                       <  2:788504/150‑1 (MQ=255)
                                       tATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGGTTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATgcgc                                                                                                      >  1:7985/1‑151 (MQ=255)
                                            gTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGa                                                                                                <  1:921572/151‑1 (MQ=255)
                                              aaTTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGAc                                                                                               >  2:701164/1‑150 (MQ=255)
                                                  aTCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAatat                                                                                          <  1:606839/151‑1 (MQ=255)
                                                        ccacGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTc                                                                                    >  1:806295/1‑151 (MQ=255)
                                                            gCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGa                                                                                 >  1:595210/1‑150 (MQ=255)
                                                            gCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGa                                                                                 >  1:537055/1‑150 (MQ=255)
                                                                aTTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCg                                                                            <  1:861425/151‑1 (MQ=255)
                                                                 ttttAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGc                                                                           >  1:380027/1‑151 (MQ=255)
                                                                      aaTGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCAGAt                                                                      <  2:378047/151‑1 (MQ=255)
                                                                             ttCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAgcg                                                               >  1:924050/1‑151 (MQ=255)
                                                                             ttCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAgcg                                                               >  1:665302/1‑151 (MQ=255)
                                                                                tttGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGccc                                                            >  2:319892/1‑151 (MQ=255)
                                                                                 ttGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGccc                                                            >  1:1001495/1‑151 (MQ=255)
                                                                                  tGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCa                                                           <  2:944410/150‑1 (MQ=255)
                                                                                       cATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCAtggctg                                                     <  1:464172/151‑1 (MQ=255)
                                                                                         tctcAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCAtggctgg                                                    >  1:39608/1‑151 (MQ=255)
                                                                                            cAGAGGGTTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCAt                                                 <  2:7985/151‑1 (MQ=255)
                                                                                             agagGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTa                                               <  2:438299/151‑1 (MQ=255)
                                                                                                ggATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTacag                                            <  2:1124970/151‑1 (MQ=255)
                                                                                                 gATTAGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTac                                              >  2:113099/1‑148 (MQ=255)
                                                                                                   ttGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTacagac                                          >  1:791470/1‑150 (MQ=255)
                                                                                                          aaGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCAt                                                          <  1:381599/127‑1 (MQ=255)
                                                                                                          aaGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCAt                                                          >  2:381588/1‑127 (MQ=255)
                                                                                                             tttGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCg                               >  1:555502/1‑151 (MQ=255)
                                                                                                             tttGGCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCg                               >  2:206395/1‑151 (MQ=255)
                                                                                                                ggCCTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGAc                            <  1:203372/151‑1 (MQ=255)
                                                                                                                   cTTTCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGtt                         >  2:373864/1‑151 (MQ=255)
                                                                                                                    tttCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTtg                        >  2:758247/1‑151 (MQ=255)
                                                                                                                     ttCATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTtg                        <  2:1098304/150‑1 (MQ=255)
                                                                                                                       cATCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTtgtg                      <  2:271507/150‑1 (MQ=255)
                                                                                                                       cATCTCGCGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCAGGGCTGGCATTACAGACAGACCTTTTCGGACGTtgtg                      >  2:271529/1‑150 (MQ=255)
                                                                                                                        aTCTCGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCt                    >  1:449670/1‑151 (MQ=255)
                                                                                                                           tcGTGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACg                 >  1:211337/1‑151 (MQ=255)
                                                                                                                              tGCAAAAAATGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACg                 >  2:1041919/1‑148 (MQ=255)
                                                                                                                                 aaaaaaTGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATgg           <  2:602434/151‑1 (MQ=255)
                                                                                                                                  aaaaaTGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATggg          >  1:972929/1‑151 (MQ=255)
                                                                                                                                  aaaaaTGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATggg          >  2:407895/1‑151 (MQ=255)
                                                                                                                                     aaTGCGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGc       >  2:93285/1‑151 (MQ=255)
                                                                                                                                         cGTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTAGGACGTTGTGCTACGAGATGGGAGCCg     >  2:939258/1‑149 (MQ=255)
                                                                                                                                          gTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGAt  >  2:149547/1‑151 (MQ=255)
                                                                                                                                          gTAATAACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGAt  >  1:1089149/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCGCAGCGCCCATGGCTGGCATTACAGACAGACCTTTTCGGACGTTGTGCTACGAGATGGGAGCCGGAT  >  NC_000913/3410096‑3410385

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