Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,205,497 T→C pseudogene (1094/1134 nt) aaaE → pseudogene, phage terminase protein A family, e14 prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,205,4970TC98.3% 389.5 / NA 116pseudogene (1094/1134 nt)aaaEpseudogene, phage terminase protein A family, e14 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (42/72);  minor base A (2/0);  total (44/72)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 116 total aligned reads displayed.

ATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGATTTATCCCCCGGCGGGTTTTGTT  >  NC_000913/1205352‑1205644
                                                                                                                                                  |                                                                                                                                                    
aTCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTg                                                                                                                                                 <  3:529399/151‑1 (MQ=255)
   ggCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATc                                                                                                                                              <  1:1323677/151‑1 (MQ=255)
      aaCGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATcaca                                                                                                                                           <  2:93509/151‑1 (MQ=255)
         gTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTg                                                                                                                                        <  4:782631/151‑1 (MQ=255)
          tGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGGCCAGCTGTACGAGAAAAAAGACACGCCGCCCGATcacact                                                                                                                                         >  4:336441/1‑147 (MQ=255)
           ggTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGa                                                                                                                                       <  1:1137130/150‑1 (MQ=255)
           ggTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGa                                                                                                                                       <  4:1660073/150‑1 (MQ=255)
                gCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTa                                                                                                                                 <  1:522309/151‑1 (MQ=255)
                 cAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTAc                                                                                                                                <  3:297076/151‑1 (MQ=255)
                  aaaaCCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACg                                                                                                                               <  4:402451/151‑1 (MQ=255)
                  aaaaCCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACg                                                                                                                               <  4:1438381/151‑1 (MQ=255)
                  aaaaCCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACg                                                                                                                               <  1:1533617/151‑1 (MQ=255)
                  aaaaCCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGGTGGGGGAGGCACGCTGTTCGAGGAAAAAAACAAGCCCCCTGATtacaatttgtcc                                                                                                                                   >  1:501877/1‑139 (MQ=255)
                   aaaCCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGTGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACgg                                                                                                                              <  4:207232/151‑1 (MQ=255)
                   aaaCCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGGCACGCGGTGCGAGAAAGAAGAAACACCGCCTGTTCACAATTAGTCCTACg                                                                                                                               >  2:507259/1‑150 (MQ=255)
                   aaaCCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACgg                                                                                                                              <  2:104083/151‑1 (MQ=255)
                    aaCCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACggg                                                                                                                             <  4:536689/151‑1 (MQ=255)
                      ccATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCTGAGCCCAGCTGTGCAGGAAAAAAGAAACGCCGCCTGTTCACATTTGGTCCttcagg                                                                                                                             >  4:1048572/1‑144 (MQ=255)
                      ccATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGAt                                                                                                                           <  1:292987/151‑1 (MQ=255)
                       cATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTTGGAGAGGCATGCTGTACGAGAAAAAAAAAACGCCGCCTGATCACATTTAGTTCTACggg                                                                                                                             >  4:1208416/1‑148 (MQ=255)
                       cATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTTGCGGGAGCAAGCTGTTCGAGAAAGAAGAAAAGCCGGCTGATTACATTTAGTCCCTCggg                                                                                                                             >  4:1371217/1‑148 (MQ=255)
                        aTTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGGAGAGGCAAGCTGTTCGAGAAAAAAAAAACACCCCCTGATCACATTTAGTCCTTCCGGGTcc                                                                                                                         >  4:1261955/1‑151 (MQ=255)
                          tCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCcg                                                                                                                        <  4:576888/150‑1 (MQ=255)
                            cGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCcgctc                                                                                                                     <  3:65831/151‑1 (MQ=255)
                            cGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCcgctc                                                                                                                     <  1:1392571/151‑1 (MQ=255)
                                 aaCGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCttt                                                                                                               <  1:698338/151‑1 (MQ=255)
                                 aaCGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCAGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCtt                                                                                                                <  3:249461/150‑1 (MQ=255)
                                     atgatgTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGTAGGGCCAAGGTGGGAGAGGAAAAAGAAAAGAAGACCGGTCACATTAAGGTCTACGGGGTCCGCGCGCCTTTaat                                                                                                            >  1:1507858/1‑150 (MQ=255)
                                     atgatgTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGGAGAGGCCTCCCGTTCGAGAAAAAAGACACGCCGCCTGATCTACTTGAGGGCTACCGGGTCCGCCCGCttt                                                                                                               >  1:234761/1‑147 (MQ=255)
                                          gtggtgAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAgg                                                                                                       <  4:1535304/150‑1 (MQ=255)
                                          gtggtgAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGt                                                                                                      <  1:1718973/151‑1 (MQ=255)
                                            ggtgAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTaa                                                                                                    <  2:200233/151‑1 (MQ=255)
                                               gAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTacgttatgtcaata                                                                                                 <  1:1596642/151‑14 (MQ=255)
                                               gAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAAtt                                                                                                  <  3:338932/150‑1 (MQ=255)
                                               gAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAAtt                                                                                                  <  3:486622/150‑1 (MQ=255)
                                               gAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTa                                                                                                 <  2:476272/151‑1 (MQ=255)
                                               gAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCAAGCTGGACGAGAAAAAAGAAACGCCGCCTGATCACATTGAGTCCTACGGGGGCCGCCCGCCTTTACTGAgg                                                                                                       >  2:149315/1‑145 (MQ=255)
                                                 aGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATg                                                                                               <  2:151704/151‑1 (MQ=255)
                                                   cccGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGAt                                                                                             <  1:1758220/151‑1 (MQ=255)
                                                     cGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCa                                                                                           <  4:1054322/151‑1 (MQ=255)
                                                     cGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGACATCCGGTACGAGGAAAAAGAAACGCCGCCTGAACACAATTAGTGCTTCCGGGGCCGCCGGCCTTTACTTAGGGAATTATg                                                                                               >  4:889771/1‑147 (MQ=255)
                                                         aaGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCt                                                                                       <  1:965577/151‑1 (MQ=255)
                                                         aaGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCt                                                                                       <  1:1357487/151‑1 (MQ=255)
                                                           ggAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTg                                                                                      <  1:44885/150‑1 (MQ=255)
                                                            gAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGGAGAGGCCTGCTGTTCGCGAAAAAAAACACGCCCGCTGTTCACATTTAGTCCTACCGGGAGCGCCCACGTTTAATGAGGGAAATATTATCAttgctct                                                                                    >  2:197997/1‑145 (MQ=255)
                                                             aGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGGAGAGCCATGCTGTAAGAGAAAAAAAAAACACCCGCTGATCAAATTTAGTTCTACGGGGTCCGCCCCCCTTAAATTTGGGAATTATGATAATTCTGAtt                                                                                   >  4:481525/1‑151 (MQ=255)
                                                               cAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGGAGAGGCATGCTGTACGAGAAAAAAAAAACGCCCGCTGATCACAATTAGTTCTACGGGGTCCGCCCGCCTTTACTTAGGGAAATATGATCATGCTGAtttt                                                                                 >  4:1190637/1‑149 (MQ=255)
                                                                  gCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGGAGAGGCCTGCTGGTCGAGGAAAAAGAAACGCCGTCTGAGCACATTAGGTCCCACGGGGTCCGCCCGCCTTAACTGAGGTAAATATTATctttgtgagt                                                                                   >  4:1552081/1‑137 (MQ=255)
                                                                  gCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTcgc                                                                              <  4:1746825/151‑1 (MQ=255)
                                                                    ggAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGGAGAGCCCTTCTGTACGAGAAAAAAAAAACCCCCTCTGCTTACATTTTGGTCTTCCGGGGCCCCCCGCCTTTTCTGAGGGGATTATGATCATACTGGtt                                                                                   >  4:1473116/1‑144 (MQ=255)
                                                                     gAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGGCGGGCCATGCCGTAGGAGGAAAAAGAAACACCGCCTGATTACATTTAGTCCCTCGGGGTCCGCCCGCCTTAACTGAGGGGATTATTATTTTGTTGATTTTTGCGcc                                                                           >  4:1900423/1‑151 (MQ=255)
                                                                      aaaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGGAGAGGCATGCTGGAAGGGGAAAAAGAAAACCCCGCTGGTCACATTTAGTCCTACGGGGGCCGCCCGCCTTTACTTAGGGGAATATTATTATGATGATTTTCGCGcc                                                                           >  1:1645289/1‑150 (MQ=255)
                                                                      aaaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGGAGAGACATGCCGGAAGAGAAAGAAGACAAGCCGGCCGTTCACATTTAGTCTTACGGGGAGCGCCCGCCTTTACTCAGGGGAATATTAAGATGCTGTTTTTCCCGcc                                                                           >  1:770003/1‑150 (MQ=255)
                                                                       aaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCctc                                                                         <  1:104924/151‑1 (MQ=255)
                                                                         aCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGGCGAGGCACGCTGTTCGCGAAAGAAGAAACACCGCCCGATCACATTTAGTCCTACGGGGTCCGCCCGCTTTTACTGAGGGAATTATGATCATGTTGATtctccc                                                                              >  1:514189/1‑142 (MQ=255)
                                                                         aCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCtg                                                                       <  1:1263994/151‑1 (MQ=255)
                                                                           aaaaTCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCCGAGCCATGCCGTACCAGAAAGAAGAAAAGCCGTCTGATCACATTTAGTCCTACCGGGCCCGCTCGCCTTTTATGAGGTAATTATTATCATGCTGATTTTTGCGcccc                                                                         >  4:1270180/1‑145 (MQ=255)
                                                                           aaaaTCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCtggt                                                                     <  1:1858094/151‑1 (MQ=255)
                                                                            aaaTCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGGGGGGCCCAGGTGGTCGAGGAAAAAAAAAAACCGCCTGGTCACACTGAGTTCTCCGGGGGCCCCTCGCGTTTACTTAGGGAAATATGATCATTGTGATTTTTgcgc                                                                            >  1:749241/1‑143 (MQ=255)
                                                                            aaaTCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCtggtg                                                                    <  3:507870/151‑1 (MQ=255)
                                                                              aTCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTggg                                                                  <  1:1684982/151‑1 (MQ=255)
                                                                                cGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGc                                                                 <  4:921053/150‑1 (MQ=255)
                                                                                        cAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGt                                                        <  1:6731/151‑1 (MQ=255)
                                                                                          gTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTgc                                                      <  4:578714/151‑1 (MQ=255)
                                                                                            tGCGCTGATTATGGCGGTTGGCGGGGCCATGCCGTTCGAGAAAAAAGAAACCCCCGCCGGTCCCATTGAGTTCTACGGGGGCCGCCCGCCTTTAATTAGGGAAATATTATTATTCTGGTTTTCCCGCCCCTGGTGGgggggcggggggcg                                                     >  1:1354262/1‑150 (MQ=255)
                                                                                              cgcTGATTATGGCGGTTGGCAATGACAACCTCTACAGTGAAAAAGACACCACGCCTGAACACATTAAGTTCTACAGGATACGCCCCCTTTTAATTAGGGGATTAATAACATGCTGATTCTCGCGCCCCTCGTGGGGGGGCTggggg                                                       >  1:662628/1‑144 (MQ=255)
                                                                                               gcTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCttt                                                 <  1:218106/151‑1 (MQ=255)
                                                                                                  gATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCtttt                                                <  1:848032/149‑1 (MQ=255)
                                                                                                  gATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCtttt                                                <  3:272522/149‑1 (MQ=255)
                                                                                                  gATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGc                                              <  3:540243/151‑1 (MQ=255)
                                                                                                      aTGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTgg                                           <  4:742570/150‑1 (MQ=255)
                                                                                                      aTGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGc                                          <  4:1114436/151‑1 (MQ=255)
                                                                                                        ggCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGt                                        >  1:475588/1‑151 (MQ=255)
                                                                                                          cGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTAt                                      <  4:21148/151‑1 (MQ=255)
                                                                                                           ggTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATg                                     <  1:1231546/151‑1 (MQ=255)
                                                                                                           ggTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATg                                     >  2:496092/1‑151 (MQ=255)
                                                                                                           ggTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATg                                     <  4:325014/151‑1 (MQ=255)
                                                                                                             ttGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGt                                   >  1:249610/1‑151 (MQ=255)
                                                                                                               ggCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGc                                 <  3:26356/151‑1 (MQ=255)
                                                                                                                gCAGAGCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGcc                                >  4:467989/1‑151 (MQ=255)
                                                                                                                  agagCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTg                              <  4:1919384/151‑1 (MQ=255)
                                                                                                                   gagCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTgg                             <  3:540740/151‑1 (MQ=255)
                                                                                                                    agCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGc                            <  3:257499/151‑1 (MQ=255)
                                                                                                                    agCCATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGc                            >  4:790162/1‑151 (MQ=255)
                                                                                                                      ccATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTg                          >  3:103722/1‑151 (MQ=255)
                                                                                                                      ccATGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTg                          >  1:1216820/1‑151 (MQ=255)
                                                                                                                        aTGCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGAt                        <  3:322306/151‑1 (MQ=255)
                                                                                                                          gCTGTACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGAttt                      >  3:97239/1‑151 (MQ=255)
                                                                                                                              tACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGATTTATcc                  >  4:1772138/1‑151 (MQ=255)
                                                                                                                              tACGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGATTTATcc                  >  4:1174809/1‑151 (MQ=255)
                                                                                                                               aCGAGAAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGATTTATcc                  <  4:1489075/150‑1 (MQ=255)
                                                                                                                                  gagaAAGAAGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGATTTATccccc               <  2:500651/151‑1 (MQ=255)
                                                                                                                                     aaagaaGACACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGATTTATCCCCcggc            >  4:1231619/1‑151 (MQ=255)
                                                                                                                                     aaagaaGACACGCCGTCTGAACACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTTCCTGGCTGATTTATCCCCcggc            <  1:1373972/151‑1 (MQ=255)
                                                                                                                                       agaagaCACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGATTTATCCCCcggcgg          <  4:794525/151‑1 (MQ=255)
                                                                                                                                         aagaCACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGATTTATCCCCCGGCGGGt        >  3:176676/1‑151 (MQ=255)
                                                                                                                                          agaCACGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGATTTATCCCCCGGCGGGtt       >  1:1159531/1‑151 (MQ=255)
                                                                                                                                             cacGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGATTTATCCCCCGGCGGGtttt     >  1:1716686/1‑150 (MQ=255)
                                                                                                                                              acGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGATTTATCCCCCGGCGGGTTttgt   >  1:1587108/1‑151 (MQ=255)
                                                                                                                                               cGCCGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGATTTATCCCCCGGCGGGTTttgtt  >  1:957871/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCTGATTCTCGCGCCTCTGGTGGGCGTGCTGGGTGCGCTTTTGCTGGCGTATGGTGCCTGGCTGATTTATCCCCCGGCGGGTTTTGTT  >  NC_000913/1205352‑1205644

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