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

GACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACA  >  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)
          cATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGG‑cttt                                                                                                                                           <  1:1314052/150‑1 (MQ=255)
              gCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGG‑cttttcat                                                                                                                                       <  1:928969/150‑4 (MQ=255)
                gATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGATTTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTg                                                                                                                                     <  1:364431/150‑1 (MQ=255)
                gATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTg                                                                                                                                     <  1:766363/150‑1 (MQ=255)
                caTGAGGTGATCACCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTg                                                                                                                                     <  2:1146999/149‑1 (MQ=255)
                caTGAGGTGAGCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTg                                                                                                                                     <  2:1152106/149‑1 (MQ=255)
                 aTGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGt                                                                                                                                    >  1:311026/1‑150 (MQ=255)
                  tGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTc                                                                                                                                   <  2:16453/150‑1 (MQ=255)
                   gAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCt                                                                                                                                  >  1:251988/1‑150 (MQ=255)
                     tttGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCtat                                                                                                                                <  1:1094319/149‑1 (MQ=255)
                       tGAGCAGCATCCCGATGGCTGGAAGCGTTTCATCCCTGAACGTTGCGGTGGCGACCGGAATTTGCTTATTTGAAGCGGGGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAg                                                                                                                              <  1:866212/150‑1 (MQ=255)
                           cAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTgg                                                                                                                          <  1:1208214/150‑1 (MQ=255)
                               aTCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAg                                                                                                                      <  2:1296130/150‑1 (MQ=255)
                                tCCCGCTGGCTGGAAGCGTTTCTTCCCGGAACGGTTCGGTTGCGACCGGAATTTGCATATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCCGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGa                                                                                                                     <  2:1082145/150‑1 (MQ=255)
                                 cccGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACAGGAATTTGCTTATATCAAGCTGTGCGCCAGCACAGCTAATTTATCAGATTTATGAAAAGCCATCCAGATTTGGATGG‑TGTTTTTTTGTCTATAGCTGGTAAGAt                                                                                                                    >  1:197574/1‑150 (MQ=255)
                                   cGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATaa                                                                                                                  >  1:723853/1‑150 (MQ=255)
                                    gATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAAt                                                                                                                 >  1:1268888/1‑150 (MQ=255)
                                      tGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                                                               <  2:99736/150‑1 (MQ=255)
                                        gCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                                                                             <  1:511109/150‑1 (MQ=255)
                                        gCTGGAAGCGTTGCTTGCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCATAATTTTGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                                                                             <  1:1230933/150‑1 (MQ=255)
                                         cTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACAGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                                                            <  1:185625/150‑1 (MQ=255)
                                               gCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                                                                                                      >  2:1021152/1‑150 (MQ=255)
                                                 gTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGaaaa                                                                                                    >  2:764593/1‑150 (MQ=255)
                                                      ttcCCTGAACGTTTCGGTTGCGACCGTAATTTGCTTATGTGACGCGGTGCGCCAGCGCAGCTAATTTCTCATAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTc                                                                                               <  1:675904/150‑1 (MQ=255)
                                                         ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                                                                            >  1:113645/1‑150 (MQ=255)
                                                         ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGCTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCAATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                                                                            >  2:139302/1‑150 (MQ=255)
                                                          cTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGTATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                                                                           <  1:667203/150‑1 (MQ=255)
                                                                gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                                                                     >  1:163444/1‑150 (MQ=255)
                                                                  ttCGGTTGCTACCGGAAATTGCTTATTTGAAGCGGTGCGCCAGCGCATCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGAAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                                                                   <  1:978439/150‑1 (MQ=255)
                                                                   tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATTAAAAGCCATCCACTTTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAatat                                                                                  <  2:215864/150‑1 (MQ=255)
                                                                     ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAttt                                                                                >  2:722162/1‑150 (MQ=255)
                                                                          cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                                                                           >  2:837753/1‑150 (MQ=255)
                                                                           gACCGGAATTTGCTTATTTGTAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACtt                                                                          <  2:849125/150‑1 (MQ=255)
                                                                            aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                                                                         <  2:1215104/150‑1 (MQ=255)
                                                                            aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                                                                         >  1:514573/1‑150 (MQ=255)
                                                                             ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTTAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTaa                                                                        <  2:251988/150‑1 (MQ=255)
                                                                               ggAATTTGCTTATTTGCAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATGGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                                                                      <  2:772940/150‑1 (MQ=255)
                                                                                gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                                                                     >  2:855257/1‑150 (MQ=255)
                                                                                 aaTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                                                                    <  2:1157935/150‑1 (MQ=255)
                                                                                 aaTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAATTAAtata                                                                    >  1:278675/1‑150 (MQ=255)
                                                                                        ttATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATTGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTTTTAATATTTAAACTTAATCTATATGGat                                                             >  1:568917/1‑149 (MQ=255)
                                                                                        ttATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa                                                             <  1:722162/150‑1 (MQ=255)
                                                                                            ttGAAGCGGTGCGCCAGCGGAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATTG‑TTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataat                                                         <  2:588107/150‑1 (MQ=255)
                                                                                            ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataat                                                         >  1:454520/1‑150 (MQ=255)
                                                                                            ttGAAGCGGAGCGCCAGCGCAGCTAATTTCTCAGAATTATGTAAAGACATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataat                                                         <  1:567609/150‑1 (MQ=255)
                                                                                              gAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat                                                       >  2:267069/1‑150 (MQ=255)
                                                                                                  cGGTGCGACAGAGAAGCTAATTTCTGAGAATTATTAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatatta                                                   <  2:1012281/150‑3 (MQ=255)
                                                                                                   ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATTAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGACAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGAc                                                  <  2:153435/150‑1 (MQ=255)
                                                                                                     tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt                                                <  2:284243/150‑1 (MQ=255)
                                                                                                       cgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATgg                                              <  2:454520/150‑1 (MQ=255)
                                                                                                          cAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTAAATCTATATGGAATAAAATATGACATGGAAt                                           >  2:820932/1‑150 (MQ=255)
                                                                                                             cgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCg                                        <  1:492858/150‑1 (MQ=255)
                                                                                                              gcAGCTAATTTCTCATAATTATGAAAAGCCATCCAGATTTGTATGG‑TTTATTTTAGACAATAGCTTTTAAGATAATTACGTATTGCAAATAATCATTTTATTAATAATCAAACTTCAGCTATATGGAATAATATATGACATGGAATCCgc                                       <  1:187251/150‑2 (MQ=255)
                                                                                                                 tcTAATTTCTAAGAATTATGAAAATCCATCCAGATTTGGATTG‑TTTTTTTTTGTCTATAGCTGGTTATATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTgg                                    <  1:1208258/149‑1 (MQ=255)
                                                                                                                 gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTgg                                    <  2:260374/150‑1 (MQ=255)
                                                                                                                 gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTgg                                    >  2:1111686/1‑150 (MQ=255)
                                                                                                                 gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTgg                                    <  1:259131/150‑1 (MQ=255)
                                                                                                                 gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTgg                                    <  1:333179/150‑1 (MQ=255)
                                                                                                                  cTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGgcg                                  >  1:449484/1‑150 (MQ=255)
                                                                                                                             gatttATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCGGGTAAGATAATTACGTATTGCAAATATTCATTTTATTTATTTTTAAACTTAATCTATATGGATTAATATATGACATGGAATCCGTTGGCGCTAGCGActg                        >  1:65587/4‑148 (MQ=255)
                                                                                                                             gAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACgg                        >  2:688931/1‑150 (MQ=255)
                                                                                                                                caaaaaaaaACCATCCAAATCTGGATGGCTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGgcg                      >  2:635865/6‑150 (MQ=255)
                                                                                                                                caaaaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGgcgc                     >  2:1291862/6‑150 (MQ=255)
                                                                                                                                caaaaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGgcgc                     >  2:699699/6‑150 (MQ=255)
                                                                                                                                   tGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTGTGTATAGCTGGTAAGATAATTACGTATTGCAAATATTAATTTTATTAATATTTAAACTTATTCTTTATGGAATAAAAAATGACAAGGAATCCGTTGGCGCTAGCGACGGCGCTg                   >  1:507297/1‑150 (MQ=255)
                                                                                                                                    taaaaGCCATCCAGATTTGGATGGTTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGc                  <  2:311152/149‑1 (MQ=255)
                                                                                                                                       aaGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAAc              <  2:1169267/150‑1 (MQ=255)
                                                                                                                                         gCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTg            >  2:1228204/1‑150 (MQ=255)
                                                                                                                                             tCCAGATTTGTATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTAcc        <  1:1186799/150‑1 (MQ=255)
                                                                                                                                             tCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTAcc        <  2:1161213/150‑1 (MQ=255)
                                                                                                                                             tCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTAcc        <  1:837753/150‑1 (MQ=255)
                                                                                                                                                 gATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGaa    >  2:918899/1‑150 (MQ=255)
                                                                                                                                                   tttGGATGG‑TTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACa  <  1:1061877/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACA  >  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: ‑