Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,465,671 C→G 5.9% pseudogene (280/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,465,6710CG 5.9% 373.7 / 13.0 119pseudogene (280/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (57/54);  new base G (3/4);  total (60/59)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

Only 100 of 138 total aligned reads displayed.

ACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAAT  >  NC_000913/1465522‑1465820
                                                                                                                                                     |                                                                                                                                                     
aCCCAGGAAAACTCAGTCCGGAAGAACTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatc                                                                                                                                                       >  2:1545309/1‑150 (MQ=255)
 cccAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatca                                                                                                                                                      <  2:1185047/150‑1 (MQ=255)
 cccAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatca                                                                                                                                                      >  1:1521846/1‑150 (MQ=255)
  ccAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcat                                                                                                                                                     <  1:1227186/150‑1 (MQ=255)
  ccAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCGCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATTACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcat                                                                                                                                                     <  1:287199/150‑1 (MQ=255)
   cAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcatc                                                                                                                                                    <  2:481555/150‑1 (MQ=255)
   cAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGACATGGCTTTCCACCGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcatc                                                                                                                                                    >  2:796964/1‑150 (MQ=255)
    aGGAAAACTCAGTCCGGAAGTGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcatca                                                                                                                                                   <  1:1452475/150‑1 (MQ=255)
     ggAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcatcat                                                                                                                                                  >  2:1520770/1‑150 (MQ=255)
     ggAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcatcat                                                                                                                                                  <  2:649478/150‑1 (MQ=255)
          aCTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATGATCATCATCg                                                                                                                                             >  2:914703/1‑150 (MQ=255)
           cTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATGATCATCATCgt                                                                                                                                            <  2:1828331/150‑2 (MQ=255)
            tCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATGATCATCATCgtc                                                                                                                                           <  1:1686931/150‑3 (MQ=255)
             cAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGcaa                                                                                                                                          >  2:1592606/1‑150 (MQ=255)
              aGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATGATCATCATCgtcgt                                                                                                                                         <  1:998326/150‑5 (MQ=255)
              aGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGcaac                                                                                                                                         <  1:808159/150‑1 (MQ=255)
               gTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGcaaca                                                                                                                                        >  1:1284303/1‑150 (MQ=255)
                 ccGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAAt                                                                                                                                      >  1:790584/1‑150 (MQ=255)
                 ccGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAAt                                                                                                                                      >  1:1569479/1‑150 (MQ=255)
                 ccGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAAt                                                                                                                                      <  2:761751/150‑1 (MQ=255)
                     aaGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATtctc                                                                                                                                  <  1:1215215/150‑1 (MQ=255)
                      agagCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTcc                                                                                                                                 <  1:782826/150‑1 (MQ=255)
                      agagCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTcc                                                                                                                                 <  1:166242/150‑1 (MQ=255)
                      agagCTTAAACGCTTACCCTCTGAATGCCCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTcc                                                                                                                                 <  2:619629/150‑1 (MQ=255)
                          cTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTc                                                                                                                             >  1:205281/1‑150 (MQ=255)
                          cTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTc                                                                                                                             <  1:426501/150‑1 (MQ=255)
                            tAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTccc                                                                                                                           >  2:1583371/1‑150 (MQ=255)
                               aCGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCAcc                                                                                                                        <  1:1454339/150‑1 (MQ=255)
                                    tACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATTATCGCAACAATTCTCCACTCCCACCGaccc                                                                                                                   >  1:693650/1‑148 (MQ=255)
                                     aCCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACAcc                                                                                                                  <  1:106773/150‑1 (MQ=255)
                                        ctctGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCCGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACAccccc                                                                                                               >  1:24400/1‑150 (MQ=255)
                                         tctGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCt                                                                                                              >  1:1334643/1‑150 (MQ=255)
                                          ctGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCCCCGACACCCCCtg                                                                                                             >  1:1166781/1‑150 (MQ=255)
                                            gAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCCGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtgat                                                                                                           >  1:1509756/1‑150 (MQ=255)
                                             aaTGCTCTCCTTTAGTCGAACAAAACCTGATGCCTTGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtgatg                                                                                                          <  1:278998/150‑1 (MQ=255)
                                             aaTGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtgatg                                                                                                          <  1:1666999/150‑1 (MQ=255)
                                             aaTGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtgatg                                                                                                          <  1:1866997/150‑1 (MQ=255)
                                             aaTGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtgatg                                                                                                          <  1:207427/150‑1 (MQ=255)
                                                 ctctCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtgatgatga                                                                                                      <  2:138868/150‑1 (MQ=255)
                                                  tctcCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGaa                                                                                                     <  1:825139/150‑1 (MQ=255)
                                                  tctcCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGaa                                                                                                     <  2:692691/150‑1 (MQ=255)
                                                    tcCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATc                                                                                                   <  1:1751814/150‑1 (MQ=255)
                                                      cTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACAACCCCTGATGATGAATCAg                                                                                                 >  2:1858301/1‑150 (MQ=255)
                                                      cTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATAACGGCCTTAGCCGTAGTGGAACTAAACTACTATGATGATCATCATCATCGCAACAATTCTCCACTACCACCGACAACCCCTGATGATGAATCAg                                                                                                 >  1:1477107/1‑150 (MQ=255)
                                                       tttAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAga                                                                                                <  1:1619703/150‑1 (MQ=255)
                                                        ttAGTCGAACAAAACCTGATGCCATGGATTTACACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATATCAACAATTCTCGACTCCCACCGACACCCCCTTATGATGAATCAgac                                                                                               <  2:1293322/150‑1 (MQ=255)
                                                           gTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACACTTCTCCACTCCCACCGACACCCCCTGATGATGAATCAgacgac                                                                                            >  2:1228650/1‑150 (MQ=255)
                                                           gTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAgacgac                                                                                            >  2:1691901/1‑150 (MQ=255)
                                                             cGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGacac                                                                                          >  2:440644/1‑150 (MQ=255)
                                                              gAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACt                                                                                         >  2:1866298/1‑150 (MQ=255)
                                                               aaCAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGCCGATGATGATCATCATCATCGCAACAATTCTCCACTCCCCCCGACACCCCCTGATGATGAAGCACACGACACTc                                                                                        >  1:1856174/1‑150 (MQ=255)
                                                               aaCAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCCCTCCCACCGACCCCCCCTGATGATGAATAAGACGacacac                                                                                        >  1:415535/1‑148 (MQ=255)
                                                               aaCAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtgatgatga                                                                                                      >  1:1070143/1‑136 (MQ=255)
                                                               aaCAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtgatgatga                                                                                                      <  2:1070143/136‑1 (MQ=255)
                                                                       cTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCCGTTccc                                                                                >  1:1378659/1‑150 (MQ=255)
                                                                          aTGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCa                                                                             <  1:1882788/150‑1 (MQ=255)
                                                                           tGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCGCTGATGATGCATCCGACGACCCTCCAGTTcccccca                                                                            >  1:309757/1‑148 (MQ=255)
                                                                            gCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAAc                                                                           >  1:1801696/1‑150 (MQ=255)
                                                                                     ttCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCgg                                                                  >  1:1746550/1‑150 (MQ=255)
                                                                                      tCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGga                                                                 <  2:1327418/150‑1 (MQ=255)
                                                                                        cacaGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGgaga                                                               <  2:1746550/150‑1 (MQ=255)
                                                                                                    tgcGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCACCCAACTCCTGGCGGAGATTAGATAATAcc                                                   >  2:22131/1‑150 (MQ=255)
                                                                                                     gcGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCg                                                  >  2:1609660/1‑150 (MQ=255)
                                                                                                      cGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCgg                                                 >  1:660444/1‑150 (MQ=255)
                                                                                                       gTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCggc                                                >  1:1839928/1‑149 (MQ=255)
                                                                                                        ttAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGAc                                               <  1:1494302/150‑1 (MQ=255)
                                                                                                         tAAACAGGGCCTTAGCGTTAGTGGAACTAAACGACGATGATGATCATAATCATCGCAACAATTCTCCCCTCCCCCCGACCCCCCCTGATGATTAATCAGACGACACTCCAGTTCCCCCAACTCCTCGCGGAGATGAGATAATACCGGACg                                              <  2:947586/150‑1 (MQ=255)
                                                                                                            tCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATc                                           >  2:1063170/1‑150 (MQ=255)
                                                                                                             cACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATcc                                          >  1:1794909/1‑150 (MQ=255)
                                                                                                             cACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCAAGTTACCCCAACTCCTGGCGGAGATCAGATAATACCGGCAGATcc                                          >  2:1445787/1‑150 (MQ=255)
                                                                                                             aaCGGCATTCGCCGTAGTGGAACTAAACGACGATTATGATCATATTCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACTACACTCCAGTTCCCCCAACTCCTGGCGGATATGAGATAATACCGGACGATcc                                          <  2:1477107/149‑1 (MQ=255)
                                                                                                                   cTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATTAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGgatga                                    <  1:1705098/150‑1 (MQ=255)
                                                                                                                    ttAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGgatgat                                   <  2:807050/150‑1 (MQ=255)
                                                                                                                       gCCGTAGTGGAACTAAACGACGATTATTATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACg                                <  1:1775565/150‑1 (MQ=255)
                                                                                                                       gCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACg                                <  1:514697/150‑1 (MQ=255)
                                                                                                                         cGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGcc                              <  2:435898/150‑1 (MQ=255)
                                                                                                                         cGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGcc                              <  2:1284303/150‑1 (MQ=255)
                                                                                                                           tAGTGGAACTAAACGACGATGATGAACATCATCATCGCAACAATTCTCCACTCCCACCGACACCGCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTa                            <  1:135317/150‑1 (MQ=255)
                                                                                                                            aGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGCGACCCCCTGATGATGAATCAGACGACACTCCAGTTGCCCCAACTCCTGGCGGTGATGAGATAATACCGGACGTTCCGGATGATCCGCCGac                           >  1:1174539/1‑150 (MQ=255)
                                                                                                                             gTGGAACTAAACGACGTTGTTGATATTAATCATCGCAACAATTCTCCACTCCAACCGACACCCCGTGATGTTGAATCAGACGACTCTCCTGTTCCGCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTact                          <  2:724396/150‑2 (MQ=255)
                                                                                                                             gTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTacac                         <  1:1056729/150‑1 (MQ=255)
                                                                                                                                gAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTGCAGTCCCACCGACACCCCCTGATGATGACTCCGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACTATCCTGATGATACGCCTACACCt                       <  2:531207/150‑1 (MQ=255)
                                                                                                                                    tAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCa                   <  2:1509756/150‑1 (MQ=255)
                                                                                                                                          gcgaTGATGATGATCATCATCGCAACAATTCTCCACTCCCCCCGACCCCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCgg             >  1:1060409/2‑150 (MQ=255)
                                                                                                                                          gcgaTGATGATGATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCGGATGATGAATAAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCgg             >  1:1137425/2‑150 (MQ=255)
                                                                                                                                          ccgaTGATGATCATCATCATCGAAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCgg             <  2:1899786/149‑1 (MQ=255)
                                                                                                                                           cgaTGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGt            >  2:817166/1‑150 (MQ=255)
                                                                                                                                             atgatgatCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGttt          <  2:1334643/150‑1 (MQ=255)
                                                                                                                                             atgatgatCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGttt          >  2:638268/1‑150 (MQ=255)
                                                                                                                                             atgatgatCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGttt          >  1:587603/1‑150 (MQ=255)
                                                                                                                                              tgatgatCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTAGCCCTCCCAACGCGGttgc         >  2:305432/1‑148 (MQ=255)
                                                                                                                                              tgatgatCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTc         >  2:122578/1‑150 (MQ=255)
                                                                                                                                               gatgatCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCg        <  1:1691901/150‑1 (MQ=255)
                                                                                                                                                atgatCATCATCATCGAAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGt       <  1:1583371/150‑1 (MQ=255)
                                                                                                                                                 tgatCATCATCATCGCATCAATTCTCCACTCCCCCCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGtt      >  1:1088927/1‑150 (MQ=255)
                                                                                                                                                 tgatCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGtt      >  1:829622/1‑150 (MQ=255)
                                                                                                                                                  tatcatcataatcGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGttt     <  1:1845039/149‑1 (MQ=255)
                                                                                                                                                    tcatcatcatcGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGCCGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTtaa   >  1:656872/1‑150 (MQ=255)
                                                                                                                                                     catcatcatcGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTtaat  <  1:1641407/150‑1 (MQ=255)
                                                                                                                                                     catcatcatcGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTtaat  <  1:1609660/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAAT  >  NC_000913/1465522‑1465820

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: