Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,957 G→A 7.7% intergenic (‑44/+29) ridA ← / ← pyrI enamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,470,9570GA 7.7% 392.1 / 14.6 130intergenic (‑44/+29)ridA/pyrIenamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base G (60/60);  new base A (6/4);  total (66/64)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

Only 100 of 135 total aligned reads displayed.

TACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATA  >  NC_000913/4470809‑4471106
                                                                                                                                                       |                                                                                                                                                     
tACCGGGATCTGACAGGAGGTGATGATCATATTGCCCTGATCAACGCCCTGTACGTAAGGACCGAAAGCTGCCGGTGCATTTTCCTTCGAGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAgg                                                                                                                                                      <  2:894688/150‑1 (MQ=255)
 aCCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGTACCGATAGCTTCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGAATTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAggg                                                                                                                                                     >  2:276887/1‑150 (MQ=255)
 aCCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAggg                                                                                                                                                     >  1:1635021/1‑150 (MQ=255)
  ccGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGa                                                                                                                                                    >  2:369358/1‑150 (MQ=255)
  ccGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGa                                                                                                                                                    >  2:3446/1‑150 (MQ=255)
   cGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAg                                                                                                                                                   >  2:377175/1‑150 (MQ=255)
    gggATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGc                                                                                                                                                  <  2:1244833/150‑1 (MQ=255)
        tCTGACAGGAGGTGATGATCATATTGCCCAGATAAACTCCCTGTACGTACTGAGCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTAATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGG‑TCCCTATAATGAGCCAGa                                                                                                                                              <  2:1006818/150‑1 (MQ=255)
         cTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGAc                                                                                                                                             <  1:1729855/150‑1 (MQ=255)
          tGACCGGAGGTGATGATCATATTGCCCAGATCATCGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGAcc                                                                                                                                            <  1:595923/150‑2 (MQ=255)
          tGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCCCTATAATGAGCCAGACt                                                                                                                                            <  1:770954/150‑1 (MQ=255)
          tGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACt                                                                                                                                            >  1:127663/1‑150 (MQ=255)
           gACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCCCTATAATGAGCCAGACtt                                                                                                                                           <  2:703943/150‑1 (MQ=255)
           gACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACtt                                                                                                                                           <  2:1421737/150‑1 (MQ=255)
             ccGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAATAAAAGTCTGGC‑TCCCTATAATGAGCCAGACtt                                                                                                                                           <  1:815744/150‑1 (MQ=38)
              cGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTtta                                                                                                                                        <  1:526206/150‑1 (MQ=255)
               ggAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTttat                                                                                                                                       <  2:389583/150‑1 (MQ=255)
                gAGGTGATGATCATATTGCCCATATCAACGCCCTGTTAGTAATGACCTATAGCTGCCGGTGAATTTTACGTCGCTATAGTTTTGATCATGATTTCTCCTTTATTACAGCTGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTttatt                                                                                                                                      <  2:324295/150‑1 (MQ=255)
                 cggTGATGATCATAGTGCCCAGGTCATCTCCCTGTACGTAAGGACCGATAGCTGCCGGTTCATTTTCCGTCGCGTTAGTTTTGCTCATGAGTTCTCCTTAATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTttatta                                                                                                                                     <  1:208496/149‑1 (MQ=255)
                 aGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTttatta                                                                                                                                     >  1:230038/1‑150 (MQ=255)
                   gTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTAcc                                                                                                                                   <  2:1329715/150‑1 (MQ=255)
                   gTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTAcc                                                                                                                                   <  1:3446/150‑1 (MQ=255)
                    tgatgaTCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCa                                                                                                                                  <  1:1348910/150‑1 (MQ=255)
                    tgatgaTCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCa                                                                                                                                  >  2:1723736/1‑150 (MQ=255)
                    tgatgaTCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCa                                                                                                                                  >  2:1485087/1‑150 (MQ=255)
                       tgatCATATTGCCCAGATCAAGTCCGTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTAGCGATAGTTTTGCTAATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCGTCATTATAGGGAGCCAGACTTTTATTACCAAc                                                                                                                                <  2:487854/150‑1 (MQ=255)
                       tgatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTCTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAAcc                                                                                                                               >  2:1593386/1‑150 (MQ=255)
                        gatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCg                                                                                                                              >  1:686341/1‑150 (MQ=255)
                        gatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCg                                                                                                                              >  1:552653/1‑150 (MQ=255)
                         ttCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGc                                                                                                                             <  2:947816/149‑1 (MQ=255)
                          tCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCa                                                                                                                            >  1:1254850/1‑150 (MQ=255)
                           cATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTTTTACCAACCGCaa                                                                                                                           >  2:886346/1‑150 (MQ=255)
                                 gCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGTGCGCCAGACTTTTATTACCTACCGCaattaatt                                                                                                                     >  1:75540/1‑150 (MQ=255)
                                 gCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                                                                                     >  1:472612/1‑150 (MQ=255)
                                  cccAGATCAACGCCCTGTACGTAAGGCCCGATAGCTGCCGGTGCTTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGCGCCAGACTTTTATTACCAACCGCATTTAATTg                                                                                                                    >  2:161161/1‑150 (MQ=255)
                                   ccAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                                                                                   >  2:1190137/1‑150 (MQ=255)
                                     aGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGcc                                                                                                                 <  1:958919/150‑1 (MQ=255)
                                        tCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGc                                                                                                              >  2:647644/1‑150 (MQ=255)
                                         cAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGca                                                                                                             <  1:920057/150‑1 (MQ=255)
                                            cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcacca                                                                                                          >  1:182896/1‑150 (MQ=255)
                                            cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcacca                                                                                                          >  1:456372/1‑150 (MQ=255)
                                             gCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAATCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcaccac                                                                                                         <  1:1573834/150‑1 (MQ=255)
                                              cccTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACcaca                                                                                                        >  2:810191/1‑150 (MQ=255)
                                                  gTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACAACATTAt                                                                                                    >  1:1292758/1‑150 (MQ=255)
                                                   ttagTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTATATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATg                                                                                                   <  1:659701/147‑1 (MQ=255)
                                                    aCGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                                                                  >  1:1408741/1‑150 (MQ=255)
                                                    aCGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                                                                  <  1:1232873/150‑1 (MQ=255)
                                                        aaGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaa                                                                                              >  1:731235/1‑150 (MQ=255)
                                                        aaGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaa                                                                                              >  2:718971/1‑150 (MQ=255)
                                                          ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAc                                                                                            >  1:715844/1‑150 (MQ=255)
                                                          ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAc                                                                                            >  2:1389199/1‑150 (MQ=255)
                                                            aCCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActc                                                                                          >  2:1094192/1‑150 (MQ=255)
                                                                  tgCTGCCGGGGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                                                                                    <  2:79364/149‑1 (MQ=255)
                                                                   gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                                                                                   <  2:1635021/150‑1 (MQ=255)
                                                                    cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                                                                                  >  2:880624/1‑150 (MQ=255)
                                                                     tGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcaca                                                                                 <  1:1603075/150‑1 (MQ=255)
                                                                          gTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAtt                                                                            >  2:396081/1‑150 (MQ=255)
                                                                            gCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                                                                          <  2:405748/150‑1 (MQ=255)
                                                                                     gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgc                                                                 <  1:1663118/150‑1 (MQ=255)
                                                                                        gcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg                                                              >  1:1684607/1‑150 (MQ=255)
                                                                                          gATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACGTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAACCACTTTTTCACATTATTTGAATTTGAGCGCGat                                                            >  2:647694/1‑150 (MQ=255)
                                                                                          gATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTAGAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGAACCACATTATGGGAAAACTCTTTTTAACAGTATTTGAATTTGAGCGCGat                                                            >  2:1500691/1‑150 (MQ=255)
                                                                                                 ttGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTg                                                     <  2:1418539/150‑1 (MQ=255)
                                                                                                    cTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGgcg                                                  >  2:1452621/1‑150 (MQ=255)
                                                                                                      cATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGcgcg                                                >  2:743358/1‑150 (MQ=255)
                                                                                                         gATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGttt                                             <  1:886346/150‑1 (MQ=255)
                                                                                                             tctcCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGc                                         >  2:1192043/1‑150 (MQ=255)
                                                                                                             tctcCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGc                                         <  2:1346021/150‑1 (MQ=255)
                                                                                                             tctcCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGc                                         >  2:181360/1‑150 (MQ=255)
                                                                                                              ctcCTTTATTACAGCGGTA‑‑AAAAGGCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAGTTGAGCGCGATAGCATTGGCGCGTTTTCGCa                                        <  1:834831/150‑1 (MQ=255)
                                                                                                               tcCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGGTTTCGCAc                                       <  2:456372/150‑1 (MQ=255)
                                                                                                                ccTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGATCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACg                                      <  1:369358/150‑1 (MQ=255)
                                                                                                                ccTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCCGACTTTTATTACCAACCGCAATTAATTGGCCCGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGAGCGATATCATTGGCGCGTTTTCGCACg                                      >  1:1473140/1‑150 (MQ=255)
                                                                                                                 cTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACgg                                     >  2:183702/1‑150 (MQ=255)
                                                                                                                  tttATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGc                                    <  1:1300435/150‑1 (MQ=255)
                                                                                                                    tattaCAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGAGCGATATCATTGGCGCGTTTTCGCACGGCaa                                  >  1:1645325/1‑150 (MQ=255)
                                                                                                                     attaCAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTATTTGGCCAGCACCACATTATGGGAAAACTCTTTATCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCaaa                                 <  1:276887/150‑1 (MQ=255)
                                                                                                                            cGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGAt                          <  2:127663/150‑1 (MQ=255)
                                                                                                                                    aaaGTCTGGC‑TCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAAcc                    >  2:1315218/1‑150 (MQ=255)
                                                                                                                                    aaaGTCTGGC‑TCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAAcc                    >  2:1530645/1‑150 (MQ=255)
                                                                                                                                     aaGTCTGGC‑TCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCg                   >  1:1186017/1‑150 (MQ=255)
                                                                                                                                      aGTCTGGC‑TCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                  >  1:55382/1‑150 (MQ=255)
                                                                                                                                      aGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                  >  1:108346/1‑150 (MQ=255)
                                                                                                                                       gTCTGGC‑TCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGt                 >  1:1731579/1‑150 (MQ=255)
                                                                                                                                       gTCTGGC‑TCCCTATAAGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTGTCGCACGGCAAAGCTGGATGAAACCGGt                 >  1:1429967/1‑150 (MQ=255)
                                                                                                                                       gTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGt                 <  1:1389199/150‑1 (MQ=255)
                                                                                                                                       gTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGt                 >  1:1081600/1‑150 (MQ=255)
                                                                                                                                         cTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTc               <  1:183702/150‑1 (MQ=255)
                                                                                                                                         cTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTc               >  1:1077253/1‑150 (MQ=255)
                                                                                                                                          tGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCg              <  1:1452621/150‑1 (MQ=255)
                                                                                                                                           ggC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCgg             >  2:430413/1‑150 (MQ=255)
                                                                                                                                            gC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGc            >  1:35675/1‑150 (MQ=255)
                                                                                                                                             c‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCa           >  1:1690520/1‑150 (MQ=255)
                                                                                                                                             c‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCa           >  1:1175743/1‑150 (MQ=255)
                                                                                                                                                cATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATg         <  1:1180407/150‑1 (MQ=255)
                                                                                                                                                 aTTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATgg        >  1:12307/1‑150 (MQ=255)
                                                                                                                                                 aTTATAGGGAGCCAGACTTATATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATgg        <  2:935170/150‑1 (MQ=255)
                                                                                                                                                  ttATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGc       <  2:788103/150‑1 (MQ=255)
                                                                                                                                                      aGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGAt   <  1:718971/150‑1 (MQ=255)
                                                                                                                                                       gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATa  <  2:46372/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
TACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGC‑TCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATA  >  NC_000913/4470809‑4471106

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑