Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 90,576 C→A 5.3% A161A (GCC→GCA rsmH → 16S rRNA m(4)C1402 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091390,5760CA 5.3% 526.3 / 15.8 151A161A (GCC→GCArsmH16S rRNA m(4)C1402 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (75/68);  new base A (4/4);  total (79/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.25e-02

Only 100 of 159 total aligned reads displayed.

TGAACGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCAT  >  NC_000913/90429‑90722
                                                                                                                                                   |                                                                                                                                                  
tGAACGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCa                                                                                                                                                   >  3:471441/1‑149 (MQ=255)
   aCGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAAc                                                                                                                                                >  1:408316/1‑149 (MQ=255)
     gTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGc                                                                                                                                              >  3:1271081/1‑149 (MQ=255)
      tGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCa                                                                                                                                             <  4:352531/149‑1 (MQ=255)
              ccTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCcgc                                                                                                                                     >  4:245258/1‑149 (MQ=255)
                          aTGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAg                                                                                                                         <  2:1283531/149‑1 (MQ=255)
                             gTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGt                                                                                                                      <  1:909107/149‑1 (MQ=255)
                              tCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTa                                                                                                                     >  1:657796/1‑149 (MQ=255)
                                cGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCaa                                                                                                                                                  <  3:391935/118‑1 (MQ=255)
                                cGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCaa                                                                                                                                                  >  4:391935/1‑118 (MQ=255)
                                      ataTGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAGGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGaa                                                                                                             <  4:1257721/147‑1 (MQ=255)
                                      aCATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGaa                                                                                                             <  2:51147/149‑1 (MQ=255)
                                         tGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAg                                                                                                          <  2:847152/149‑1 (MQ=255)
                                          gCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGc                                                                                                         <  1:1021222/149‑1 (MQ=255)
                                             tATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGa                                                                                                      <  2:124441/149‑1 (MQ=255)
                                                 gACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCAt                                                                                                                                <  1:1208761/119‑1 (MQ=255)
                                                  aCCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGGAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGAcc                                                                                                 >  2:125111/1‑149 (MQ=255)
                                                  aCCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGc                                                                                                                                         >  1:1107132/1‑109 (MQ=255)
                                                    ccAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCg                                                                                               >  3:51238/1‑149 (MQ=255)
                                                    ccAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCg                                                                                               >  2:1095229/1‑149 (MQ=255)
                                                    ccAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTGGCAAAACGctc                                                                                                                                            <  1:33761/104‑3 (MQ=255)
                                                    ccAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTGGCAAAACGctc                                                                                                                                            >  2:33761/1‑102 (MQ=255)
                                                      aaCCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGa                                                                                                                          >  1:259610/1‑120 (MQ=255)
                                                      aaCCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGa                                                                                                                          <  2:259610/120‑1 (MQ=255)
                                                      aaCCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGt                                                                                                                      >  4:625411/1‑124 (MQ=255)
                                                           gTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCaaa                                                                                        <  2:657796/149‑1 (MQ=255)
                                                           gTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCaaa                                                                                        <  4:537385/149‑1 (MQ=255)
                                                            tGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAg                                                                                       >  1:1066340/1‑149 (MQ=255)
                                                                      gCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCaaa                                                                                        >  2:1415701/1‑138 (MQ=255)
                                                                      gCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGa                                                                             >  3:74822/1‑149 (MQ=255)
                                                                       ccGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGaa                                                                            <  4:1010294/149‑1 (MQ=255)
                                                                          cTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTc                                                                         >  4:891969/1‑149 (MQ=255)
                                                                                ggCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGgct                                                                   >  2:218651/1‑149 (MQ=255)
                                                                                   tACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAAc                                                                                    <  3:86817/129‑1 (MQ=255)
                                                                                   tACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAAc                                                                                    >  4:86817/1‑129 (MQ=255)
                                                                                   tACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGgctgct                                                                >  3:228495/1‑149 (MQ=255)
                                                                                    aCAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCaaa                                                                                        <  4:661750/124‑1 (MQ=255)
                                                                                    aCAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCaaa                                                                                        >  3:661750/1‑124 (MQ=255)
                                                                                    aCAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGaa                                                                                     <  3:1254283/127‑1 (MQ=255)
                                                                                    aCAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGaa                                                                                     >  4:1254283/1‑127 (MQ=255)
                                                                                       aaCCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTgg                                                                                 >  3:53140/1‑128 (MQ=255)
                                                                                       aaCCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCACTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTgg                                                                                 <  4:53140/128‑1 (MQ=255)
                                                                                              gaagaaGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGaa                                                                                     <  3:1126905/117‑1 (MQ=255)
                                                                                               aagaagCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGAc                                                                                                  >  2:1331052/1‑103 (MQ=255)
                                                                                               aagaagCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGAc                                                                                                  <  1:1331052/103‑1 (MQ=255)
                                                                                                 gaagCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGaa                                                  >  1:371628/1‑149 (MQ=255)
                                                                                                  aagCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGaaa                                                 <  4:74822/149‑1 (MQ=255)
                                                                                                  aagCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGaaa                                                 <  1:1275140/149‑1 (MQ=255)
                                                                                                   agCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAg                                                <  3:891969/149‑1 (MQ=255)
                                                                                                    gCCGATATCGTCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAAct                                                                                                           >  4:372830/1‑88 (MQ=255)
                                                                                                    gCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACa                                                                                                           <  3:372830/89‑1 (MQ=255)
                                                                                                    gCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCaaa                                                                                        <  3:28287/108‑1 (MQ=255)
                                                                                                    gCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGa                                               >  3:31229/1‑149 (MQ=255)
                                                                                                     ccGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGaaa                                                 >  1:1399976/1‑146 (MQ=255)
                                                                                                     ccGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGAt                                              >  4:443629/1‑149 (MQ=255)
                                                                                                       gATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATaa                                            <  4:1273612/149‑1 (MQ=255)
                                                                                                         tatCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGt                                          <  2:408316/149‑1 (MQ=255)
                                                                                                             gCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGa                                                                                                      <  1:567891/85‑1 (MQ=255)
                                                                                                             gCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTaa                                      >  4:491171/1‑149 (MQ=255)
                                                                                                              ccTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTaaa                                     <  2:1066340/149‑1 (MQ=255)
                                                                                                              ccTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTaaa                                     <  3:43411/149‑1 (MQ=255)
                                                                                                              ccTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTaaa                                     <  1:1095229/149‑1 (MQ=255)
                                                                                                              ccTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTaaa                                     <  4:228495/149‑1 (MQ=255)
                                                                                                               cTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCaaa                                                                                        >  1:315908/1‑97 (MQ=255)
                                                                                                               cTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCaaa                                                                                        <  2:315908/97‑1 (MQ=255)
                                                                                                               cTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCGCGCCATTGTCGa                                                                                                                          <  2:546834/62‑1 (MQ=255)
                                                                                                                   gTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATcc                                <  2:371628/149‑1 (MQ=255)
                                                                                                                    tATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGaaa                                                 <  1:852101/131‑1 (MQ=255)
                                                                                                                     aTTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCcg                              >  4:744946/1‑149 (MQ=255)
                                                                                                                      ttGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTc                                                                         >  2:876771/1‑105 (MQ=255)
                                                                                                                        gAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCg                                                       <  3:626325/121‑1 (MQ=255)
                                                                                                                        gAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGa                           >  2:503477/1‑149 (MQ=255)
                                                                                                                         aaaaCCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATa                                             >  1:336184/1‑130 (MQ=255)
                                                                                                                          aaaCCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGaa                                                                            <  3:573643/98‑1 (MQ=255)
                                                                                                                           aaCCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGAccc                        <  1:153517/149‑1 (MQ=255)
                                                                                                                              cTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACa                                                                                                           <  3:1054964/63‑1 (MQ=255)
                                                                                                                              cTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACa                                                                                                           >  4:1054964/1‑63 (MQ=255)
                                                                                                                                   gTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGt                                                                       <  1:511950/94‑1 (MQ=255)
                                                                                                                                   gTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTc                >  1:458458/1‑149 (MQ=255)
                                                                                                                                       aatgCGTTTGGCAAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGaa                                                                            <  2:929174/82‑1 (MQ=255)
                                                                                                                                       aatgCGTTTGGCAAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGaa                                                                            >  1:929174/4‑85 (MQ=255)
                                                                                                                                       aatgCGTTTGGCAAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTa                     <  1:1205410/137‑1 (MQ=255)
                                                                                                                                       aatgCGTTTGGCAAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTa                     >  2:1205410/4‑140 (MQ=255)
                                                                                                                                        gagCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGc           >  2:325835/1‑149 (MQ=255)
                                                                                                                                           cGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGACCAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGt        <  3:1240043/149‑1 (MQ=255)
                                                                                                                                           cGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGt                                          <  4:144524/115‑1 (MQ=255)
                                                                                                                                           cGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGt                                          >  3:144524/1‑115 (MQ=255)
                                                                                                                                           cGTTTGGCAAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTgg                                                                     <  3:578718/88‑1 (MQ=255)
                                                                                                                                           cGTTTGGCAAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTgg                                                                     >  4:578718/1‑88 (MQ=255)
                                                                                                                                            gTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTg       <  3:491171/149‑1 (MQ=255)
                                                                                                                                             ttttGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGt                                                                       >  2:537979/1‑84 (MQ=255)
                                                                                                                                             ttttGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGt                                                                       <  1:537979/84‑1 (MQ=255)
                                                                                                                                                 gCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACGTCCCGCGACCCGTACCTTCCAGGCGGTGCGCAt  >  4:965724/1‑149 (MQ=255)
                                                                                                                                                  ccAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGgctgct                                                                >  2:761714/1‑86 (MQ=255)
                                                                                                                                                  ccAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGgctgct                                                                <  1:761714/86‑1 (MQ=255)
                                                                                                                                                   cAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGc                                                                                <  4:448124/69‑1 (MQ=255)
                                                                                                                                                   cAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGc                                                                                >  3:448124/1‑69 (MQ=255)
                                                                                                                                                   cAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCgg         <  4:1189250/140‑1 (MQ=255)
                                                                                                                                                   cAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCgg         >  1:900618/1‑140 (MQ=255)
                                                                                                                                                   cAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCgg         <  2:900618/140‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGAACGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCAT  >  NC_000913/90429‑90722

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