Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 29,839 C→T 5.6% V63V (GTC→GTT carA → carbamoyl phosphate synthetase small subunit, glutamine amidotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091329,8390CT 5.6% 239.4 / 3.3 71V63V (GTC→GTTcarAcarbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base T (2/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

AGTTTCACGGTCGGGCCATAGGGGCAACAGGTTCGGCGGTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGC  >  NC_000913/29691‑29982
                                                                                                                                                    |                                                                                                                                               
aGTTTCACGGTCGGGCCATAGGGGCAACAGGTTCGGCGGTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTc                                                                                                                                                 >  4:259178/1‑149 (MQ=255)
    tCACGGTCGGGCCATAGGGGCAACAGGTTCGGCGGTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCa                                                                                                                                             >  2:234612/1‑149 (MQ=255)
            gggCCATAGGGGCAACAGGTTCGGCGGTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAAt                                                                                                                                        <  3:248511/146‑1 (MQ=255)
            gggCCATAGGGGCAACAGGTTCGGCGGTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAAt                                                                                                                                        >  4:248511/1‑146 (MQ=255)
            gggCCATAGGGGCAACAGGTTCGGCGGTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGAc                                                                                                                                     <  2:15652/149‑1 (MQ=255)
                 aTAGGGGCAACAGGTTCGGCGGTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGa                                                                                                                                <  2:455658/149‑1 (MQ=255)
                     gggCAACAGGTTCGGCGGTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATgaa                                                                                                                            >  3:339899/1‑149 (MQ=255)
                     gggCAACAGGTTCGGCGGTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATgaa                                                                                                                            <  1:234612/149‑1 (MQ=255)
                     gggCAACAGGTTCGGCGGTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATgaa                                                                                                                            <  2:147687/149‑1 (MQ=255)
                         aaCAGGTTCGGCGGTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAAt                                                                                                                        <  4:187740/149‑1 (MQ=255)
                         aaCAGGTTCGGCGGTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAAt                                                                                                                        >  1:98109/1‑149 (MQ=255)
                                  ggcggTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAgg                                                                                                               <  2:82076/149‑1 (MQ=255)
                                      gTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACa                                                                                                           <  3:303753/149‑1 (MQ=255)
                                               gTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAgg                                                                                                  >  3:18127/1‑149 (MQ=255)
                                                   ttttCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGt                                                                                                 <  2:37942/146‑1 (MQ=255)
                                                   ttttCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGt                                                                                                 >  1:37942/1‑146 (MQ=255)
                                                       cAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGa                                                                                          <  4:339899/149‑1 (MQ=255)
                                                        aaTACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGAt                                                                                         >  4:54131/1‑149 (MQ=255)
                                                             ttCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTa                                                                                                             <  1:443740/124‑1 (MQ=255)
                                                             ttCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTa                                                                                                             >  2:443740/1‑124 (MQ=255)
                                                             ttCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTcgcg                                                                                    >  3:141760/1‑149 (MQ=255)
                                                             ttCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCCTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTcgcg                                                                                    >  2:255833/1‑149 (MQ=255)
                                                                    aCCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGcc                                                                             >  1:422709/1‑149 (MQ=255)
                                                                        gTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACAt                                                                                                          <  3:180089/116‑1 (MQ=255)
                                                                        gTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACAt                                                                                                          >  4:180089/1‑116 (MQ=255)
                                                                        gTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTg                                                                         <  4:141760/149‑1 (MQ=255)
                                                                         ttATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGa                                                                        <  2:260246/149‑1 (MQ=255)
                                                                          tATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGAt                                                                       >  3:8789/1‑149 (MQ=255)
                                                                           aTCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGAtt                                                                      >  4:407857/1‑149 (MQ=255)
                                                                              aaGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGt                                                                                                              >  1:84524/1‑106 (MQ=255)
                                                                              aaGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGt                                                                                                              <  2:84524/106‑1 (MQ=255)
                                                                               aGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCa                                                                  <  3:335939/149‑1 (MQ=255)
                                                                               aGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCa                                                                  <  2:422709/149‑1 (MQ=255)
                                                                               aGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCa                                                                  <  3:407857/149‑1 (MQ=255)
                                                                               aGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCa                                                                  >  2:128281/1‑149 (MQ=255)
                                                                                gAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTgg                                                                                             <  2:142505/121‑1 (MQ=255)
                                                                                gAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTgg                                                                                             >  1:142505/1‑121 (MQ=255)
                                                                                   aTCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCaa                                                              <  3:54131/149‑1 (MQ=255)
                                                                                       tCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGc                                                                    <  4:306603/139‑1 (MQ=255)
                                                                                       tCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGc                                                                    >  3:306603/1‑139 (MQ=255)
                                                                                              tccttcctATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATa                                                   <  2:98109/149‑1 (MQ=255)
                                                                                                          tcGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAAt                                                                                                                        >  1:272980/1‑68 (MQ=255)
                                                                                                          tcGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAAt                                                                                                                        <  2:272980/68‑1 (MQ=255)
                                                                                                            gTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTTGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACctctct                                     >  3:444611/1‑149 (MQ=255)
                                                                                                            gTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTTGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACctctct                                     >  1:69848/1‑149 (MQ=255)
                                                                                                            gTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGc                                                                    >  2:171342/1‑118 (MQ=255)
                                                                                                            gTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGc                                                                    <  1:171342/118‑1 (MQ=255)
                                                                                                             tCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCtctt                                    >  2:2513/1‑149 (MQ=255)
                                                                                                                 aTCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTc                                                                                       >  2:338467/1‑94 (MQ=255)
                                                                                                                 aTCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTc                                                                                       <  1:338467/94‑1 (MQ=255)
                                                                                                                 aTCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTa                                <  4:18127/149‑1 (MQ=255)
                                                                                                                  tCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTAc                               <  1:91255/149‑1 (MQ=255)
                                                                                                                    gTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAAtc                                                                                                                       <  3:143507/59‑1 (MQ=255)
                                                                                                                    gTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAAtc                                                                                                                       >  4:143507/1‑59 (MQ=255)
                                                                                                                        ctctTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGaaa                         >  3:137378/1‑149 (MQ=255)
                                                                                                                        ctctTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGaaa                         <  2:56571/149‑1 (MQ=255)
                                                                                                                         tctTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGa                                                                        <  1:159045/101‑1 (MQ=255)
                                                                                                                         tctTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGa                                                                        >  2:159045/1‑101 (MQ=255)
                                                                                                                         tctTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAAc                        >  1:379110/1‑149 (MQ=255)
                                                                                                                                   ccccATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACAt              >  3:468432/1‑149 (MQ=255)
                                                                                                                                    cccATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATc             >  3:161946/1‑149 (MQ=255)
                                                                                                                                     ccATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGt                                                                                                              >  4:140229/1‑51 (MQ=255)
                                                                                                                                     ccATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGt                                                                                                              <  3:140229/51‑1 (MQ=255)
                                                                                                                                        tatTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTa                                <  3:434991/126‑1 (MQ=255)
                                                                                                                                        tatTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTa                                >  4:434991/1‑126 (MQ=255)
                                                                                                                                            ggCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTg                                                                                              >  4:156138/1‑60 (MQ=255)
                                                                                                                                            ggCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTg                                                                                              <  3:156138/60‑1 (MQ=255)
                                                                                                                                               aaTGTTGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGc  <  2:69848/149‑1 (MQ=255)
                                                                                                                                               aaTGTTGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGc  <  4:444611/149‑1 (MQ=255)
                                                                                                                                               aaTGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGc  <  4:8789/149‑1 (MQ=255)
                                                                                                                                                 tGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGa                                            <  2:258977/105‑1 (MQ=255)
                                                                                                                                                 tGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGa                                            >  1:258977/1‑105 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AGTTTCACGGTCGGGCCATAGGGGCAACAGGTTCGGCGGTTGGGGAAGTCGTTTTCAATACTTCAATGACCGGTTATCAAGAAATCCTCACTGATCCTTCCTATTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGC  >  NC_000913/29691‑29982

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