Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,321,427 A→T 7.5% L507Q (CTG→CAG)  trpE ← component I of anthranilate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,321,4270AT 7.5% 376.5 / 16.4 107L507Q (CTG→CAG) trpEcomponent I of anthranilate synthase
Reads supporting (aligned to +/- strand):  ref base A (51/48);  new base T (3/5);  total (54/53)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 107 total aligned reads displayed.

GGTTGCGGTAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCC  >  NC_000913/1321285‑1321569
                                                                                                                                              |                                                                                                                                              
ggTTGCGGTAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAg                                                                                                                                          >  3:369447/1‑149 (MQ=255)
ggTTGCGGTAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAg                                                                                                                                          >  1:599743/1‑149 (MQ=255)
  ttGCGGTAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAgcg                                                                                                                                        >  4:143806/1‑149 (MQ=255)
   tGCGGTAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGACAATAGCGCGCAGTACAgcgc                                                                                                                                       <  3:612701/149‑1 (MQ=255)
    gCGGTAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGcgcg                                                                                                                                      <  4:393910/149‑1 (MQ=255)
     cGGTAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCTGTACTGcgcgc                                                                                                                                     <  3:370469/149‑2 (MQ=255)
        tAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCTGTACTGcgcgctat                                                                                                                                  <  3:585356/149‑5 (MQ=255)
        tAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCt                                                                                                                                  <  2:599743/149‑1 (MQ=255)
          aaTCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCttt                                                                                                                                <  3:460098/149‑1 (MQ=255)
           aTCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTg                                                                                                                               <  2:515079/149‑1 (MQ=255)
             caccacGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGtt                                                                                                                             <  3:531254/149‑1 (MQ=255)
                cacGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACg                                                                                                                          >  4:514592/1‑149 (MQ=255)
                    ttATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGtt                                                                                                                      >  3:745064/1‑149 (MQ=255)
                     tATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGttt                                                                                                                     >  4:570834/1‑149 (MQ=255)
                         cccATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTtcgtc                                                                                                                 >  1:225034/1‑149 (MQ=255)
                         cccATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTtcgtc                                                                                                                 >  2:6625/1‑149 (MQ=255)
                          ccATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTtcgtcg                                                                                                                >  2:559969/1‑149 (MQ=255)
                          ccATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTtcgtcg                                                                                                                >  1:506488/1‑149 (MQ=255)
                            aTTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGcgcg                                                                                                                                      <  3:339746/125‑1 (MQ=255)
                            aTTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGcgcg                                                                                                                                      >  4:339746/1‑125 (MQ=255)
                               gctgcGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACgg                                                                                                                         >  1:343157/1‑135 (MQ=255)
                               gctgcGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACgg                                                                                                                         <  2:343157/135‑1 (MQ=255)
                               gctgcGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTc                                                                                                           <  3:33349/149‑1 (MQ=255)
                                 tgcGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCg                                                                                                         >  3:665670/1‑149 (MQ=255)
                                  gcgcAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCGTTGTTACGGGTTTCGTCGGCTTCCGa                                                                                                        <  3:496538/149‑1 (MQ=255)
                                      aaCTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGc                                                                                                    <  3:514592/149‑1 (MQ=255)
                                      aaCTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGc                                                                                                    >  1:809965/1‑149 (MQ=255)
                                      aaCTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGc                                                                                                    <  1:559969/149‑1 (MQ=255)
                                      aaCTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGc                                                                                                    <  2:225034/149‑1 (MQ=255)
                                       aCTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCg                                                                                                   <  1:330657/149‑1 (MQ=255)
                                         tGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGa                                                                                                 >  4:335942/1‑149 (MQ=255)
                                         tGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGa                                                                                                 >  2:281928/1‑149 (MQ=255)
                                               gCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGaa                                                                                           >  3:176488/1‑149 (MQ=255)
                                                 cAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATc                                                                                         >  4:491252/1‑149 (MQ=255)
                                                   ggTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCaa                                                                                       >  1:587876/1‑149 (MQ=255)
                                                   ggTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCaa                                                                                       <  1:538609/149‑1 (MQ=255)
                                                           gTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGa                                                                                    >  3:405885/1‑144 (MQ=255)
                                                           gTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGa                                                                                    <  4:405885/144‑1 (MQ=255)
                                                            tAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCaa                                                                                       >  3:362624/1‑140 (MQ=255)
                                                            tAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCaa                                                                                       <  4:362624/140‑1 (MQ=255)
                                                            tAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTacac                                                                              <  3:335942/149‑1 (MQ=255)
                                                            tAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTacac                                                                              >  4:281823/1‑149 (MQ=255)
                                                             aaaaGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCg                                                                                                   >  1:190243/1‑127 (MQ=255)
                                                             aaaaGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCg                                                                                                   <  2:190243/127‑1 (MQ=255)
                                                             aaaaGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACAcc                                                                             >  3:746149/1‑149 (MQ=255)
                                                                 gagTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGa                                                                                    >  2:206042/1‑138 (MQ=255)
                                                                 gagTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGa                                                                                    <  1:206042/138‑1 (MQ=255)
                                                                 gagTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCa                                                                         <  2:506488/149‑1 (MQ=255)
                                                                    tCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCAccc                                                                      >  4:680943/1‑149 (MQ=255)
                                                                    tCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCAccc                                                                      >  1:134766/1‑149 (MQ=255)
                                                                       atatTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGa                                                                                    <  3:555059/132‑1 (MQ=255)
                                                                       atatTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGa                                                                                    >  4:555059/1‑132 (MQ=255)
                                                                           tatCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCa                                                               >  3:615477/1‑149 (MQ=255)
                                                                             tCGAGCATCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACg                                                             >  2:440458/1‑149 (MQ=255)
                                                                             tCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGa                                                                                    <  2:375360/126‑1 (MQ=255)
                                                                             tCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGa                                                                                    >  1:375360/1‑126 (MQ=255)
                                                                                  cagcagAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGc                                                        <  4:665670/149‑1 (MQ=255)
                                                                                   agcagAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCg                                                       >  2:344587/1‑149 (MQ=255)
                                                                                    gcagAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGa                                                      <  3:491252/149‑1 (MQ=255)
                                                                                       gAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAg                                                                                      >  1:207174/1‑114 (MQ=255)
                                                                                       gAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAg                                                                                      <  4:724107/114‑1 (MQ=255)
                                                                                            tCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGttt                                              <  4:531432/149‑1 (MQ=255)
                                                                                            tCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGttt                                              >  4:111025/1‑149 (MQ=255)
                                                                                            tCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGttt                                              <  2:809965/149‑1 (MQ=255)
                                                                                             cAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATTAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGtttt                                             >  3:217876/1‑149 (MQ=255)
                                                                                             cAGCCATCAGAAAGGCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGtttt                                             <  4:745064/149‑1 (MQ=255)
                                                                                              aGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTc                                            >  1:566401/1‑149 (MQ=255)
                                                                                              aGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTc                                            <  2:134766/149‑1 (MQ=255)
                                                                                              aGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTc                                            >  1:771452/1‑149 (MQ=255)
                                                                                                ccATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTcca                                          <  1:281928/149‑1 (MQ=255)
                                                                                                   tCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTccacca                                       >  3:27100/1‑149 (MQ=255)
                                                                                                   tCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTccacca                                       <  4:176488/149‑1 (MQ=255)
                                                                                                    cAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGttt                                              >  4:426987/1‑141 (MQ=255)
                                                                                                    cAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGttt                                              <  3:426987/141‑1 (MQ=255)
                                                                                                         aGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCg                                 <  4:256976/149‑1 (MQ=255)
                                                                                                          gTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGa                                >  2:686626/1‑149 (MQ=255)
                                                                                                           tctcCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAg                               <  1:445085/149‑1 (MQ=255)
                                                                                                           tctcCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAg                               >  1:269022/1‑149 (MQ=255)
                                                                                                             tcCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTcc                                           >  4:335402/1‑135 (MQ=255)
                                                                                                             tcCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTcc                                           <  3:335402/135‑1 (MQ=255)
                                                                                                               cTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGa                                                                                    <  3:314713/92‑1 (MQ=255)
                                                                                                               cTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGa                                                                                    >  4:314713/1‑92 (MQ=255)
                                                                                                                tgtgCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTcc                                                                                                          >  4:246803/1‑69 (MQ=255)
                                                                                                                tgtgCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTcc                                                                                                          <  3:246803/69‑1 (MQ=255)
                                                                                                                tgtgCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGAt                          >  3:713437/1‑149 (MQ=255)
                                                                                                                     atgatgCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACaa                     <  3:281823/149‑1 (MQ=255)
                                                                                                                               gTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGgtggc           <  1:686626/149‑3 (MQ=255)
                                                                                                                               gTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTc           <  3:143806/149‑1 (MQ=255)
                                                                                                                                tGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCg          >  1:530650/1‑149 (MQ=255)
                                                                                                                                 ggCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCga         <  1:344587/149‑1 (MQ=255)
                                                                                                                                    aaTAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGAt      <  4:412941/149‑1 (MQ=255)
                                                                                                                                    aaTAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGAt      >  1:743883/1‑149 (MQ=255)
                                                                                                                                     aTAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATc     >  1:211282/1‑149 (MQ=255)
                                                                                                                                     aTAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATc     <  3:570834/149‑1 (MQ=255)
                                                                                                                                      gcccgcgcTGTACTGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACCGAATCAAGGa                                                                                    <  3:49460/66‑1 (MQ=38)
                                                                                                                                      gcccgcgcTGTACTGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGttt                                              >  2:517358/4‑107 (MQ=255)
                                                                                                                                      gcccgcgcTGTACTGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGttt                                              <  1:517358/104‑1 (MQ=255)
                                                                                                                                      gcccgcgcTGTACTGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGAt      >  4:116517/4‑147 (MQ=255)
                                                                                                                                      gcccgcgcTGTACTGCGCGGGCTATGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGa                                                                                    >  4:49460/4‑69 (MQ=38)
                                                                                                                                       aGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGc   <  1:6625/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
GGTTGCGGTAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCC  >  NC_000913/1321285‑1321569

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: