Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 644,420 A→C 17.6% V195G (GTG→GGG)  rna ← ribonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913644,4200AC17.6% 115.2 / 2.0 62V195G (GTG→GGG) rnaribonuclease I
Reads supporting (aligned to +/- strand):  ref base A (16/33);  new base C (10/3);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 9.39e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCG  >  NC_000913/644192‑644639
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                            
gtgttttAATACCCCCTTTTAACAATAAAAAAGGTTTGGCCACAGTTCCCTGGGGGAGGTTGGGGTAAGAATGATGTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCa                                                                                                                                                                                                                              <  2:724606/224‑1 (MQ=255)
      taataaCCCGCTTTATCAATCCAAAAGGTTTTGCCACAGTTACCTGGGTGGGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGAGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGttt                                                                                                                                                                                                                       <  2:262471/229‑1 (MQ=255)
        ataaCCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGGGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGCCGGCTCACTGtttt                                                                                                                                                                                                                      >  1:582883/1‑228 (MQ=255)
           aCCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTTcc                                                                                                                                                                                                                   >  1:597726/1‑228 (MQ=255)
               gCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGGCAAAGTCACGTCGGCCTCCTGttttttccgaaa                                                                                                                                                                                                              >  2:430974/1‑221 (MQ=255)
                cTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCCCTTTTCCCCAGCTTTTGGCAAAGCCGGCGTCAAGGTCACGGCGGCTACCCGTTTTTCTCTAAtt                                                                                                                                                                                                             >  1:298918/1‑229 (MQ=255)
                cTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGAGGAGCGTTGATGGCGTCAGCTTTGATTGAGATCTGAATTTCAGTCAAATACGCATGGTTACCCTGGCACTTTAGCTTAACTGCTGTAAAGTTCGCTTTTCCCCAGCTTTTGGCAAAGTCTGCGCCAAAGTCACGGCGGTTCAATACTTTTCCGGAAtt                                                                                                                                                                                                             >  2:237887/1‑229 (MQ=255)
                  ttAACAATCAACAAGGTTTTGCCAAAGTTACCTGGGTGGGGTTGTGGCAAGAATGAGGTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTAt                                                                                                                                                                                                           <  2:597725/229‑1 (MQ=255)
                      caatccCAAAGGTTTTGCCCCAGTTACCGGGGGGAGGTTGGGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGc                                                                                                                                                                                                       <  2:118147/229‑1 (MQ=255)
                        atcaCAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTTATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTGGGCAAAGGCGGGGTCACAGTCACGGCGGCTCCTTTTTTTTCCGCAATTTTGCGCaa                                                                                                                                                                                                     >  1:291058/1‑229 (MQ=255)
                             aaaGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCCGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGTTCCCCGTCTTTCTGTAATATTCCGCCAGGaaa                                                                                                                                                                                                >  2:170051/1‑228 (MQ=255)
                                 gTTTTGCCACAGTTCCCTGGGTGAGGGTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTc                                                                                                                                                                                            <  1:468183/229‑1 (MQ=255)
                                   tttGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGGTCACCGTTTTACCGCAATTATTCGCAAGACATTTTcc                                                                                                                                                                                           >  1:124156/1‑228 (MQ=255)
                                    ttGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCCCCGTTTTTCCGTAATTTTCCGCAAGGAATTTTcc                                                                                                                                                                                           >  1:275760/1‑227 (MQ=255)
                                    ttGCCACAGTTACCTGGGTGAGGATGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGGTCACTGGTTTTCCGTCATTATCCGCAAGGAATTTTcc                                                                                                                                                                                           >  1:712718/1‑227 (MQ=255)
                                     tGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAATGGCGGCGTCAAAGTCACGGCGGCTCCCCTGTTTTACGGCATTTATCGGCAGGGATTTTACGGc                                                                                                                                                                                        >  2:670224/1‑229 (MQ=255)
                                      gCCAGAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGAGCGAGATCTGAATTGCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCt                                                                                                                                                                                       <  1:170051/229‑1 (MQ=255)
                                       ccACAGCTACCTCGGCGTGGATGTGGCCAGCATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATAGAGATCTGAATTTCAGTCAAATACGCAGGGGTACCCTGGCACGTTAGTTTAACTCCTTTCACGTTCCTTTTTCCCCAGCTTTTGGCAGAGGCGGCGTCCAAGCCAGGGCGGCTTACCGTTTTACCGTAACAATCCGCAAGGGACTTTcgcgccg                                                                                                                                                                                      >  2:12710/1‑223 (MQ=255)
                                         acaGTTACCTGGGTGAGGTTGTGGGAAGAATGAGGTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCg                                                                                                                                                                                    <  1:70578/229‑1 (MQ=255)
                                           aGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCCCTGTTTTTCCGGAATTATCCCCAAGGAATTTTCCGGCCTCGct                                                                                                                                                                                  >  1:522778/1‑229 (MQ=255)
                                                        aGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCt                                                                                                                                                                     <  1:231337/229‑1 (MQ=255)
                                                             gtgGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGAtt                                                                                                                                                                <  1:313038/229‑1 (MQ=255)
                                                                   agagTGAGTTTTCAAAAAGGGGAGCGGTGAAGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATCCGCAGGGTTACCCTGGCAGGTTAGCCAAACTGCTTTCACGTTCTCTGTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCg                                                                                                                                                          <  2:360185/229‑1 (MQ=255)
                                                                   aGAATGAGTTTGGAAGAAGCGGAGCGTTGAGGGCGTCAGCTTTGATCGAGATCTGAATTTCAGGCAAATTCGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCg                                                                                                                                                          <  1:248198/229‑1 (MQ=255)
                                                                    gAATGGGTTTGCAGAAAGGGGAGGGTTGATGGCGTCGGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCg                                                                                                                                                          <  2:47346/228‑1 (MQ=255)
                                                                    gAATGAGTTTGCAGCAAGCGGAGCGTTGGTGGGGTCAACTGTGATCGAAACCTGAATTTCAGTCAAATCCGCAGGGTTGCCCGTGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCCTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCg                                                                                                                                                          <  1:237887/228‑1 (MQ=255)
                                                                     aaTGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTa                                                                                                                                                        >  1:272359/1‑229 (MQ=255)
                                                                      aTGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTa                                                                                                                                                        <  1:430974/228‑1 (MQ=255)
                                                                      aTGAGTTTGCAAAAAGTGGAGCGTTGTGGGCGTCAGCTTTGATCGAGATCTTAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTGTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTa                                                                                                                                                        <  2:140294/228‑1 (MQ=255)
                                                                        ggggttGAAGAAAGGGGAGGGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATTCGCAGGGTTCCCCTGGCATGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCa                                                                                                                                                     <  2:582882/225‑1 (MQ=255)
                                                                        gggTTTGCAGAAAGCGGAGCGTTGATGGGGTCAGCTTTGACCGAGATCTGAATTTCAGTCAAATCCGCAGGGTTCACCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCa                                                                                                                                                     <  1:284608/227‑1 (MQ=255)
                                                                         aGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGCAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTTAGTGGTAccct                                                                                                                                                    >  1:676037/1‑227 (MQ=255)
                                                                              tcaaaaaGCGGGGGGTTGGTGGCGTCAGCTTTGAACGGTATTTGAATTTCCGTCAAATACGACGGGTTACCCTGGCCAGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCAGCGTAc                                                                                                                                               <  2:298918/225‑1 (MQ=255)
                                                                                   aaGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTCTACCGTCATTATCCGCCAGGAATTTTCCGGCTTCCCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGGAcccgag                                                                                                                                          >  1:216379/1‑225 (MQ=255)
                                                                                      cGGAGCGTTGGTGGCGTCAGCTTTTGTCGTGATTTGAATTTCAGTCAAATACTCAGGGTTACCCTGGAACGTTAGCGTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTAt                                                                                                                                       <  2:243500/229‑1 (MQ=255)
                                                                                      cGGAGCGTTGATGGCTGCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTAt                                                                                                                                       <  1:428507/229‑1 (MQ=255)
                                                                                             ttGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCGAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTc                                                                                                                                  <  1:168924/227‑1 (MQ=255)
                                                                                                 tggggtccacTTTTACCGAGATCTGAATTTCAGTCAAATCCGCAGGGTTCCCCTGGGACGTTAGGTTAACTGGTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGAt                                                                                                                             <  2:272359/219‑1 (MQ=255)
                                                                                                 tGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGAt                                                                                                                             <  1:77930/228‑1 (MQ=255)
                                                                                                          cTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCCCCCTGTTTTTCCCTAATTATCCGCCCGGGATTTTCCGCCTTCCCTCTCTTTAATTTCTTTATTCATCCCTACCATCGTACCGACGAATCCGTCCCGATCAAAACCaaa                                                                                                                   >  2:335217/1‑229 (MQ=255)
                                                                                                            ttGAACGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTGAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGc                                                                                                                 <  1:743091/229‑1 (MQ=255)
                                                                                                                 tgggaTCTGAATTTCCGTCAATTTCGCACGGTTACCCTGGCACGTTAGCTTAACTCTTTTCCCGTTCTCTTTTCCCCATCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAgggt                                                                                                            <  2:87272/226‑3 (MQ=255)
                                                                                                                 cGATATCTGAAGTTGAGTCAATTTCGTAGGGTTCCCCTGTGCCGGTAGCTTAACTTCTTTCACGTTCTTGTTTTCCCAGCTTTTGGCAAAGGCTGCGTCAAAGTCACGGCGGCTCCCTGTTTTACCGTAATTAGCCGCATGGAATTTTCCGGCTTCGCTCACTTTAATGTCTTGATTCAGGCGTTCCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCa                                                                                                            <  2:64907/229‑1 (MQ=255)
                                                                                                                 cGAGATCTGAATTTCAGTCAAATACGCAGGGTTAACCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCa                                                                                                            <  1:256943/229‑1 (MQ=255)
                                                                                                                    gaTCTGAATTTCAGTCAAATACGCAGGGTTACCCTGTCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATTCGCCAAGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCATGAGTGACATCGTAACACCGCTACAGTCCGGATCAAAACCCAAACCGGCAcca                                                                                                         >  1:444885/1‑228 (MQ=255)
                                                                                                                    gaTCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCGCCAGCTTTTGGCAAAGGCGGCGGCAAAGTCACGGCGGTTGCCTTTTTTACCGTAATTATCCGCCAGGGATTTTCCGCGTACCGCCTCTTTAATTTTTTGTTTCACTCGTACCACCGTACCGAAGGATTCCTCCGGTTCCAAACCCAAACAgagccag                                                                                                         >  2:537848/1‑223 (MQ=255)
                                                                                                                          aaTTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCCAGGCATTTTCCGGCTTCGCACTCTTTAATTTTTTGTTTCAGGCGTAACATCGTACCCAAGTATGCGTCCGGTTCAAAACCAAACACGGCACAGGGTTTg                                                                                                   >  2:172822/1‑229 (MQ=255)
                                                                                                                                  tCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCCCTGTTTTACCGTAATTATCCGCCAGGAATTTTTCGGCTTCGCTCTCTTTAATTTCTTGATTCCGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAACCAGGCACCCGGTTTGGCATCttt                                                                                           >  1:45318/1‑228 (MQ=255)
                                                                                                                                     aagtcgcgGGTTTTCCTTGGCAGGTTGGGTTACCTGCTTTCACGTGCTCTTTTCCCCATCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGt                                                                                         <  2:687150/224‑1 (MQ=255)
                                                                                                                                      aTACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTa                                                                                        <  1:172822/228‑1 (MQ=255)
                                                                                                                                            aGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTTCCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTTATTTAGCCGTACCATCGTACCCAAGAATTCTCCCGGATTCAAACCAAAACAGGCCCCCGGTTTGGCCTTTTCTTAGGCTTcc                                                                                 >  2:6419/1‑229 (MQ=255)
                                                                                                                                              ggTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGGTTCGCTCTCTTTAATTTCTTGATTCAAGTGTACCATCGTACCGAAGTATGCGTCCGGGTCAAAACCAAAGAAGTGACCGTGTTTGGCATATTCGTAGCGTTCca                                                                                >  1:114727/1‑228 (MQ=255)
                                                                                                                                                              tAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAAGCGTACCATCGTACCGAAGTATGCGTCCGGGTCAAAACCAAAGCAGGCACCGTGTTTGGCAATTTCGTAGCGTTCCAGGCAGGAAGGTCCGCCa                                                               >  1:258429/1‑229 (MQ=255)
                                                                                                                                                                 cTTAACTGTTTTCCCGTTCTCTTTGCCCCCGTTTTTGGACAAGGCGCCGTCAAAGTCACGGCGGCTCACGGTTTTACCGTAATTATCCTCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCt                                                            <  2:114727/229‑1 (MQ=255)
                                                                                                                                                                 cTTAACTGCTTTCACGTTCCCTTTTCCCCAGCTTTTGGCAGAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTGCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCt                                                            <  1:12710/229‑1 (MQ=255)
                                                                                                                                                                       tGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGGTTTTCCGTAATTATCCGCCAGGAATTTTCCGGGTTCGCTCTCTTTAATTTCTTGATTCATGCGTACCATCGTACCGAAGTATGTGTCCGGATCAACACCAAAGCACGCACCGTGCTTGGCATATTCGTAGCGTTACAAGCAGGAACGTCCGCCAGCTCCTg                                                        >  1:447456/1‑227 (MQ=255)
                                                                                                                                                                               cGTTCTCTTTTCCCCAGCTTTTGGGAAAGGCGGCGTCAAAGCCCCGGCGGCTCACTGTTTTACCGTAATTATCCGCTAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCCTCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTc                                              <  2:447456/229‑1 (MQ=255)
                                                                                                                                                                                   tttttttCCCCAGCTTTTGGCAATGGCCGCGTCAAAGTCACGGCGGTCCATGTTGTTCCCGTAATTAGCCGCAAGGAATTTTGCGGCTTCGCTCTCTGTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGGGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACtt                                          <  2:13659/226‑1 (MQ=255)
                                                                                                                                                                                   ctctTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCCAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGCCGTACCATCCTACCGAAGTATGCGTCCGGGTCACAACCCAAGCAGGCACCGTGTTTGGGATATTCGTATCGTTCCAGGCAGGAACGCCCCCCACCCCCCGCCTTGACTTCAcct                                          >  2:515531/1‑227 (MQ=255)
                                                                                                                                                                                                       cAAAGGCGGCGTCAAAGTCACGGCGGCTTCATGTTTTACCGTAATTAGCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAg                      <  1:77886/229‑1 (MQ=255)
                                                                                                                                                                                                                         gccggcggCTCACTGGTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTGCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCg     <  1:503489/226‑1 (MQ=255)
                                                                                                                                                                                                                         ctcggcggCTTACTTTTTTACCGTAATTACCCGGAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCg     <  1:526046/226‑1 (MQ=255)
                                                                                                                                                                                                                           cggcggCTCACTTTTTTACAGTAATTAGCCGCAAGGAATTTTCCGGTTTCGCTCTCTTTAATTGCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTGTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCg  <  2:275760/229‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                            
GATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCTGGCACGTTAGCTTAACTGCTTTCACGTTCTCTTTTCCCCAGCTTTTGGCAAAGGCGGCGTCAAAGTCACGGCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCG  >  NC_000913/644192‑644639

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