Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,650,147 G→A 7.0% pseudogene (70/921 nt) ydfE → Qin prophage; protein YdfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,650,1470GA 7.0% 303.6 / 8.8 99pseudogene (70/921 nt)ydfEQin prophage; protein YdfE
Reads supporting (aligned to +/- strand):  ref base G (49/43);  new base A (3/4);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

TAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCG  >  NC_000913/1649999‑1650294
                                                                                                                                                    |                                                                                                                                                   
taaataTCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAATTCTTTATTTGCTCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGCCTGt                                                                                                                                                    >  2:740946/1‑150 (MQ=255)
taaataTCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGt                                                                                                                                                    >  2:596061/1‑150 (MQ=255)
 aaataTCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTa                                                                                                                                                   <  2:591214/150‑1 (MQ=255)
  aataTCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAg                                                                                                                                                  <  1:1164806/150‑1 (MQ=255)
       ccGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCcaca                                                                                                                                             <  1:1290025/150‑1 (MQ=255)
          gTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGcc                                                                                                                                          <  1:1004260/150‑1 (MQ=255)
          gTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGcc                                                                                                                                          <  2:1147499/150‑1 (MQ=255)
          gTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGcc                                                                                                                                          <  1:1283070/150‑1 (MQ=255)
          gTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGATGTGGCTGTAGCCACAGca                                                                                                                                          >  1:1052002/1‑149 (MQ=255)
           tttCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgc                                                                                                                                         <  1:964794/150‑3 (MQ=255)
            ttCATTCTTATATGCCAGCAATGTCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgcc                                                                                                                                        <  2:463650/150‑4 (MQ=255)
               aTTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGa                                                                                                                                     >  2:137947/1‑150 (MQ=255)
               aTTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGa                                                                                                                                     <  1:301816/150‑1 (MQ=255)
               aTTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgccttc                                                                                                                                     <  1:1012544/150‑7 (MQ=255)
                ttCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGaa                                                                                                                                    <  2:314195/150‑1 (MQ=255)
                     atatGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTgg                                                                                                                               <  1:349038/150‑1 (MQ=255)
                     atatGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTgg                                                                                                                               >  1:125307/1‑150 (MQ=255)
                        tGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcga                                                                                                                            >  2:961563/1‑150 (MQ=255)
                        tGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcga                                                                                                                            >  2:621680/1‑150 (MQ=255)
                        tGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcga                                                                                                                            >  1:980994/1‑150 (MQ=255)
                        tGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcga                                                                                                                            >  2:1122657/1‑150 (MQ=255)
                         gCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCACTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcgac                                                                                                                           >  1:235445/1‑150 (MQ=255)
                          ccAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcgacg                                                                                                                          >  2:716658/1‑150 (MQ=255)
                           cAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcgacga                                                                                                                         >  1:913996/1‑150 (MQ=255)
                            aGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGaa                                                                                                                        >  2:368767/1‑150 (MQ=255)
                              cAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACg                                                                                                                      <  2:662088/150‑1 (MQ=255)
                                    cAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGc                                                                                                                >  2:97874/1‑150 (MQ=255)
                                    cAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGc                                                                                                                >  1:842807/1‑150 (MQ=255)
                                      gggATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACAAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCAACAGCCATTGAAGCTGGCGACGAACGCCGTGCTc                                                                                                              >  1:1412308/1‑150 (MQ=255)
                                       ggATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCg                                                                                                             <  1:97874/150‑1 (MQ=255)
                                       ggATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCACTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCg                                                                                                             <  2:235445/150‑1 (MQ=255)
                                        gATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGACGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGa                                                                                                            >  1:1254575/1‑150 (MQ=255)
                                         aTTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAg                                                                                                           <  1:512874/150‑1 (MQ=255)
                                         aTTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAg                                                                                                           <  1:594974/150‑1 (MQ=255)
                                          tttGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGc                                                                                                          <  2:980994/150‑1 (MQ=255)
                                             gTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCaaa                                                                                                       >  2:508704/1‑150 (MQ=255)
                                                 aTCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAAttt                                                                                                   >  1:828944/1‑150 (MQ=255)
                                                   ccTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCa                                                                                                 >  2:400423/1‑150 (MQ=255)
                                                       aaaTCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGATGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTgg                                                                                             >  2:1051354/1‑150 (MQ=255)
                                                       aaaTCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTgg                                                                                             <  2:61389/150‑1 (MQ=255)
                                                           cTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCg                                                                                                                   <  2:1278421/124‑1 (MQ=255)
                                                           cTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCg                                                                                                                   >  1:1278421/1‑124 (MQ=255)
                                                              tCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTcc                                                                                      <  1:508704/150‑1 (MQ=255)
                                                                    ggTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAAc                                                                                <  2:809204/150‑1 (MQ=255)
                                                                    ggTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACCAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAAc                                                                                >  2:740580/1‑150 (MQ=255)
                                                                    ggTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAACCTGGCGACGAACGCCGTCCTCGAGCAAAATTTCACTCGCAATTCCTGGAAc                                                                                >  2:296302/1‑150 (MQ=255)
                                                                     gTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACa                                                                               >  1:896972/1‑150 (MQ=255)
                                                                       tAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACAtt                                                                             >  1:94396/1‑150 (MQ=255)
                                                                         aaaCAAATGAGTAAAGTCTTTATTTGCGCCGCAATTCCTGACGAACTGCCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGACATTGAAGCTGGCGACGAACGCGGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTAt                                                                           >  1:348733/1‑150 (MQ=255)
                                                                            cAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCg                                                                        >  2:597568/1‑150 (MQ=255)
                                                                                 gAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGTAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTg                                                                                              <  2:332350/123‑1 (MQ=255)
                                                                                 gAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGa                                                                                  >  1:1057406/1‑135 (MQ=255)
                                                                                 gAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGa                                                                                  <  2:1057406/135‑1 (MQ=255)
                                                                                 gAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGAAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTg                                                                                              >  1:332350/1‑123 (MQ=255)
                                                                                    tAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCa                                                                <  1:1118161/150‑1 (MQ=255)
                                                                                    tAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCa                                                                <  2:400039/150‑1 (MQ=255)
                                                                                        ttCTTTATTTGCTCCGCTATTCCTGACGAACTGGCAACAATGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGATGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc                                                            <  1:895875/149‑1 (MQ=255)
                                                                                           tttATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgc                                                         <  2:913996/150‑1 (MQ=255)
                                                                                            ttATTTGCGCCGCTATTCCTGACGAACTGGCAAGAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgcg                                                        <  2:79582/150‑1 (MQ=255)
                                                                                                 tGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTtat                                                   >  1:155439/1‑150 (MQ=255)
                                                                                                   cgccgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATaa                                                 <  2:183123/150‑1 (MQ=255)
                                                                                                   cgccgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATaa                                                 >  1:414801/1‑150 (MQ=255)
                                                                                                     ccgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAAt                                               <  1:528063/150‑1 (MQ=255)
                                                                                                     ccgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAAt                                               >  2:241237/1‑150 (MQ=255)
                                                                                                      gggtATTCCTGAAGATCTGGCTACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTTGCAATTCCTGGAACATTATCCGTCTGCTCAGTACTGCGCTTATAAAtt                                              <  1:740946/147‑1 (MQ=255)
                                                                                                      cgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAAtt                                              <  1:596061/150‑1 (MQ=255)
                                                                                                       gcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAAttt                                             >  1:816188/1‑150 (MQ=255)
                                                                                                           ttCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTg                                         >  2:1116362/1‑150 (MQ=255)
                                                                                                             ccTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATaa                                                 >  1:704367/1‑140 (MQ=255)
                                                                                                             ccTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATaa                                                 <  2:704367/140‑1 (MQ=255)
                                                                                                              cTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCt                                      >  2:472383/1‑150 (MQ=255)
                                                                                                               tGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTg                                     >  1:905512/1‑150 (MQ=255)
                                                                                                                gACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGc                                    <  1:472383/150‑1 (MQ=255)
                                                                                                                gACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGc                                    >  2:325623/1‑150 (MQ=255)
                                                                                                                   aaaCTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAg                                 <  2:1412308/149‑1 (MQ=255)
                                                                                                                    aaCTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAgg                                <  2:1395732/150‑1 (MQ=255)
                                                                                                                     aCTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGa                               <  1:961563/150‑1 (MQ=255)
                                                                                                                        ggCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATaa                            >  1:1151031/1‑150 (MQ=255)
                                                                                                                           aacaaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAAcc                         <  1:137947/150‑1 (MQ=255)
                                                                                                                             caaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAAATTCACTGGCAATTCCTGGAACAATATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCAGCGAGCATAAAccag                       >  2:1157338/1‑148 (MQ=255)
                                                                                                                                gggAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGgctgc                                                                   <  2:75988/103‑1 (MQ=255)
                                                                                                                                gggAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGgctgc                                                                   >  1:75988/1‑103 (MQ=255)
                                                                                                                                 ggAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCAt                   >  1:68883/1‑150 (MQ=255)
                                                                                                                                 ggAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTAAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCAt                   >  1:63679/1‑150 (MQ=255)
                                                                                                                                    agaagGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATAcc                <  2:896972/150‑1 (MQ=255)
                                                                                                                                          ggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgccc          >  2:568254/2‑150 (MQ=255)
                                                                                                                                          cgcTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgccc          <  1:716658/150‑1 (MQ=255)
                                                                                                                                          cgcTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgccc          >  2:117749/1‑150 (MQ=255)
                                                                                                                                           gcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTa                                            <  2:275195/115‑1 (MQ=255)
                                                                                                                                           gcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTa                                            >  1:275195/1‑115 (MQ=255)
                                                                                                                                           gcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgccct         >  1:1176227/1‑150 (MQ=255)
                                                                                                                                              gtgGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgccctgcc      <  1:491333/150‑1 (MQ=255)
                                                                                                                                                 gCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTc   >  1:1293890/1‑150 (MQ=255)
                                                                                                                                                 gCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTc   >  2:45770/1‑150 (MQ=255)
                                                                                                                                                  cTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCg  <  1:1051354/150‑1 (MQ=255)
                                                                                                                                                   tGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAgg                                                              >  1:256944/1‑89 (MQ=255)
                                                                                                                                                   tGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAgg                                                              <  2:256944/89‑1 (MQ=255)
                                                                                                                                                   tGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgccctg        >  1:251459/1‑143 (MQ=255)
                                                                                                                                                   tGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgccctg        <  2:251459/143‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TAAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCG  >  NC_000913/1649999‑1650294

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