Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 169,748 T→A 5.5% intergenic (+21/‑30) fhuA → / → fhuC ferrichrome outer membrane transporter/iron(3+)‑hydroxamate import ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913169,7480TA 5.5% 388.6 / 7.0 134intergenic (+21/‑30)fhuA/fhuCferrichrome outer membrane transporter/iron(3+)‑hydroxamate import ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (61/65);  new base A (3/5);  total (64/70)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

GGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGTGACC  >  NC_000913/169599‑169897
                                                                                                                                                     |                                                                                                                                                     
 gctCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAtt                                                                                                                                                      >  1:1579457/1‑150 (MQ=255)
    ccAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTcc                                                                                                                                                   >  1:663335/1‑150 (MQ=255)
       aCGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTg                                                                                                                                                <  1:937526/150‑1 (MQ=255)
       aCGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTCATGGCTGCTTCTGGGGCGCTGAACGTCTGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTCGGGCTCGGATTACCGTg                                                                                                                                                >  2:1466753/1‑150 (MQ=255)
         gTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                                                                                                                                              <  1:593022/150‑1 (MQ=255)
         gTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                                                                                                                                              >  1:1491664/1‑150 (MQ=255)
           ggCGCTGCATGTTAACAACCTGTTCGAGCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCa                                                                                                                                            <  1:1579591/150‑1 (MQ=255)
            gcgcTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCAt                                                                                                                                           <  2:48953/150‑1 (MQ=255)
             cgcTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGATTTAACACTAATGGCTGGTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGTCACGGATTTCCGTGCCCAtt                                                                                                                                          >  1:1021368/1‑150 (MQ=255)
              gcTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAAATCCGTGccccaaa                                                                                                                                         <  2:1657524/150‑5 (MQ=255)
               cTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTc                                                                                                                                        <  2:2015374/150‑1 (MQ=255)
               cTGCATGTTAACAACCTGTTCGATCGTGAATAAGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGTCTCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAAATCCGTGccccaaaa                                                                                                                                        <  1:1811194/150‑6 (MQ=255)
                tGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGTGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCAGTTTCTAATTTCTCTTTTGGGGCACGGAAATCCGTGccccaaaag                                                                                                                                       <  1:940222/150‑7 (MQ=255)
                tGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAAATCCGTGccccaaaag                                                                                                                                       <  2:1805189/150‑7 (MQ=255)
                 gCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTcac                                                                                                                                      >  1:113902/1‑150 (MQ=255)
                  cATGTTAACAACCTGTTCGATCGTGAATACCTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTcaca                                                                                                                                     <  1:677845/150‑1 (MQ=255)
                    tGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCACCCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCAGGGATTTACGTGCCAATTTAACAAg                                                                                                                                   >  1:1058476/1‑150 (MQ=255)
                    tGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTATGGGGCACGGATTTCCGTGCCCATTTCACAAg                                                                                                                                   <  2:1223478/150‑1 (MQ=255)
                        aacaacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTTACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTgg                                                                                                                               <  1:210397/150‑1 (MQ=255)
                        aacaacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTgg                                                                                                                               >  2:1040351/1‑150 (MQ=255)
                          aaacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATATCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCt                                                                                                                             <  2:366249/149‑1 (MQ=255)
                           aatcTGTTCGATCGTTAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTg                                                                                                                            <  1:1878290/147‑1 (MQ=255)
                            acCTGTTCGCTCGGTAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCCCGGATTTCCGTGCCCATTTCACAAGTTGGCTGt                                                                                                                           <  1:1202194/150‑1 (MQ=255)
                            acCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGt                                                                                                                           <  1:1020927/150‑1 (MQ=255)
                             ccTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGtt                                                                                                                          <  2:1457590/150‑1 (MQ=255)
                              cTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTa                                                                                                                         <  2:663335/150‑1 (MQ=255)
                              cTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTa                                                                                                                         >  2:1523650/1‑150 (MQ=255)
                                gTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAACTTCGCTGTTATg                                                                                                                       >  1:311780/1‑150 (MQ=255)
                                   cGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAg                                                                                                                    <  1:1054423/150‑1 (MQ=255)
                                     aTCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGa                                                                                                                  >  1:695140/1‑150 (MQ=255)
                                     aTCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGa                                                                                                                  >  1:687993/1‑150 (MQ=255)
                                       cGTTAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAAt                                                                                                                <  2:229842/150‑1 (MQ=255)
                                          gAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAATTTGGCTGTTATGCAGGAATaaa                                                                                                             >  2:1821033/1‑148 (MQ=255)
                                                   gCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCAtt                                                                                                    >  1:1489558/1‑150 (MQ=255)
                                                     cAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTcc                                                                                                  <  1:1859958/150‑1 (MQ=255)
                                                        ctgctTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTGTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGAt                                                                                               >  2:1319967/1‑150 (MQ=255)
                                                        ctgctTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGAt                                                                                               <  1:828364/150‑1 (MQ=255)
                                                        ctgctTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGAt                                                                                               >  1:62296/1‑150 (MQ=255)
                                                        ctgctTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCAGTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGAt                                                                                               >  1:1062918/1‑150 (MQ=255)
                                                          gctTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATAc                                                                                             >  1:1977635/1‑150 (MQ=255)
                                                          gctTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATAc                                                                                             >  1:1786730/1‑150 (MQ=255)
                                                            tttAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCa                                                                                           >  2:32228/1‑150 (MQ=255)
                                                               aaCACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACtt                                                                                        >  2:909898/1‑150 (MQ=255)
                                                               aaCACTTATGGCTGCATCTGGGTCGCAGAACGTCAGGTAGTTGCAACCGCATCCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACtt                                                                                        <  1:228870/150‑1 (MQ=255)
                                                                 cacTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGAACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTACGATACCACtttt                                                                                      >  1:170827/1‑150 (MQ=255)
                                                                     tATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCAc                                                                                  <  1:1938075/150‑1 (MQ=255)
                                                                     tATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCAc                                                                                  <  1:1310263/150‑1 (MQ=255)
                                                                            gctTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTaa                                                                           <  2:1631897/150‑1 (MQ=255)
                                                                              ttCTGGGGCGCAGAACGTCAGTTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAata                                                                         <  1:1808112/150‑1 (MQ=255)
                                                                              ttCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAata                                                                         >  2:1492004/1‑150 (MQ=255)
                                                                              ttCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAata                                                                         <  2:1499097/150‑1 (MQ=255)
                                                                                cTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAatata                                                                       >  2:1747177/1‑149 (MQ=255)
                                                                                cTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATAtc                                                                       <  2:650196/150‑1 (MQ=255)
                                                                                   gggCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTcc                                                                    >  1:1314718/1‑150 (MQ=255)
                                                                                    ggCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCt                                                                   >  2:364205/1‑150 (MQ=255)
                                                                                      cgcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCttt                                                                 >  2:1187322/1‑150 (MQ=255)
                                                                                       gcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTc                                                                <  2:1685419/150‑1 (MQ=255)
                                                                                       gcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTc                                                                <  1:1204541/150‑1 (MQ=255)
                                                                                       gcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTc                                                                >  1:1425753/1‑150 (MQ=255)
                                                                                       gcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTc                                                                >  2:1448753/1‑150 (MQ=255)
                                                                                         aGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCgt                                                              >  2:647257/1‑150 (MQ=255)
                                                                                          gAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCgtg                                                             <  1:1910848/150‑1 (MQ=255)
                                                                                                  ggTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCggg                                                     >  2:1874051/1‑150 (MQ=255)
                                                                                                        tGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACg                                               >  1:811125/1‑150 (MQ=255)
                                                                                                            accgcaaccTTCCGTTTCTAATTTCTCTTTTGGGCCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGCGCCCGGGCACACgaaag                                           <  1:369046/150‑5 (MQ=255)
                                                                                                             ccgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtttt                                          >  1:450707/1‑150 (MQ=255)
                                                                                                             ccgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGAAACCACTATTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtttt                                          >  1:1539741/1‑150 (MQ=255)
                                                                                                             ccgcaaccTTCCGTTTCTAATTTCGCTTTTGCGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCATTCGTGTGCCCGGGCGCACGCtttt                                          <  1:571302/150‑1 (MQ=255)
                                                                                                              cgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTg                                         <  2:1348297/150‑1 (MQ=255)
                                                                                                              cgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTg                                         <  1:436677/150‑1 (MQ=255)
                                                                                                              cgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGGGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTg                                         <  2:425055/150‑1 (MQ=255)
                                                                                                              cgcaaccTTCCGTTGCTAATTGCTCTTTTGGGGCACCGATGTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGCGCCCGGGCACACgaaagga                                         <  1:1540262/150‑7 (MQ=255)
                                                                                                               gcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGc                                        <  2:1251942/150‑1 (MQ=255)
                                                                                                                caaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCa                                       <  1:1298000/150‑1 (MQ=255)
                                                                                                                     ttCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGc                                  <  1:658679/150‑1 (MQ=255)
                                                                                                                      tCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCt                                 <  1:1040351/150‑1 (MQ=255)
                                                                                                                       ccGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTg                                >  2:133299/1‑150 (MQ=255)
                                                                                                                             cTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTa                          >  2:1145032/1‑150 (MQ=255)
                                                                                                                             cTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTa                          >  1:1725696/1‑150 (MQ=255)
                                                                                                                               aaTTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAAc                        >  2:1715751/1‑150 (MQ=255)
                                                                                                                               aaTTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTGTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAAc                        >  2:1707506/1‑150 (MQ=255)
                                                                                                                               aaTTTCTCTTTGGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATGATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAAc                        <  1:1722372/150‑1 (MQ=255)
                                                                                                                                    ctctTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTc                   <  2:1579457/150‑1 (MQ=255)
                                                                                                                                    ctctTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTc                   >  1:1657004/1‑150 (MQ=255)
                                                                                                                                     tctTTTGGGGCACGGATTTCCGTGCCCATTTCACATGTTGGCTGTTATGCAGGAATACACGATTCATTCCGATAGCACTTTTGCACTGCGTAATATCTGCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTcc                  <  2:935586/150‑1 (MQ=255)
                                                                                                                                     tctTTTGGGGCACGGATTTCCGTGCCCAGTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTcc                  >  2:1274489/1‑150 (MQ=255)
                                                                                                                                        tttGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGc               <  2:113902/150‑1 (MQ=255)
                                                                                                                                         ttGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGcc              <  1:1120770/150‑1 (MQ=255)
                                                                                                                                           ggggCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCTATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCgg            <  2:1786730/150‑1 (MQ=255)
                                                                                                                                           ggggCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCgg            <  1:1492004/150‑1 (MQ=255)
                                                                                                                                           ggggCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTAGCATCCGATGTCGTTAACCTTTCCTGCCgg            >  1:153678/1‑150 (MQ=255)
                                                                                                                                            gggCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCggg           <  1:1876653/150‑1 (MQ=255)
                                                                                                                                            gggCACGGAAATCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACGCGAATCCTTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCggg           >  2:1097835/1‑150 (MQ=255)
                                                                                                                                             ggCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGa          >  1:1188105/1‑150 (MQ=255)
                                                                                                                                             ggCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGa          >  1:794554/1‑150 (MQ=255)
                                                                                                                                                 cGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGt      >  2:1652457/1‑150 (MQ=255)
                                                                                                                                                 cGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGt      >  2:495025/1‑150 (MQ=255)
                                                                                                                                                 cGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGt      >  1:1432848/1‑150 (MQ=255)
                                                                                                                                                  ggATTTCCGTGCCCATTTCACAAGTTTGCTGTTTTGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCCGGCAGGAAAggt     <  1:1703062/150‑3 (MQ=255)
                                                                                                                                                     tttCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGTGAcc  <  1:364205/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGTGACC  >  NC_000913/169599‑169897

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: