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‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATA  >  NC_000913/4470809‑4471106
                                                                                                                                                       |                                                                                                                                                     
tACCGGGATCTGACAGGAGGTGATGATCATATTGCCCTGATCAACGCCCTGTACGTAAGGACCGAAAGCTGCCGGTGCATTTTCCTTCGAGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAgg                                                                                                                                                      <  2:894688/150‑1 (MQ=255)
 aCCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGTACCGATAGCTTCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGAATTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAggg                                                                                                                                                     >  2:276887/1‑150 (MQ=255)
 aCCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAggg                                                                                                                                                     >  1:1635021/1‑150 (MQ=255)
  ccGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGa                                                                                                                                                    >  2:369358/1‑150 (MQ=255)
  ccGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGa                                                                                                                                                    >  2:3446/1‑150 (MQ=255)
   cGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAg                                                                                                                                                   >  2:377175/1‑150 (MQ=255)
    gggATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGc                                                                                                                                                  <  2:1244833/150‑1 (MQ=255)
        tCTGACAGGAGGTGATGATCATATTGCCCAGATAAACTCCCTGTACGTACTGAGCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTAATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGGTCCCTATAATGAGCCAGa                                                                                                                                              <  2:1006818/150‑1 (MQ=255)
         cTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGAc                                                                                                                                             <  1:1729855/150‑1 (MQ=255)
          tGACCGGAGGTGATGATCATATTGCCCAGATCATCGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGAcc                                                                                                                                            <  1:595923/150‑2 (MQ=255)
          tGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGACt                                                                                                                                            <  1:770954/150‑1 (MQ=255)
          tGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACt                                                                                                                                            >  1:127663/1‑150 (MQ=255)
           gACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCCCTATAATGAGCCAGACtt                                                                                                                                           <  2:703943/150‑1 (MQ=255)
           gACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACtt                                                                                                                                           <  2:1421737/150‑1 (MQ=255)
             ccGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAATAAAAGTCTGGCTCCCTATAATGAGCCAGACtt                                                                                                                                           <  1:815744/150‑1 (MQ=38)
              cGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTtta                                                                                                                                        <  1:526206/150‑1 (MQ=255)
               ggAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTttat                                                                                                                                       <  2:389583/150‑1 (MQ=255)
                gAGGTGATGATCATATTGCCCATATCAACGCCCTGTTAGTAATGACCTATAGCTGCCGGTGAATTTTACGTCGCTATAGTTTTGATCATGATTTCTCCTTTATTACAGCTGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatt                                                                                                                                      <  2:324295/150‑1 (MQ=255)
                 cggTGATGATCATAGTGCCCAGGTCATCTCCCTGTACGTAAGGACCGATAGCTGCCGGTTCATTTTCCGTCGCGTTAGTTTTGCTCATGAGTTCTCCTTAATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatta                                                                                                                                     <  1:208496/149‑1 (MQ=255)
                 aGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatta                                                                                                                                     >  1:230038/1‑150 (MQ=255)
                   gTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAcc                                                                                                                                   <  2:1329715/150‑1 (MQ=255)
                   gTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAcc                                                                                                                                   <  1:3446/150‑1 (MQ=255)
                    tgatgaTCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCa                                                                                                                                  <  1:1348910/150‑1 (MQ=255)
                    tgatgaTCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCa                                                                                                                                  >  2:1723736/1‑150 (MQ=255)
                    tgatgaTCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCa                                                                                                                                  >  2:1485087/1‑150 (MQ=255)
                       tgatCATATTGCCCAGATCAAGTCCGTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTAGCGATAGTTTTGCTAATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCGTCATTATAGGGAGCCAGACTTTTATTACCAAc                                                                                                                                <  2:487854/150‑1 (MQ=255)
                       tgatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTCTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAcc                                                                                                                               >  2:1593386/1‑150 (MQ=255)
                        gatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCg                                                                                                                              >  1:686341/1‑150 (MQ=255)
                        gatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCg                                                                                                                              >  1:552653/1‑150 (MQ=255)
                         ttCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGc                                                                                                                             <  2:947816/149‑1 (MQ=255)
                          tCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCa                                                                                                                            >  1:1254850/1‑150 (MQ=255)
                           cATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTTTTACCAACCGCaa                                                                                                                           >  2:886346/1‑150 (MQ=255)
                                 gCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGTGCGCCAGACTTTTATTACCTACCGCaattaatt                                                                                                                     >  1:75540/1‑150 (MQ=255)
                                 gCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                                                                                     >  1:472612/1‑150 (MQ=255)
                                  cccAGATCAACGCCCTGTACGTAAGGCCCGATAGCTGCCGGTGCTTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGCGCCAGACTTTTATTACCAACCGCATTTAATTg                                                                                                                    >  2:161161/1‑150 (MQ=255)
                                   ccAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                                                                                   >  2:1190137/1‑150 (MQ=255)
                                     aGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGcc                                                                                                                 <  1:958919/150‑1 (MQ=255)
                                        tCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGc                                                                                                              >  2:647644/1‑150 (MQ=255)
                                         cAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGca                                                                                                             <  1:920057/150‑1 (MQ=255)
                                            cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcacca                                                                                                          >  1:182896/1‑150 (MQ=255)
                                            cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcacca                                                                                                          >  1:456372/1‑150 (MQ=255)
                                             gCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAATCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcaccac                                                                                                         <  1:1573834/150‑1 (MQ=255)
                                              cccTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACcaca                                                                                                        >  2:810191/1‑150 (MQ=255)
                                                  gTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACAACATTAt                                                                                                    >  1:1292758/1‑150 (MQ=255)
                                                   ttagTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTATATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATg                                                                                                   <  1:659701/147‑1 (MQ=255)
                                                    aCGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                                                                  >  1:1408741/1‑150 (MQ=255)
                                                    aCGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                                                                  <  1:1232873/150‑1 (MQ=255)
                                                        aaGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaa                                                                                              >  1:731235/1‑150 (MQ=255)
                                                        aaGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaa                                                                                              >  2:718971/1‑150 (MQ=255)
                                                          ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAc                                                                                            >  1:715844/1‑150 (MQ=255)
                                                          ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAc                                                                                            >  2:1389199/1‑150 (MQ=255)
                                                            aCCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActc                                                                                          >  2:1094192/1‑150 (MQ=255)
                                                                  tgCTGCCGGGGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                                                                                    <  2:79364/149‑1 (MQ=255)
                                                                   gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                                                                                   <  2:1635021/150‑1 (MQ=255)
                                                                    cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                                                                                  >  2:880624/1‑150 (MQ=255)
                                                                     tGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcaca                                                                                 <  1:1603075/150‑1 (MQ=255)
                                                                          gTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAtt                                                                            >  2:396081/1‑150 (MQ=255)
                                                                            gCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                                                                          <  2:405748/150‑1 (MQ=255)
                                                                                     gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgc                                                                 <  1:1663118/150‑1 (MQ=255)
                                                                                        gcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg                                                              >  1:1684607/1‑150 (MQ=255)
                                                                                          gATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACGTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAACCACTTTTTCACATTATTTGAATTTGAGCGCGat                                                            >  2:647694/1‑150 (MQ=255)
                                                                                          gATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTAGAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGAACCACATTATGGGAAAACTCTTTTTAACAGTATTTGAATTTGAGCGCGat                                                            >  2:1500691/1‑150 (MQ=255)
                                                                                                 ttGCTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTg                                                     <  2:1418539/150‑1 (MQ=255)
                                                                                                    cTCATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGgcg                                                  >  2:1452621/1‑150 (MQ=255)
                                                                                                      cATGATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGcgcg                                                >  2:743358/1‑150 (MQ=255)
                                                                                                         gATTTCTCCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGttt                                             <  1:886346/150‑1 (MQ=255)
                                                                                                             tctcCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGc                                         >  2:1192043/1‑150 (MQ=255)
                                                                                                             tctcCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGc                                         <  2:1346021/150‑1 (MQ=255)
                                                                                                             tctcCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGc                                         >  2:181360/1‑150 (MQ=255)
                                                                                                              ctcCTTTATTACAGCGGTA‑‑AAAAGGCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAGTTGAGCGCGATAGCATTGGCGCGTTTTCGCa                                        <  1:834831/150‑1 (MQ=255)
                                                                                                               tcCTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGGTTTCGCAc                                       <  2:456372/150‑1 (MQ=255)
                                                                                                                ccTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGATCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACg                                      <  1:369358/150‑1 (MQ=255)
                                                                                                                ccTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCCGACTTTTATTACCAACCGCAATTAATTGGCCCGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGAGCGATATCATTGGCGCGTTTTCGCACg                                      >  1:1473140/1‑150 (MQ=255)
                                                                                                                 cTTTATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACgg                                     >  2:183702/1‑150 (MQ=255)
                                                                                                                  tttATTACAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGc                                    <  1:1300435/150‑1 (MQ=255)
                                                                                                                    tattaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGAGCGATATCATTGGCGCGTTTTCGCACGGCaa                                  >  1:1645325/1‑150 (MQ=255)
                                                                                                                     attaCAGCGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTATTTGGCCAGCACCACATTATGGGAAAACTCTTTATCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCaaa                                 <  1:276887/150‑1 (MQ=255)
                                                                                                                            cGGTA‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGAt                          <  2:127663/150‑1 (MQ=255)
                                                                                                                                    aaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAAcc                    >  2:1315218/1‑150 (MQ=255)
                                                                                                                                    aaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAAcc                    >  2:1530645/1‑150 (MQ=255)
                                                                                                                                     aaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCg                   >  1:1186017/1‑150 (MQ=255)
                                                                                                                                      aGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                  >  1:55382/1‑150 (MQ=255)
                                                                                                                                      aGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                  >  1:108346/1‑150 (MQ=255)
                                                                                                                                       gTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGt                 >  1:1731579/1‑150 (MQ=255)
                                                                                                                                       gTCTGGCTCCCTATAAGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTGTCGCACGGCAAAGCTGGATGAAACCGGt                 >  1:1429967/1‑150 (MQ=255)
                                                                                                                                       gTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGt                 <  1:1389199/150‑1 (MQ=255)
                                                                                                                                       gTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGt                 >  1:1081600/1‑150 (MQ=255)
                                                                                                                                         cTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTc               <  1:183702/150‑1 (MQ=255)
                                                                                                                                         cTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTc               >  1:1077253/1‑150 (MQ=255)
                                                                                                                                          tGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCg              <  1:1452621/150‑1 (MQ=255)
                                                                                                                                           ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCgg             >  2:430413/1‑150 (MQ=255)
                                                                                                                                            gCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGc            >  1:35675/1‑150 (MQ=255)
                                                                                                                                             cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCa           >  1:1690520/1‑150 (MQ=255)
                                                                                                                                             cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCa           >  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‑‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATA  >  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: