Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,371,401 C→T A126V (GCT→GTT)  groL → chaperonin GroEL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,371,4010CT98.4% 208.9 / NA 62A126V (GCT→GTT) groLchaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (29/32);  minor base G (0/1);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TGCACCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAG  >  NC_000913/4371159‑4371645
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
tGCACCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGCCAAAGCGGTTACCGTTGCAGTTg                                                                                                                                                                                                                                              >  1:17614/1‑251 (MQ=255)
  cACCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTgaa                                                                                                                                                                                                                                            <  2:14626/251‑1 (MQ=255)
   aCCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTga                                                                                                                                                                                                                                             >  2:298325/1‑249 (MQ=255)
        ccATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTg                                                                                                                                                                                                                                      >  1:61730/1‑251 (MQ=255)
        ccATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTg                                                                                                                                                                                                                                      <  1:326559/251‑1 (MQ=255)
              ccAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAGACGACGCTGCAGGCGACGGTACCATCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAGAGCTGTTGCTGGGGGCATGAACCCGATGGACCTGAAACGGGGTATCGACAGAGCGGTTACCGTTGCAGTTGAAGAACTGAAAggg                                                                                                                                                                                                                                >  1:389448/1‑249 (MQ=255)
              ccAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGAGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAgcg                                                                                                                                                                                                                                >  2:370245/1‑251 (MQ=255)
                 aaGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTTAAAGCGCTg                                                                                                                                                                                                                             <  2:78027/251‑1 (MQ=255)
                                  gCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGACAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCACTg                                                                                                                                                                                                             >  2:56791/1‑250 (MQ=255)
                                   cTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGAc                                                                                                                                                                                                           >  1:125650/1‑251 (MQ=255)
                                   cTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGAc                                                                                                                                                                                                           >  1:302464/1‑251 (MQ=255)
                                      gTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGActct                                                                                                                                                                                                        <  1:46791/251‑1 (MQ=255)
                                      gTGAAATCGAACGGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGActct                                                                                                                                                                                                        <  2:383018/251‑1 (MQ=255)
                                           aTCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTAACACCACTGCAACCGTACTGGCTCCGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGGCAAAGCGGTTACCGTTGCAGTTGCAGAACTGAAAGCGCTGTCCGTAcc                                                                                                                                                                                                                     >  2:94278/1‑233 (MQ=255)
                                           aTCGAACTGGAAGAAAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTAcc                                                                                                                                                                                                                     <  1:94278/233‑1 (MQ=255)
                                                      aGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCGGACTCTAAAGCGATTGCTCAgg                                                                                                                                                                                        <  1:127879/251‑1 (MQ=255)
                                                         cAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACTCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGATCTGAAAGCGCTGTCCGTACCATGCGCTGACTCTAAAGCGATTGCTCAGGTTg                                                                                                                                                                                     <  2:28081/251‑1 (MQ=255)
                                                          aaGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCTCCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTgg                                                                                                                                                                                    <  1:404095/251‑1 (MQ=255)
                                                                     tatGGGTGCGCAGCTGGTGCAAGAAGTTGACTCTAAAGCACACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATAATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGAGATTGCTCAGGTTGGTACGATCTCCg                                                                                                                                                                         >  2:140960/1‑251 (MQ=255)
                                                                        gggTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGATCTGAAAGCGCTGTCCGTGCCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCt                                                                                                                                                                       <  1:51149/250‑1 (MQ=255)
                                                                              gcAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCGGCGGGCGACGGTAACACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCGGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCGGACTCTAAAGCGAGTGCGCAGGTTGGTACCATCGCCGCTAACTCcg                                                                                                                                                                <  2:91058/251‑1 (MQ=255)
                                                                                aGATGGTGAAAGAAGTTGCCTCTACAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCcgac                                                                                                                                                              >  2:198132/1‑251 (MQ=255)
                                                                                      tGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAAcc                                                                                                                                                        >  2:385439/1‑251 (MQ=255)
                                                                                        aaaGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGt                                                                                                                                                      <  1:86855/251‑1 (MQ=255)
                                                                                         aagaagTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGt                                                                                                                                                      >  2:189515/1‑250 (MQ=255)
                                                                                             atttGCCTCTAATGCAAGCGACGCTGCATGCGACGGTACCTCCACTGCTACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGATCCCGATGGACCTTAAACGTGGTATCGACTAAGCGGTTACCGTTGCAGTTGAATAACTGAAAGCGCTGTCCGTTCCATGCTCTGACTCTAAAGCGATTGCTCATGTTGTTACTATcgccgcggaatccttcgagaccttaggta                                                                                                                                                 <  2:264764/249‑24 (MQ=255)
                                                                                                  ccTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGAACCATGCTCTGACTCTAAAGCGATTGCTCAGGATGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTg                                                                                                                                            >  1:174336/1‑251 (MQ=255)
                                                                                                  ccTCTAAAGCAAACGACACTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCGCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTg                                                                                                                                            <  2:307461/251‑1 (MQ=255)
                                                                                                           cAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAGCCGTCGGTAAACTGATccatga                                                                                                                                    >  2:138962/1‑245 (MQ=255)
                                                                                                           cAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGaa                                                                                                                                   <  1:385788/251‑1 (MQ=255)
                                                                                                           cAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGATCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGCAGAACTGCAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGa                                                                                                                                    >  2:31310/1‑250 (MQ=255)
                                                                                                               cgacgCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTAGCAGCACTGAAGGTCTGAAAGCTGTTGCGGCGGGCATGAACCCGATGGACCTGAAACGTGGGATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCAGGCTCTGACTCTAAAGCGAGTGCTCAGGTTGGTACCAGCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGa                                                                                                                               <  2:204798/251‑1 (MQ=255)
                                                                                                                gacgCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGAt                                                                                                                              >  2:160551/1‑251 (MQ=255)
                                                                                                                  cgCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATgg                                                                                                                            <  1:312871/251‑1 (MQ=255)
                                                                                                                              cGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTa                                                                                                                <  1:91595/251‑1 (MQ=255)
                                                                                                                                   ccaccaCTGCAACCGTACTGGCTCAGGCTATCATCACTGAATGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGCTCGCTGAAGCGATGTACAAAGTCGGTAaagaa                                                                                                           <  2:293192/251‑1 (MQ=255)
                                                                                                                                      ccacTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTCCCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGc                                                                                                        >  1:347675/1‑251 (MQ=255)
                                                                                                                                      ccacTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAggg                                                                                                        >  1:51518/1‑250 (MQ=255)
                                                                                                                                           gCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTa                                                                                                    <  2:253896/250‑1 (MQ=255)
                                                                                                                                             aaCCGTACTGGATCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGCCAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGATGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATc                                                                                                  >  2:351811/1‑250 (MQ=255)
                                                                                                                                                 gTACTGGCTCAGGCTATCATCACTGAAGGGCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGt                                                                                             >  2:97617/1‑251 (MQ=255)
                                                                                                                                                     tGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGGTGGTACCATCTCCGCTAACTACGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGa                                                                                          >  2:1405/1‑250 (MQ=255)
                                                                                                                                                       gCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATTGTCCTGAAACGTGGTATCGACAAAGCGGTTACCTGTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTTCTCAGGTTGGTACTATCTCCGCTAACTCCTACGAAACCTTAGGTAAACTGATCGCTGAAGCGATGGACAAAGACGGTAAATAATTCGTTATCACCGGTGAAGa                                                                                       <  1:319089/251‑1 (MQ=255)
                                                                                                                                                                         gAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGa                                                                     <  1:342121/251‑1 (MQ=255)
                                                                                                                                                                          aaGGTCTGAAAGCTGTTGCTGCGGGCATGAACCTGCTGGACCTGAACCGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGGACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGAc                                                                    >  2:65303/1‑251 (MQ=255)
                                                                                                                                                                                 gAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTgg                                                             >  2:190282/1‑251 (MQ=255)
                                                                                                                                                                                    aGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCTGTCTGCAGGACTAACTGGACg                                                          <  1:326033/251‑1 (MQ=255)
                                                                                                                                                                                     gCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGt                                                         <  1:190282/251‑1 (MQ=255)
                                                                                                                                                                                          tgctgcGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGtt                                                     <  1:328643/250‑1 (MQ=255)
                                                                                                                                                                                           gctgctGGCATGATCCCGATGGGCTTGAAACGTGGTATCGACAAAGCGGTGACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTTCCATGCTCTGACTCTAAAGCGATTGCTCAGGTTTGTACCATCTCCGCGATCTCCTACGAAACCGTAGGTAAACTTACCGCTGAAGGGATGTACAGAGTCGGTAAAGAAGTCGTGCTCGCCGGTGAAGACTGTACCGGTCTGTGGGACGAGATGGACGTGGTTg                                                    <  2:401028/250‑1 (MQ=255)
                                                                                                                                                                                                gggCATGAACCCGATGGACCTGAAACGTGGCATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCACTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGt                                               >  1:114591/1‑250 (MQ=255)
                                                                                                                                                                                                  gCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGACGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGCAGACGGTACCGGTCTGCAGGACGAACTGGACGTGTTTGAAGTTATg                                            >  2:144022/1‑251 (MQ=255)
                                                                                                                                                                                                      gAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAg                                         <  1:198132/250‑1 (MQ=255)
                                                                                                                                                                                                      gAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAg                                         <  1:351811/250‑1 (MQ=255)
                                                                                                                                                                                                          ccGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGCAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCg                                     >  1:379240/1‑250 (MQ=255)
                                                                                                                                                                                                           cGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCg                                                                                                                                                                         >  1:323270/1‑117 (MQ=255)
                                                                                                                                                                                                           cGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAGCTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGAc                                   >  1:70276/1‑251 (MQ=255)
                                                                                                                                                                                                           cGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGGTGGTACCATCTCCg                                                                                                                                                                         <  2:323270/117‑1 (MQ=255)
                                                                                                                                                                                                                        aCGTGGTATCGACAAAGCGGGTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGt                      <  2:51518/251‑1 (MQ=255)
                                                                                                                                                                                                                         cGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGtc                     <  1:370245/251‑1 (MQ=255)
                                                                                                                                                                                                                          gTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGtct                    <  1:87896/251‑1 (MQ=255)
                                                                                                                                                                                                                                             tACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAg  <  1:226793/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
TGCACCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAG  >  NC_000913/4371159‑4371645

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