Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 52,332 A→C I33M (ATT→ATG rsmA ← 16S rRNA m(6)A1518, m(6)A1519 dimethyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091352,3320AC96.5% 331.9 / 1.9 87I33M (ATT→ATGrsmA16S rRNA m(6)A1518, m(6)A1519 dimethyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base C (16/68);  total (16/71)
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.61e-01

GGCCTAAGAATGGATGCGTTTGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGTAATAAAA  >  NC_000913/52184‑52478
                                                                                                                                                    |                                                                                                                                                  
ggCCTAAGAATGGATGCGTTTGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCAt                                                                                                                                                  <  2:926215/151‑1 (MQ=255)
ggCCTAAGAATGGATGCGTTTGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAAt                                                                                                                                                  <  1:346781/151‑1 (MQ=255)
   cTAAGAATGGATCCGTTTGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGATGGTCGAGACGTTCGCCGGCCGGTTCGGTCAATGCCGCCTGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATgg                                                                                                                                                <  1:474082/150‑1 (MQ=255)
     aaGAATGGATGCGTTTGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGc                                                                                                                                               <  2:129060/149‑1 (MQ=255)
      aGAATGGATGCGTTTGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGc                                                                                                                                               >  2:360156/1‑148 (MQ=255)
         aTGGATGCGTTTGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGacac                                                                                                                                         >  1:574377/1‑151 (MQ=255)
         aTGGATGCGTTTGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGacac                                                                                                                                         <  1:205039/151‑1 (MQ=255)
            gATGCGTTTGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAAt                                                                                                                                      <  2:945893/151‑1 (MQ=255)
             aTGCGTTTGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATa                                                                                                                                     <  2:382046/151‑1 (MQ=255)
               gCGTTTGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACt                                                                                                                                   <  2:1297397/151‑1 (MQ=255)
                  tttGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTc                                                                                                                                >  1:1176118/1‑151 (MQ=255)
                   ttGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTc                                                                                                                                <  1:760500/150‑1 (MQ=255)
                   ttGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCg                                                                                                                               >  2:1236491/1‑151 (MQ=255)
                   ttGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTc                                                                                                                                <  1:273571/150‑1 (MQ=255)
                    tGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCg                                                                                                                               <  2:866324/150‑1 (MQ=255)
                    tGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCg                                                                                                                               <  2:692999/150‑1 (MQ=255)
                    tGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCg                                                                                                                               <  2:18998/150‑1 (MQ=255)
                    tGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCg                                                                                                                               <  1:1054357/150‑1 (MQ=255)
                       aGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCa                                                                                                                           <  2:1083317/151‑1 (MQ=255)
                       aGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCAATc                                                                                                                            >  1:425788/1‑150 (MQ=255)
                          cGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACg                                                                                                                         >  2:227578/1‑150 (MQ=255)
                          cGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGa                                                                                                                        <  2:831151/151‑1 (MQ=255)
                          cGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGa                                                                                                                        <  1:542287/151‑1 (MQ=255)
                           gTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGaa                                                                                                                       <  2:1227980/151‑1 (MQ=255)
                             gccgccAGATCGCGGTCAATTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTGTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACt                                                                                                                     <  2:1270656/151‑1 (MQ=255)
                             gccgccAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACt                                                                                                                     <  2:959622/151‑1 (MQ=255)
                             gccgccAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACt                                                                                                                     <  1:1030511/151‑1 (MQ=255)
                             gccgccAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACt                                                                                                                     <  2:715175/151‑1 (MQ=255)
                             gccgccAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACt                                                                                                                     <  1:441859/151‑1 (MQ=255)
                                gccAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGCCCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGAt                                                                                                                  >  2:45959/1‑151 (MQ=255)
                                gccAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGAt                                                                                                                  <  2:1232895/151‑1 (MQ=255)
                                 ccAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATc                                                                                                                 <  2:1232845/151‑1 (MQ=255)
                                 ccAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATc                                                                                                                 <  1:1297104/151‑1 (MQ=255)
                                  cAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCg                                                                                                                <  1:1082617/151‑1 (MQ=255)
                                           gTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGaa                                                                                                         >  2:550485/1‑149 (MQ=255)
                                             cAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGt                                                                                                      <  2:224848/150‑1 (MQ=255)
                                                     aTGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGa                                                                                             <  1:877127/151‑1 (MQ=255)
                                                       gACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAg                                                                                           <  1:928959/151‑1 (MQ=255)
                                                        aCCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGc                                                                                          <  2:726209/151‑1 (MQ=255)
                                                          cGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGt                                                                                        <  2:953350/151‑1 (MQ=255)
                                                              aGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTAc                                                                                    <  2:213220/151‑1 (MQ=255)
                                                                        gACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTg                                                                          <  1:806584/151‑1 (MQ=255)
                                                                         aCGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTgg                                                                         <  1:67205/151‑1 (MQ=255)
                                                                           gTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGcc                                                                       <  2:569516/151‑1 (MQ=255)
                                                                            ttCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGc                                                                        <  1:406328/149‑1 (MQ=255)
                                                                              cGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCtg                                                                    >  1:252192/1‑151 (MQ=255)
                                                                              cGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCt                                                                     <  2:1240176/150‑1 (MQ=255)
                                                                                ccGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCtgg                                                                   <  2:808526/150‑1 (MQ=255)
                                                                                    ccGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGAc                                                              >  1:253847/1‑151 (MQ=255)
                                                                                     cGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACt                                                             <  1:1259645/151‑1 (MQ=255)
                                                                                       gTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTc                                                            <  1:523560/150‑1 (MQ=255)
                                                                                           ggTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGatta                                                       <  1:292576/151‑1 (MQ=255)
                                                                                           ggTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGatta                                                       <  2:345762/151‑1 (MQ=255)
                                                                                           ggTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGatt                                                        >  1:71511/1‑150 (MQ=255)
                                                                                            gTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGattat                                                      <  2:148542/151‑1 (MQ=255)
                                                                                              cAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGattatt                                                     <  2:78057/150‑1 (MQ=255)
                                                                                              cAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGattatt                                                     <  2:905502/150‑1 (MQ=255)
                                                                                                aTGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCAt                                                  <  2:203870/151‑1 (MQ=255)
                                                                                                    cgccAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTggg                                              <  1:1249275/151‑1 (MQ=255)
                                                                                                    cgccAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTggg                                              <  2:930819/151‑1 (MQ=255)
                                                                                                       cAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGgtgt                                           <  1:900725/151‑1 (MQ=255)
                                                                                                            cGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACa                                      <  1:972855/151‑1 (MQ=255)
                                                                                                               ggCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATc                                   <  1:1251635/151‑1 (MQ=255)
                                                                                                                 ccGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCAt                                 <  1:1198235/151‑1 (MQ=255)
                                                                                                                   gATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCAttt                               >  2:206528/1‑151 (MQ=255)
                                                                                                                   gATTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCAttt                               <  1:866767/151‑1 (MQ=255)
                                                                                                                    aTTTCGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCAtttt                              <  1:121417/151‑1 (MQ=255)
                                                                                                                      ttCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTg                             <  1:418399/150‑1 (MQ=255)
                                                                                                                        cGACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGAt                           <  1:1216815/150‑1 (MQ=255)
                                                                                                                         gACCATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATgg                         <  1:118180/151‑1 (MQ=255)
                                                                                                                           ccATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCg                       <  2:1300681/151‑1 (MQ=255)
                                                                                                                            cATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCga                      <  1:1205811/151‑1 (MQ=255)
                                                                                                                            cATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCga                      <  1:1218256/151‑1 (MQ=255)
                                                                                                                            cATCGCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCg                       >  1:1032962/1‑150 (MQ=255)
                                                                                                                                gCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGAtt                  >  2:216632/1‑151 (MQ=255)
                                                                                                                                gCCTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGAtt                  <  1:916498/151‑1 (MQ=255)
                                                                                                                                 ccTGGCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTa                 <  2:208925/151‑1 (MQ=255)
                                                                                                                                    ggCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAg               <  2:438577/150‑1 (MQ=255)
                                                                                                                                     gCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAgcg             >  2:843442/1‑151 (MQ=255)
                                                                                                                                     gCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAgcg             <  1:589530/151‑1 (MQ=255)
                                                                                                                                     gCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAgcg             <  1:796941/151‑1 (MQ=255)
                                                                                                                                     gCCCTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAg               <  2:498406/149‑1 (MQ=255)
                                                                                                                                      cccTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCGCGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAgcgc            <  1:1193979/151‑1 (MQ=255)
                                                                                                                                       ccTTTTGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGcc           >  1:212808/1‑151 (MQ=255)
                                                                                                                                           ttGCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGta        <  2:958845/150‑1 (MQ=255)
                                                                                                                                             gCGGGTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGtaata     <  2:880835/151‑1 (MQ=255)
                                                                                                                                                ggTTCATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGTAATaaaa  <  1:493534/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGCCTAAGAATGGATGCGTTTGCAGACGTGCCGCCAGATCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCTAAGTGGCCCTGGTGGACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGTAATAAAA  >  NC_000913/52184‑52478

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