Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,465,671 C→G 5.1% 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.1% 335.2 / 8.0 98pseudogene (280/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (51/42);  new base G (2/3);  total (53/45)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

Only 100 of 109 total aligned reads displayed.

CCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAA  >  NC_000913/1465523‑1465819
                                                                                                                                                    |                                                                                                                                                    
cccAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatc                                                                                                                                                      <  2:268279/149‑1 (MQ=255)
 ccAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatca                                                                                                                                                     <  8:235060/149‑1 (MQ=255)
 ccAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatca                                                                                                                                                     >  5:301930/1‑149 (MQ=255)
  cAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcat                                                                                                                                                    <  2:245768/149‑1 (MQ=255)
           tCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCg                                                                                                                                            >  3:141745/1‑148 (MQ=255)
           tCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCg                                                                                                                                            >  2:282948/1‑148 (MQ=255)
             aGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATGATCATCATCgtc                                                                                                                                          <  8:467859/148‑3 (MQ=255)
                 cGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAAt                                                                                                                                     >  4:285346/1‑149 (MQ=255)
                   gAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATtc                                                                                                                                   <  4:278797/149‑1 (MQ=255)
                     agagCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATtctc                                                                                                                                 <  1:142570/149‑1 (MQ=255)
                     agagCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATtctc                                                                                                                                 <  8:406499/149‑1 (MQ=255)
                      gagCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATtctc                                                                                                                                 >  7:291252/1‑148 (MQ=255)
                        gCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCAc                                                                                                                              <  5:427037/149‑1 (MQ=255)
                          ttAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTc                                                                                                                            >  8:440924/1‑149 (MQ=255)
                             aaCGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCa                                                                                                                         >  5:27506/1‑149 (MQ=255)
                                  ttACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCgcc                                                                                                                    >  2:467181/1‑147 (MQ=255)
                                     cccTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGacac                                                                                                                  <  6:301930/148‑1 (MQ=255)
                                       ctctGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACAcccc                                                                                                               >  3:326486/1‑149 (MQ=255)
                                       ctctGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACAcccc                                                                                                               >  2:255289/1‑149 (MQ=255)
                                       ctctGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACAcccc                                                                                                               <  8:291252/149‑1 (MQ=255)
                                       ctctGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACAcccc                                                                                                               >  2:238599/1‑149 (MQ=255)
                                       ctctGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATAATCGCAACAATTCTCCACTCCCACCGACAcccc                                                                                                               >  5:201096/1‑149 (MQ=255)
                                        tctGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGATCTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACAccccc                                                                                                              >  1:343967/1‑149 (MQ=255)
                                        tctGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGcaa                                                                                                                                         <  7:349841/122‑1 (MQ=255)
                                        tctGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGcaa                                                                                                                                         >  8:349841/1‑122 (MQ=255)
                                            aaTGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtgat                                                                                                          >  8:114815/1‑149 (MQ=255)
                                            aaTGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtgat                                                                                                          >  5:288813/1‑149 (MQ=255)
                                            aaTGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtga                                                                                                           <  4:326486/148‑1 (MQ=255)
                                            aaTGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtg                                                                                                            <  1:429906/147‑1 (MQ=255)
                                                 tctcCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTGTCCACTCCCACCGACCCGCCCtgaagatg                                                                                                      >  6:292582/1‑148 (MQ=255)
                                                 tctcCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtgatgatg                                                                                                      >  7:157276/1‑148 (MQ=255)
                                                 tctcCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtgatgatg                                                                                                      >  8:446649/1‑148 (MQ=255)
                                                   tcCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtgatgatga                                                                                                     <  1:2068/147‑1 (MQ=255)
                                                     cTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCa                                                                                                 <  7:446649/149‑1 (MQ=255)
                                                     cTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCa                                                                                                 <  8:157276/149‑1 (MQ=255)
                                                      tttAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGCCACCCCCTGATGATGAAt                                                                                                   >  8:228737/1‑146 (MQ=255)
                                                      tttAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAAt                                                                                                   >  7:210089/1‑146 (MQ=255)
                                                      tttAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAAACACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACACCCTGATGATGAAt                                                                                                   >  8:167251/1‑146 (MQ=255)
                                                        tAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAga                                                                                               <  7:231724/148‑1 (MQ=255)
                                                          gTCGAACAAAACCTGCTGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAgacga                                                                                            <  6:288813/149‑1 (MQ=255)
                                                           tCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGGCGATGATGATGATCATCATCgtc                                                                                                                                          >  1:154407/1‑100 (MQ=255)
                                                           tCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGGCGATGATGATGATCATCATCgtc                                                                                                                                          <  2:154407/102‑3 (MQ=255)
                                                           tCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCCCCGACCCCCCCTGAGGATGAATCAgacgc                                                                                            >  4:1981/1‑147 (MQ=255)
                                                           tCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAgacgac                                                                                           >  4:415880/1‑149 (MQ=255)
                                                             gAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcatcatcatc                                                                                                                                             >  1:248761/1‑97 (MQ=255)
                                                             gAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGatcatcatcatc                                                                                                                                             <  2:248761/97‑1 (MQ=255)
                                                                 aaaaCCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCa                                                                                     <  6:27506/149‑1 (MQ=255)
                                                                      cTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTcc                                                                                <  1:255289/149‑1 (MQ=255)
                                                                      cTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTcc                                                                                <  1:336289/149‑1 (MQ=255)
                                                                      cTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTcc                                                                                <  6:34826/149‑1 (MQ=255)
                                                                       tGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACt                                                                                                                             >  8:36468/1‑103 (MQ=255)
                                                                         aTGCCATGGCTTTCCACAGGCGCTGCTTCGTTAATCACGGCCTTAGCCGTAGTGTAACTAAACGACGATGATGATCATCATCATCGCAACATTTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCATTTccccc                                                                             <  3:140914/149‑1 (MQ=255)
                                                                            ccATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAAc                                                                          >  8:437172/1‑149 (MQ=255)
                                                                             cATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTcc                                                                                <  8:199107/142‑1 (MQ=255)
                                                                             cATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTcc                                                                                >  7:199107/1‑142 (MQ=255)
                                                                              aTGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTc                                                                        <  3:100180/149‑1 (MQ=255)
                                                                              aTGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTc                                                                        >  2:91669/1‑149 (MQ=255)
                                                                               tGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTcc                                                                       >  5:399552/1‑149 (MQ=255)
                                                                                ggCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCt                                                                      >  5:411861/1‑149 (MQ=255)
                                                                                ggCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCt                                                                      >  7:465846/1‑149 (MQ=255)
                                                                                 gCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTg                                                                     <  1:238599/149‑1 (MQ=255)
                                                                                       cacaGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCg                                                                  >  6:239974/1‑146 (MQ=255)
                                                                                        acaGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGgaga                                                              <  3:285346/149‑1 (MQ=255)
                                                                                         caGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGgagat                                                             <  3:347208/149‑1 (MQ=255)
                                                                                          aGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTc                                                                                       <  5:45359/122‑1 (MQ=255)
                                                                                          aGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGgagatg                                                            <  7:114815/149‑1 (MQ=255)
                                                                                            gcgcTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGgagatgag                                                          <  3:118359/149‑1 (MQ=255)
                                                                                              gctgctGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGgagatgagat                                                        <  3:415880/149‑1 (MQ=255)
                                                                                                 gctgcGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGataa                                                      >  5:276870/1‑148 (MQ=255)
                                                                                                  ctgcGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGataa                                                      <  7:437796/147‑1 (MQ=255)
                                                                                                     tgTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGga                                                                <  4:87613/133‑1 (MQ=255)
                                                                                                     cGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATAc                                                   <  4:141745/147‑1 (MQ=255)
                                                                                                     cGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCg                                                 >  1:438377/1‑149 (MQ=255)
                                                                                                      gTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGaca                                                                                          >  1:102000/1‑107 (MQ=255)
                                                                                                      gTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGaca                                                                                          <  2:102000/107‑1 (MQ=255)
                                                                                                      gTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGaca                                                                                          <  3:338423/107‑1 (MQ=255)
                                                                                                      gTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGaca                                                                                          >  4:338423/1‑107 (MQ=255)
                                                                                                       ttAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGa                                               <  7:440924/149‑1 (MQ=255)
                                                                                                         aaTCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGgagat                                                             >  7:72763/1‑133 (MQ=255)
                                                                                                             aCGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATcc                                         >  7:316452/1‑149 (MQ=255)
                                                                                                             aCGGCCTTAGCCGTAGGGGAACTAAACGCGGATGATTATCATCATCTTCTCAACAATTCTCCCCTCCCACCGACACCCCCTGATGATGATTCAGACGACACTCGAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATc                                          <  2:466616/148‑1 (MQ=255)
                                                                                                              cGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCg                                        <  6:399552/149‑1 (MQ=255)
                                                                                                                 tcatAGCCGTAGTGGAAATAAAAGCCGATGTTGATCATCATCATCGCTACATTTCACCCCTCCCACCGACACCCCCTGATGATGATTCAGACGACACTCCAGTTCCCCCAACTCCTTGCGGAGATTAGATAATACCGGACGCTCCGgat                                     <  2:301196/146‑1 (MQ=255)
                                                                                                                  cTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGgat                                     >  3:49729/1‑148 (MQ=255)
                                                                                                                         gTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGcc                             >  7:336969/1‑149 (MQ=255)
                                                                                                                            gTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCg                                                 >  8:138438/1‑126 (MQ=255)
                                                                                                                            gTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCg                                                 <  7:138438/126‑1 (MQ=255)
                                                                                                                               gAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCCCCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGCTACGCCTccacg                       >  1:307992/1‑148 (MQ=255)
                                                                                                                                  cTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTcc                    <  1:91669/149‑1 (MQ=255)
                                                                                                                                   tAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGCCCCCCCCTGAAGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTCGCGGAGATGAGAAAATACCGGACGATCCGCATGATACGCCTACACCTcc                    >  5:264916/1‑148 (MQ=255)
                                                                                                                                   attgttgcgaTGATGATGATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTcc                    <  8:256015/141‑1 (MQ=255)
                                                                                                                                   attgttgcgaTGATGATGATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTcc                    >  7:256015/8‑148 (MQ=255)
                                                                                                                                        gacgaTGATGATCATCATCATCGCAACAATTCTCCCCTCCCACCGACCCCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAccc              >  5:228846/1‑149 (MQ=255)
                                                                                                                                           gatgatgatCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGt           >  7:93298/1‑149 (MQ=255)
                                                                                                                                            atgatgatCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGtt          >  4:37384/1‑149 (MQ=255)
                                                                                                                                             tgatgatCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGttt         >  3:242433/1‑149 (MQ=255)
                                                                                                                                                 gatgatCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGtt     >  3:185361/1‑149 (MQ=255)
                                                                                                                                                  atcatcatcatcGCAACAATTCTCCACTCCCACCGACACCCCCTGATTATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGtt     <  2:343967/148‑1 (MQ=255)
                                                                                                                                                  atcatcatcatcGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGtt     <  4:49729/148‑1 (MQ=255)
                                                                                                                                                    catcatcatcGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGCCACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTGTCGGTtaa  >  3:355776/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAA  >  NC_000913/1465523‑1465819

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