Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,310,974 T→A 6.1% T3S (ACG→TCG) ‡ alsA ← D‑allose ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,310,9740TA 6.1% 272.1 / 4.8 91T3S (ACG→TCG) ‡alsAD‑allose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (38/47);  new base A (2/4);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGATAAAACTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGG  >  NC_000913/4310825‑4311122
                                                                                                                                                     |                                                                                                                                                    
ggATAAAACTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGt                                                                                                                                                      <  1:841988/150‑1 (MQ=255)
  aTAAAACTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTgg                                                                                                                                                    <  1:159639/150‑1 (MQ=255)
   tAAAACTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGc                                                                                                                                                   <  1:913316/150‑1 (MQ=255)
    aaaaCTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGAGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGcc                                                                                                                                                  >  1:1473794/1‑150 (MQ=255)
     aaaCTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCTACTTTAATGCGTGAACCCGACCAAAGGACTTGCCGATCACCGCCATCGATATATATGGCGTGGCCa                                                                                                                                                 >  2:871328/1‑150 (MQ=255)
     aaaCTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCa                                                                                                                                                 <  2:1152500/150‑1 (MQ=255)
     aaaCTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCa                                                                                                                                                 >  1:1120625/1‑150 (MQ=255)
      aaCTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCat                                                                                                                                                >  2:56116/1‑150 (MQ=255)
      aaCTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTACTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCat                                                                                                                                                <  2:859626/150‑1 (MQ=255)
       aCTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCata                                                                                                                                               <  2:21393/150‑1 (MQ=255)
        cTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCataa                                                                                                                                              >  1:706761/1‑150 (MQ=255)
         tttCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATTTATGGCCACGCCatata                                                                                                                                             <  1:1366115/150‑3 (MQ=255)
         tttCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCCACGCCatata                                                                                                                                             <  1:158435/150‑3 (MQ=255)
          ttCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCataaat                                                                                                                                            >  2:564584/1‑150 (MQ=255)
              ttAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAAcc                                                                                                                                        <  1:467223/150‑1 (MQ=255)
                aGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACctc                                                                                                                                      >  2:520298/1‑150 (MQ=255)
                aGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCCACGCCatatatatcgat                                                                                                                                      <  2:1327448/150‑6 (MQ=255)
                  cGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATTTATGGCCACGCCatatatatcgatgg                                                                                                                                    <  1:1500057/150‑8 (MQ=255)
                  cGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACctctc                                                                                                                                    >  2:757142/1‑150 (MQ=255)
                    tGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTg                                                                                                                                  >  2:600321/1‑150 (MQ=255)
                    tGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTg                                                                                                                                  <  2:101086/150‑1 (MQ=255)
                       aTTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCAc                                                                                                                               <  1:1477693/150‑1 (MQ=255)
                         ttACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGa                                                                                                                             <  2:229360/150‑1 (MQ=255)
                          tACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGaa                                                                                                                            >  1:307739/1‑150 (MQ=255)
                          tACCCGCGCCAGTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGaa                                                                                                                            >  2:1313626/1‑150 (MQ=255)
                           aCCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTTAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGaaa                                                                                                                           >  2:916156/1‑150 (MQ=255)
                           aCCCGCGCCA‑TTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAat                                                                                                                          >  1:37458/1‑150 (MQ=255)
                                cgcCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCg                                                                                                                      <  2:1472232/150‑1 (MQ=255)
                                 gcCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGt                                                                                                                     >  2:482561/1‑150 (MQ=255)
                                   cATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGtat                                                                                                                   <  1:1372308/150‑1 (MQ=255)
                                     ttttCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGtatta                                                                                                                 <  2:896777/150‑1 (MQ=255)
                                       ttCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTattatt                                                                                                               <  1:1462598/150‑1 (MQ=255)
                                        tctcCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTg                                                                                                              <  2:260747/150‑1 (MQ=255)
                                          tcCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGcc                                                                                                            <  2:198796/150‑1 (MQ=255)
                                           tcTATTATTGCTTTTATTTCACCAGGATAAACCTTTATTTTAACGGACTTTAATGCGTGAACCTGACCATAGGACTTGCCGATCCCCGCCTTCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACTAAATATCTTATTATTGctt                                                                                                           <  2:1473794/149‑3 (MQ=255)
                                             tAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGa                                                                                                         >  2:1478491/1‑150 (MQ=255)
                                                   tGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTa                                                                                                   >  1:1442458/1‑150 (MQ=255)
                                                     cATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACa                                                                                                 >  2:630947/1‑150 (MQ=255)
                                                      aTGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACAt                                                                                                >  1:1032758/1‑150 (MQ=255)
                                                         tATTTCAGCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTa                                                                                             <  2:896100/150‑1 (MQ=255)
                                                         tATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTa                                                                                             >  2:74072/1‑150 (MQ=255)
                                                         tATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTa                                                                                             <  2:1316371/150‑1 (MQ=255)
                                                              cACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAgg                                                                                        <  1:757142/150‑1 (MQ=255)
                                                                  aGGATAAACCGTTAAATTAACCGACTTTTATGCGTGAACCGTACCAAAGGACTTGCCGATCACCTCCATCGATATATATTGCGTGGCCATAAATAAGCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTa                                                                                    <  2:933908/150‑1 (MQ=255)
                                                                             ttAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTc                                                                               <  2:893153/144‑1 (MQ=255)
                                                                             ttAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTc                                                                               >  1:893153/1‑144 (MQ=255)
                                                                              tAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTAt                                                                        >  1:1112671/1‑150 (MQ=255)
                                                                               aaaTTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTAtt                                                                       >  1:1142754/1‑150 (MQ=255)
                                                                               aaaTTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGAATTATTGCATGATGGCTACATTTATCAGGTCTACGTTCAGGTTAtt                                                                       >  2:478851/1‑150 (MQ=255)
                                                                                aaTTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTAttt                                                                      >  2:1421015/1‑150 (MQ=255)
                                                                                 aTTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACAttt                                                                                              <  2:1508249/125‑1 (MQ=255)
                                                                                 aTTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACAttt                                                                                              >  1:1508249/1‑125 (MQ=255)
                                                                                    aaCCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTttgtt                                                                  <  1:195674/150‑1 (MQ=255)
                                                                                     aCCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTttgttg                                                                 <  1:396072/150‑1 (MQ=255)
                                                                                      ccGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAgg                                                                                        <  1:75712/126‑1 (MQ=255)
                                                                                      ccGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAgg                                                                                        >  2:75712/1‑126 (MQ=255)
                                                                                       cGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGc                                                               <  2:733164/150‑1 (MQ=255)
                                                                                          cTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCgg                                                            <  2:225960/150‑1 (MQ=255)
                                                                                          cTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCgg                                                            <  1:1024050/150‑1 (MQ=255)
                                                                                            ttAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCg                                                             <  1:425612/147‑1 (MQ=255)
                                                                                            ttAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCg                                                             >  2:425612/1‑147 (MQ=255)
                                                                                            ttAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGAt                                                          <  2:307739/150‑1 (MQ=255)
                                                                                                 gCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGac                                                     <  1:564584/150‑1 (MQ=255)
                                                                                                 gCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGac                                                     >  2:1012713/1‑150 (MQ=255)
                                                                                                    tGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACt                                                  >  1:106827/1‑150 (MQ=255)
                                                                                                         cGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCAtt                                             <  2:1120625/150‑1 (MQ=255)
                                                                                                          ggACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTa                                            <  2:634845/150‑1 (MQ=255)
                                                                                                           gACCAAAGGACTTTCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAg                                           <  2:566588/150‑1 (MQ=255)
                                                                                                           gACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAg                                           <  1:317654/150‑1 (MQ=255)
                                                                                                             ccAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGca                                         >  1:934161/1‑150 (MQ=255)
                                                                                                               aaaGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGcatc                                       >  1:816456/1‑150 (MQ=255)
                                                                                                                aaGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAAACTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGcatca                                      >  2:869484/1‑150 (MQ=255)
                                                                                                                 aGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGcatcat                                     <  2:37458/150‑1 (MQ=255)
                                                                                                                 aGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGcatcat                                     <  1:56116/150‑1 (MQ=255)
                                                                                                                        gCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCgg                              <  2:1009091/150‑1 (MQ=255)
                                                                                                                         ccGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCTTATTATTGCCTGATGGCTACATTTATCAGGGCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGc                             <  1:452540/150‑1 (MQ=255)
                                                                                                                           gATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGcat                           >  1:1108948/1‑150 (MQ=255)
                                                                                                                             gccacGCCATATATATCGATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGcatca                         >  2:224815/5‑150 (MQ=255)
                                                                                                                               cccGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCAtt                       <  1:1523163/150‑1 (MQ=255)
                                                                                                                                     aTCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGttt                 <  2:1112671/150‑1 (MQ=255)
                                                                                                                                        gttatttatGGCCACGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTAt              >  2:1123583/3‑150 (MQ=255)
                                                                                                                                          tatttatGGCCACGCCATATATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTg            >  2:737293/1‑150 (MQ=255)
                                                                                                                                          tatatatGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTg            <  1:600321/150‑1 (MQ=255)
                                                                                                                                          tatatatGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTg            <  2:1304215/150‑1 (MQ=255)
                                                                                                                                          tatatatGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTg            <  1:482561/150‑1 (MQ=255)
                                                                                                                                            tatatGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAg          >  1:1493015/1‑150 (MQ=255)
                                                                                                                                            tatatGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAg          <  1:81780/150‑1 (MQ=255)
                                                                                                                                             atatGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGt         >  1:905144/1‑150 (MQ=255)
                                                                                                                                              tatGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTg        <  1:869484/150‑1 (MQ=255)
                                                                                                                                                 ggCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATAAGCATCATCATCCGGCATCATTCAGGTTTATTGAGTGAcc     >  2:1444856/1‑150 (MQ=255)
                                                                                                                                                   cGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAg   >  2:1058126/1‑150 (MQ=255)
                                                                                                                                                    gTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAgg  <  2:816456/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGATAAAACTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGG  >  NC_000913/4310825‑4311122

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