Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 892,457 T→A 9.0% I164N (ATT→AAT)  rimK → ribosomal protein S6 modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913892,4570TA 9.0% 248.7 / 12.1 78I164N (ATT→AAT) rimKribosomal protein S6 modification protein
Reads supporting (aligned to +/- strand):  ref base T (31/39);  new base A (2/5);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAA  >  NC_000913/892309‑892604
                                                                                                                                                     |                                                                                                                                                   
cctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGAt                                                                                                                                                     <  4:410173/149‑1 (MQ=255)
cctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGAt                                                                                                                                                     <  3:103812/149‑1 (MQ=255)
 ctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGAtt                                                                                                                                                    <  3:156961/149‑1 (MQ=255)
 ctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGAtt                                                                                                                                                    <  3:336634/149‑1 (MQ=255)
   gTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTg                                                                                                                                                   <  2:543504/148‑1 (MQ=255)
    tCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGa                                                                                                                                                  <  1:26445/148‑1 (MQ=255)
        gggCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCtt                                                                                                                                             >  4:140725/1‑149 (MQ=255)
         ggCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTTCAGACCGCGGAAAGCGTCAATCACGCtt                                                                                                                                             <  4:770542/149‑1 (MQ=255)
         ggCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTCAATCACGCttt                                                                                                                                            <  4:182958/149‑1 (MQ=255)
              tGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCcgcg                                                                                                                                       <  4:849306/149‑1 (MQ=255)
                cgcATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGt                                                                                                                                     <  4:733835/149‑1 (MQ=255)
                 gcATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGACCGCGGAAAGCGTCAATCACGCTTTCcgccgc                                                                                                                                     <  1:17277/148‑4 (MQ=39)
                  cATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTc                                                                                                                                    >  1:804274/1‑148 (MQ=255)
                          ggATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACgcgc                                                                                                                           <  3:208229/149‑1 (MQ=255)
                               atACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATAtt                                                                                                                      <  3:472848/149‑1 (MQ=255)
                                   cAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTgg                                                                                                                  >  4:298312/1‑149 (MQ=255)
                                     gCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGt                                                                                                                 >  2:336704/1‑148 (MQ=255)
                                     gCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGt                                                                                                                 >  2:805947/1‑148 (MQ=255)
                                      cGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTg                                                                                                                >  2:112194/1‑148 (MQ=255)
                                       gATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGc                                                                                                               <  1:79717/148‑1 (MQ=255)
                                       gATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCa                                                                                                              <  3:841082/149‑1 (MQ=255)
                                         tttAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAg                                                                                                             >  2:359144/1‑148 (MQ=255)
                                            aaTCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGaa                                                                                                          <  1:765216/148‑1 (MQ=255)
                                             aTCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGATTTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAat                                                                                                         >  2:146683/1‑148 (MQ=255)
                                             aTCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAAcgcg                                                                                                                            <  1:630979/129‑1 (MQ=255)
                                             aTCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAAcgcg                                                                                                                            >  2:630979/1‑129 (MQ=255)
                                             aTCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAata                                                                                                        <  3:569168/149‑1 (MQ=255)
                                              tCGACATGGTCGGTGGTGCACCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTCAATCACGCtt                                                                                                                                             >  1:722685/1‑111 (MQ=255)
                                              tCGACATGGTCGGTGGTGCACCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTCAATCACGCtt                                                                                                                                             <  2:722685/111‑1 (MQ=255)
                                                  cATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCa                                                                                                   >  3:30720/1‑149 (MQ=255)
                                                    tGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCaaa                                                                                                 <  4:799330/149‑1 (MQ=255)
                                                    tGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCaaa                                                                                                 <  3:326223/149‑1 (MQ=255)
                                                    tGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCaa                                                                                                  <  1:203021/148‑1 (MQ=255)
                                                      gTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGTAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCCGAACGCGCATATTCTGGTGCAGGAATATATCAAaga                                                                                               <  4:807069/149‑1 (MQ=255)
                                                      gTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAag                                                                                                <  1:836692/148‑1 (MQ=255)
                                                                 gccgcTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGt                                                                                    >  3:68464/1‑149 (MQ=255)
                                                                           gtCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCg                                                                          <  4:724202/149‑1 (MQ=255)
                                                                             cAAGGTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGc                                                                         <  1:753475/148‑1 (MQ=255)
                                                                                gttggttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGt                                                                      <  1:336704/148‑1 (MQ=255)
                                                                                 ttggttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCt                                                                    >  4:148356/1‑149 (MQ=255)
                                                                                 ttggttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCt                                                                    >  3:162966/1‑149 (MQ=255)
                                                                                    gttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgt                                                                 <  3:522799/149‑1 (MQ=255)
                                                                                     ttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgt                                                                 <  1:112194/148‑1 (MQ=255)
                                                                                      tgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGGGGAAAGCGTGAGTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgtt                                                                >  2:828352/1‑148 (MQ=255)
                                                                                      tgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCt                                                                                                                    <  2:264166/96‑1 (MQ=255)
                                                                                      tgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCt                                                                                                                    >  1:264166/1‑96 (MQ=255)
                                                                                       gAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgttg                                                               >  1:128568/1‑148 (MQ=255)
                                                                                        aaGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgttgtt                                                             <  4:30720/149‑1 (MQ=255)
                                                                                        aaGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgttgtt                                                             <  4:68464/149‑1 (MQ=255)
                                                                                          ggCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGttgttg                                                            >  1:33929/1‑148 (MQ=255)
                                                                                           gCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTgg                                                           >  1:943550/1‑148 (MQ=255)
                                                                                               gCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATg                                                      >  3:822719/1‑149 (MQ=255)
                                                                                                 aGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAaga                                                                                               >  3:170511/1‑106 (MQ=255)
                                                                                                 aGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAaga                                                                                               <  4:170511/106‑1 (MQ=255)
                                                                                                  gggAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGaa                                                    <  1:134194/148‑1 (MQ=255)
                                                                                                       ttGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGt                                               >  1:888539/1‑148 (MQ=255)
                                                                                                           cgtcgtGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCt                                           >  1:861064/1‑148 (MQ=255)
                                                                                                             tcgtGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGc                                         >  2:534204/1‑148 (MQ=255)
                                                                                                             tcgtGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCg                                        >  3:544275/1‑149 (MQ=255)
                                                                                                                  cTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTg                                    <  2:458941/148‑1 (MQ=255)
                                                                                                                   tGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGaa                                  <  4:268784/149‑1 (MQ=255)
                                                                                                                   tGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGa                                   <  1:453066/148‑1 (MQ=255)
                                                                                                                    ggCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGaa                                  >  1:661642/1‑148 (MQ=255)
                                                                                                                         agaCGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAAcggcg                             <  1:932888/148‑1 (MQ=255)
                                                                                                                           aCGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACggcgggc                          <  3:127648/149‑1 (MQ=255)
                                                                                                                               g‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGTGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGaaa                      <  4:646934/149‑1 (MQ=255)
                                                                                                                                 tCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGCCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGaaa                      >  2:695863/1‑148 (MQ=255)
                                                                                                                                   aGGCGGCGGAAAGCGTGATTGACGCTTTGCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAaga                    <  1:519898/148‑1 (MQ=255)
                                                                                                                                    ggcggcGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAagag                   >  2:877376/1‑148 (MQ=255)
                                                                                                                                    ggcggcGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAgg                  <  3:761910/149‑1 (MQ=255)
                                                                                                                                     gcggcgGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAggg                 >  3:431141/1‑149 (MQ=255)
                                                                                                                                     gcggcgGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAgg                  <  2:47445/148‑1 (MQ=255)
                                                                                                                                            aaaGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGAttt           >  1:909448/1‑148 (MQ=255)
                                                                                                                                             aaGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgtt                                                                >  4:426901/1‑94 (MQ=255)
                                                                                                                                             aaGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgtt                                                                <  3:426901/94‑1 (MQ=255)
                                                                                                                                             aaGCGTCAATCACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGt                                                  <  2:770168/108‑1 (MQ=255)
                                                                                                                                             aaGCGTCAATCACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGt                                                  >  1:770168/1‑108 (MQ=255)
                                                                                                                                               gCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGt       <  4:481469/149‑1 (MQ=255)
                                                                                                                                                cGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGtt      <  4:544275/149‑1 (MQ=255)
                                                                                                                                                cGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGt       <  2:42515/148‑1 (MQ=255)
                                                                                                                                                   gATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCa   >  4:941410/1‑149 (MQ=255)
                                                                                                                                                    aTTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCaa  >  4:270735/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCG‑TCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAA  >  NC_000913/892309‑892604

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