Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,410,037 T→C 8.0% intergenic (+31/‑96) rlmB → / → yjfI 23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,410,0370TC 8.0% 208.5 / 7.3 74intergenic (+31/‑96)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (26/42);  new base C (3/3);  total (29/45)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

GACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACA  >  NC_000913/4409889‑4410185
                                                                                                                                                    |                                                                                                                                                     
gACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGAttt                                                                                                                                                      >  2:582652/1‑150 (MQ=255)
 aCTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGATGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTg                                                                                                                                                     <  1:250948/150‑1 (MQ=255)
    cGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGAt                                                                                                                                                  <  2:610432/150‑1 (MQ=255)
          cATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttt                                                                                                                                           <  1:1314052/150‑1 (MQ=255)
              gCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                                                                       <  1:928969/150‑4 (MQ=255)
                gATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGATTTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTg                                                                                                                                     <  1:364431/150‑1 (MQ=255)
                gATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTg                                                                                                                                     <  1:766363/150‑1 (MQ=255)
                caTGAGGTGATCACCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTg                                                                                                                                     <  2:1146999/149‑1 (MQ=255)
                caTGAGGTGAGCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTg                                                                                                                                     <  2:1152106/149‑1 (MQ=255)
                 aTGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGt                                                                                                                                    >  1:311026/1‑150 (MQ=255)
                  tGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTc                                                                                                                                   <  2:16453/150‑1 (MQ=255)
                   gAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCt                                                                                                                                  >  1:251988/1‑150 (MQ=255)
                     tttGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCtat                                                                                                                                <  1:1094319/149‑1 (MQ=255)
                       tGAGCAGCATCCCGATGGCTGGAAGCGTTTCATCCCTGAACGTTGCGGTGGCGACCGGAATTTGCTTATTTGAAGCGGGGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAg                                                                                                                              <  1:866212/150‑1 (MQ=255)
                           cAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTgg                                                                                                                          <  1:1208214/150‑1 (MQ=255)
                               aTCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAg                                                                                                                      <  2:1296130/150‑1 (MQ=255)
                                tCCCGCTGGCTGGAAGCGTTTCTTCCCGGAACGGTTCGGTTGCGACCGGAATTTGCATATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCCGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGa                                                                                                                     <  2:1082145/150‑1 (MQ=255)
                                 cccGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACAGGAATTTGCTTATATCAAGCTGTGCGCCAGCACAGCTAATTTATCAGATTTATGAAAAGCCATCCAGATTTGGATGGTGTTTTTTTGTCTATAGCTGGTAAGAt                                                                                                                    >  1:197574/1‑150 (MQ=255)
                                   cGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATaa                                                                                                                  >  1:723853/1‑150 (MQ=255)
                                    gATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAAt                                                                                                                 >  1:1268888/1‑150 (MQ=255)
                                      tGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                                                               <  2:99736/150‑1 (MQ=255)
                                        gCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                                                                             <  1:511109/150‑1 (MQ=255)
                                        gCTGGAAGCGTTGCTTGCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCATAATTTTGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                                                                             <  1:1230933/150‑1 (MQ=255)
                                         cTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACAGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                                                            <  1:185625/150‑1 (MQ=255)
                                               gCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                                                                                                      >  2:1021152/1‑150 (MQ=255)
                                                 gTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGaaaa                                                                                                    >  2:764593/1‑150 (MQ=255)
                                                      ttcCCTGAACGTTTCGGTTGCGACCGTAATTTGCTTATGTGACGCGGTGCGCCAGCGCAGCTAATTTCTCATAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTc                                                                                               <  1:675904/150‑1 (MQ=255)
                                                         ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                                                                            >  1:113645/1‑150 (MQ=255)
                                                         ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGCTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCAATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                                                                            >  2:139302/1‑150 (MQ=255)
                                                          cTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGTATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                                                                           <  1:667203/150‑1 (MQ=255)
                                                                gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                                                                     >  1:163444/1‑150 (MQ=255)
                                                                  ttCGGTTGCTACCGGAAATTGCTTATTTGAAGCGGTGCGCCAGCGCATCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGAAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                                                                   <  1:978439/150‑1 (MQ=255)
                                                                   tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATTAAAAGCCATCCACTTTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAatat                                                                                  <  2:215864/150‑1 (MQ=255)
                                                                     ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAttt                                                                                >  2:722162/1‑150 (MQ=255)
                                                                          cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                                                                           >  2:837753/1‑150 (MQ=255)
                                                                           gACCGGAATTTGCTTATTTGTAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACtt                                                                          <  2:849125/150‑1 (MQ=255)
                                                                            aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                                                                         <  2:1215104/150‑1 (MQ=255)
                                                                            aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                                                                         >  1:514573/1‑150 (MQ=255)
                                                                             ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTTAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTaa                                                                        <  2:251988/150‑1 (MQ=255)
                                                                               ggAATTTGCTTATTTGCAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATGGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                                                                      <  2:772940/150‑1 (MQ=255)
                                                                                gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                                                                     >  2:855257/1‑150 (MQ=255)
                                                                                 aaTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                                                                    <  2:1157935/150‑1 (MQ=255)
                                                                                 aaTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAATTAAtata                                                                    >  1:278675/1‑150 (MQ=255)
                                                                                        ttATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATTGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTTTTAATATTTAAACTTAATCTATATGGat                                                             >  1:568917/1‑149 (MQ=255)
                                                                                        ttATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa                                                             <  1:722162/150‑1 (MQ=255)
                                                                                            ttGAAGCGGTGCGCCAGCGGAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATTGTTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataat                                                         <  2:588107/150‑1 (MQ=255)
                                                                                            ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataat                                                         >  1:454520/1‑150 (MQ=255)
                                                                                            ttGAAGCGGAGCGCCAGCGCAGCTAATTTCTCAGAATTATGTAAAGACATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataat                                                         <  1:567609/150‑1 (MQ=255)
                                                                                              gAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat                                                       >  2:267069/1‑150 (MQ=255)
                                                                                                  cGGTGCGACAGAGAAGCTAATTTCTGAGAATTATTAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatatta                                                   <  2:1012281/150‑3 (MQ=255)
                                                                                                   ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATTAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGACAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGAc                                                  <  2:153435/150‑1 (MQ=255)
                                                                                                     tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt                                                <  2:284243/150‑1 (MQ=255)
                                                                                                       cgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATgg                                              <  2:454520/150‑1 (MQ=255)
                                                                                                          cAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTAAATCTATATGGAATAAAATATGACATGGAAt                                           >  2:820932/1‑150 (MQ=255)
                                                                                                             cgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCg                                        <  1:492858/150‑1 (MQ=255)
                                                                                                              gcAGCTAATTTCTCATAATTATGAAAAGCCATCCAGATTTGTATGGTTTATTTTAGACAATAGCTTTTAAGATAATTACGTATTGCAAATAATCATTTTATTAATAATCAAACTTCAGCTATATGGAATAATATATGACATGGAATCCgc                                       <  1:187251/150‑2 (MQ=255)
                                                                                                                 tcTAATTTCTAAGAATTATGAAAATCCATCCAGATTTGGATTGTTTTTTTTTGTCTATAGCTGGTTATATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTgg                                    <  1:1208258/149‑1 (MQ=255)
                                                                                                                 gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTgg                                    <  2:260374/150‑1 (MQ=255)
                                                                                                                 gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTgg                                    >  2:1111686/1‑150 (MQ=255)
                                                                                                                 gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTgg                                    <  1:259131/150‑1 (MQ=255)
                                                                                                                 gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTgg                                    <  1:333179/150‑1 (MQ=255)
                                                                                                                  cTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGgcg                                  >  1:449484/1‑150 (MQ=255)
                                                                                                                             gatttATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCGGGTAAGATAATTACGTATTGCAAATATTCATTTTATTTATTTTTAAACTTAATCTATATGGATTAATATATGACATGGAATCCGTTGGCGCTAGCGActg                        >  1:65587/4‑148 (MQ=255)
                                                                                                                             gAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACgg                        >  2:688931/1‑150 (MQ=255)
                                                                                                                                caaaaaaaaACCATCCAAATCTGGATGGCTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGgcg                      >  2:635865/6‑150 (MQ=255)
                                                                                                                                caaaaaaaaACCATCCAAATCTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGgcgc                     >  2:1291862/6‑150 (MQ=255)
                                                                                                                                caaaaaaaaACCATCCAAATCTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGgcgc                     >  2:699699/6‑150 (MQ=255)
                                                                                                                                   tGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTGTGTATAGCTGGTAAGATAATTACGTATTGCAAATATTAATTTTATTAATATTTAAACTTATTCTTTATGGAATAAAAAATGACAAGGAATCCGTTGGCGCTAGCGACGGCGCTg                   >  1:507297/1‑150 (MQ=255)
                                                                                                                                    taaaaGCCATCCAGATTTGGATGGTTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGc                  <  2:311152/149‑1 (MQ=255)
                                                                                                                                       aaGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAAc              <  2:1169267/150‑1 (MQ=255)
                                                                                                                                         gCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTg            >  2:1228204/1‑150 (MQ=255)
                                                                                                                                             tCCAGATTTGTATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTAcc        <  1:1186799/150‑1 (MQ=255)
                                                                                                                                             tCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTAcc        <  2:1161213/150‑1 (MQ=255)
                                                                                                                                             tCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTAcc        <  1:837753/150‑1 (MQ=255)
                                                                                                                                                 gATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGaa    >  2:918899/1‑150 (MQ=255)
                                                                                                                                                   tttGGATGGTTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACa  <  1:1061877/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACA  >  NC_000913/4409889‑4410185

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: