Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,694,475 C→A 6.6% E357* (GAG→TAG)  purL ← phosphoribosylformyl‑glycineamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,694,4750CA 6.6% 215.7 / 5.2 61E357* (GAG→TAG) purLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base A (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

GCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGC  >  NC_000913/2694331‑2694622
                                                                                                                                                |                                                                                                                                                   
gcAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGc                                                                                                                                                 <  4:412422/149‑1 (MQ=255)
    ctctTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTAc                                                                                                                                             >  3:368251/1‑149 (MQ=255)
       ttCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTAcc                                                                                                                                            >  1:149362/1‑147 (MQ=255)
       ttCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTAcc                                                                                                                                            <  2:790813/147‑1 (MQ=255)
       ttCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTAcc                                                                                                                                            >  1:790813/1‑147 (MQ=255)
       ttCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTAcc                                                                                                                                            <  2:149362/147‑1 (MQ=255)
       ttCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGa                                                                                                                                          >  2:605597/1‑149 (MQ=255)
         cGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGaaa                                                                                                                                        <  1:424714/149‑1 (MQ=255)
                 tgGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCttcttc                                                                                                                                <  1:605597/149‑1 (MQ=255)
                   gCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTccc                                                                                                                              <  1:231644/149‑1 (MQ=255)
                       ttCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACgg                                                                                                                          <  2:677807/149‑1 (MQ=255)
                        tCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGc                                                                                                                         <  3:668435/149‑1 (MQ=255)
                         cACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCt                                                                                                                        <  2:66940/149‑1 (MQ=255)
                           cTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGt                                                                                                                      >  3:58116/1‑149 (MQ=255)
                            tttttCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGtt                                                                                                                     <  3:662201/149‑1 (MQ=255)
                            tttttCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGtt                                                                                                                     <  4:463264/149‑1 (MQ=255)
                            tttttCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGtt                                                                                                                     >  4:555480/1‑149 (MQ=255)
                                tcttcATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGaa                                                                                                                 <  2:463195/149‑1 (MQ=255)
                                  ttcATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGc                                                                                                               >  2:307739/1‑149 (MQ=255)
                                  ttcATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGc                                                                                                               >  4:645247/1‑149 (MQ=255)
                                     aTAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAg                                                                                                            <  3:555480/149‑1 (MQ=255)
                                     aTAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAg                                                                                                            <  4:76056/149‑1 (MQ=255)
                                     aTAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAg                                                                                                            >  1:630952/1‑149 (MQ=255)
                                     aTAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAg                                                                                                            >  2:197836/1‑149 (MQ=255)
                                         gTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAAt                                                                                                        >  4:686212/1‑149 (MQ=255)
                                           aCGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTc                                                                                                      >  3:61975/1‑149 (MQ=255)
                                             ggAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTAAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGc                                                                                                    >  4:36840/1‑149 (MQ=255)
                                             ggAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTAAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGc                                                                                                    >  4:112085/1‑149 (MQ=255)
                                                  tAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGtt                                                                                               >  2:532077/1‑149 (MQ=255)
                                                  tAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGtt                                                                                               <  3:645247/149‑1 (MQ=255)
                                                     ccGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGtt                                                                                               >  3:606929/1‑146 (MQ=255)
                                                     ccGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGtt                                                                                               <  4:606929/146‑1 (MQ=255)
                                                        ttCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAAc                                                                                         >  3:605864/1‑149 (MQ=255)
                                                          cAGTGCCGGACGACCGAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCATGTTGGAAACgg                                                                                       <  2:47289/149‑1 (MQ=255)
                                                          cAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACgg                                                                                       <  2:262516/149‑1 (MQ=255)
                                                               ccGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGaaa                                                                                  <  1:544327/149‑1 (MQ=255)
                                                                  gacgacCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAaacca                                                                               >  3:350543/1‑149 (MQ=255)
                                                                   acgacCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTAAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAaaccaa                                                                              <  3:36840/149‑1 (MQ=255)
                                                                   acgacCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTAAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAaaccaa                                                                              <  3:112085/149‑1 (MQ=255)
                                                                    cgacCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAaaccaac                                                                             <  1:532077/149‑1 (MQ=255)
                                                                        cAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGa                                                                         >  4:481390/1‑149 (MQ=255)
                                                                                ttgttAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCTCACGGCTGTTCGAAGcc                                                                                                              >  3:483828/1‑104 (MQ=255)
                                                                                ttgttAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCTCACGGCTGTTCGAAGcc                                                                                                              <  4:483828/104‑1 (MQ=255)
                                                                                 tgttAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTc                                                                <  2:630952/149‑1 (MQ=255)
                                                                                   ttAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCgg                                                              >  4:290633/1‑149 (MQ=255)
                                                                                     aaaCGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCt                                                            >  4:677367/1‑149 (MQ=255)
                                                                                     aaaCGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCt                                                            >  2:603176/1‑149 (MQ=255)
                                                                                        cgccgcgccgCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTg                                                         <  4:58116/149‑1 (MQ=255)
                                                                                         gccgcgccgcCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTgc                                                        >  2:180643/1‑149 (MQ=255)
                                                                                         gccgcgccgcCCAGCGGGCCTCCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTgc                                                        >  4:11420/1‑149 (MQ=255)
                                                                                                cgccCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTgcgccgcgc                                                 >  2:439302/1‑149 (MQ=255)
                                                                                                      gCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTg                                           <  3:290633/149‑1 (MQ=255)
                                                                                                           ccTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCAcc                                      <  4:350543/149‑1 (MQ=255)
                                                                                                                       atgTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGAt                          >  4:136786/1‑149 (MQ=255)
                                                                                                                           ccAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTg                                           <  4:290435/128‑1 (MQ=255)
                                                                                                                           ccAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTg                                           >  3:290435/1‑128 (MQ=255)
                                                                                                                           ccAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCa                      <  3:136786/149‑1 (MQ=255)
                                                                                                                            cAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCAc                     >  4:223664/1‑149 (MQ=255)
                                                                                                                                 cgGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCAccgccg                <  3:504056/149‑1 (MQ=255)
                                                                                                                                        aaTGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCgg         >  2:164483/1‑149 (MQ=255)
                                                                                                                                               tCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTcgccgc  <  1:164483/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGC  >  NC_000913/2694331‑2694622

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