Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,775,276 T→C 6.1% D197D (GAT→GAC aroD → 3‑dehydroquinate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,775,2760TC 6.1% 213.6 / 4.0 65D197D (GAT→GACaroD3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (26/35);  new base C (2/2);  total (28/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 = 4.46e-01

GCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAA  >  NC_000913/1775129‑1775424
                                                                                                                                                   |                                                                                                                                                    
gCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATc                                                                                                                                                     <  3:195749/149‑1 (MQ=255)
gCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATc                                                                                                                                                     <  3:226518/149‑1 (MQ=255)
gCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTCACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATc                                                                                                                                                     >  1:203044/1‑149 (MQ=255)
       gCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAAt                                                                                                                                              <  3:516140/149‑1 (MQ=255)
          gaagaaATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTAt                                                                                                                                           >  1:658328/1‑149 (MQ=255)
               aaTCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGa                                                                                                                                      <  2:153672/149‑1 (MQ=255)
               aaTCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGa                                                                                                                                      >  2:504410/1‑149 (MQ=255)
                aTCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGAt                                                                                                                                     <  3:630184/149‑1 (MQ=255)
                  cATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGt                                                                                                                                   >  2:239158/1‑149 (MQ=255)
                  cATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGACGCCCCACATTAACAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCACTATGCCGATCGTCCAATTATCACGATGt                                                                                                                                   <  2:495060/149‑1 (MQ=255)
                   aTTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTc                                                                                                                                  >  3:32546/1‑149 (MQ=255)
                      gCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGTAGTATGCCGATCGTCCAATTATCACGATGTCGAt                                                                                                                               >  1:199584/1‑149 (MQ=255)
                      gCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGAt                                                                                                                               <  1:164578/149‑1 (MQ=255)
                      gCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGACGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGAt                                                                                                                               <  4:136253/149‑1 (MQ=255)
                        ccGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATgg                                                                                                                             >  1:509653/1‑149 (MQ=255)
                                  aaaaTGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTgg                                                                                                                   >  1:236252/1‑149 (MQ=255)
                                     aTGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGACCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGt                                                                                                                >  1:255763/1‑149 (MQ=255)
                                     aTGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGACCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGt                                                                                                                >  1:207683/1‑149 (MQ=255)
                                         aaTCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTaa                                                                                                              >  2:256939/1‑147 (MQ=255)
                                         aaTCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTaa                                                                                                              <  1:256939/147‑1 (MQ=255)
                                         aaTCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAAtt                                                                                                            <  2:236252/149‑1 (MQ=255)
                                             cTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGa                                                                                                                                <  2:471779/125‑1 (MQ=255)
                                             cTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGa                                                                                                                                >  1:471779/1‑125 (MQ=255)
                                              ttCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCg                                                                                                       >  4:44761/1‑149 (MQ=255)
                                                 gACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTct                                                                                                    >  1:555009/1‑149 (MQ=255)
                                                 gACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTct                                                                                                    <  4:32546/149‑1 (MQ=255)
                                                    gCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctggc                                                                                                 <  1:173969/149‑1 (MQ=255)
                                                     ccGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctggct                                                                                                <  2:533162/149‑1 (MQ=255)
                                                     ccGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctggct                                                                                                <  1:391868/149‑1 (MQ=255)
                                                       gATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctggctgg                                                                                              <  1:504410/149‑1 (MQ=255)
                                                           ttCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGAAGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGaa                                                                                          <  3:123892/149‑1 (MQ=255)
                                                             ccTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGt                                                                                        <  2:509653/149‑1 (MQ=255)
                                                              cTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTa                                                                                       >  2:130662/1‑149 (MQ=255)
                                                               tAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTAt                                                                                      >  2:253927/1‑149 (MQ=255)
                                                               tAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTAt                                                                                      >  4:278145/1‑149 (MQ=255)
                                                                aaGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTAtt                                                                                     <  2:555009/149‑1 (MQ=255)
                                                                aaGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGACCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTAtt                                                                                     <  2:255763/149‑1 (MQ=255)
                                                                aaGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGACCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTAtt                                                                                     <  2:207683/149‑1 (MQ=255)
                                                                 aGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTAttt                                                                                    <  2:643431/149‑1 (MQ=255)
                                                                   aTTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTgg                                                                                  <  2:658328/149‑1 (MQ=255)
                                                                    ttGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGTAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGc                                                                                 <  2:199584/149‑1 (MQ=255)
                                                                              gCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCaa                                                                       >  1:218232/1‑149 (MQ=255)
                                                                                cGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACt                                                                     >  1:85369/1‑149 (MQ=255)
                                                                                 gCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACtt                                                                    >  2:39593/1‑149 (MQ=255)
                                                                                                    cTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGc                                                 <  1:620426/149‑1 (MQ=255)
                                                                                                              ttgcCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTAttt                                                                                    <  4:653069/104‑1 (MQ=255)
                                                                                                              ttgcCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTAttt                                                                                    >  3:653069/1‑104 (MQ=255)
                                                                                                                       cccTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAAtc                              <  3:44761/149‑1 (MQ=255)
                                                                                                                         cTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAAtctc                            >  3:266510/1‑149 (MQ=255)
                                                                                                                          tGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTgc                                           <  4:233885/133‑1 (MQ=255)
                                                                                                                          tGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTgc                                           >  3:233885/1‑133 (MQ=255)
                                                                                                                          tGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCg                           <  1:582904/149‑1 (MQ=255)
                                                                                                                          tGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCg                           >  2:87486/1‑149 (MQ=255)
                                                                                                                           ggAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCgg                          >  4:568868/1‑149 (MQ=255)
                                                                                                                            gagaTGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCa                                       >  3:640603/1‑135 (MQ=255)
                                                                                                                            gagaTGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCa                                       <  4:640603/135‑1 (MQ=255)
                                                                                                                            gagaTGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGt                         <  1:253927/149‑1 (MQ=255)
                                                                                                                                   aGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGAt                  >  1:193761/1‑149 (MQ=255)
                                                                                                                                    ggAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGAtt                 <  2:203044/149‑1 (MQ=255)
                                                                                                                                     gAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGAttt                <  1:238864/149‑1 (MQ=255)
                                                                                                                                     gAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGAttt                <  4:435366/149‑1 (MQ=255)
                                                                                                                                     gAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGAttt                <  1:518791/149‑1 (MQ=255)
                                                                                                                                     gAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGAttt                >  4:81153/1‑149 (MQ=255)
                                                                                                                                       gCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTgc              <  3:278145/149‑1 (MQ=255)
                                                                                                                                            aTGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACg         <  1:239158/149‑1 (MQ=255)
                                                                                                                                                   tCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTaa  <  1:130662/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAA  >  NC_000913/1775129‑1775424

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