Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 579,540 T→C 100% N19D (AAC→GAC)  ybcV ← DLP12 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913579,5400TC98.6% 14051.0 / 53.2 4910N19D (AAC→GAC) ybcVDLP12 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (33/42);  new base C (2384/2427);  total (2433/2477)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-02
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 4912 total aligned reads displayed.

GAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGA  >  NC_000913/579392‑579688
                                                                                                                                                    |                                                                                                                                                    
      ttAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATc                                                                                                                                                >  2:67375/1‑149 (MQ=255)
        aaCATTTACAACCTTTTTAAGTCCTTATATTAACACGGTGTTATCGTTTTCTAACACGAAGTGAATATTATCTGTGGCTAGATATTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCtt                                                                                                                                              >  2:301053/1‑149 (MQ=255)
         aCATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCttt                                                                                                                                             >  2:253237/1‑149 (MQ=255)
         aCATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGGTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCttt                                                                                                                                             <  2:276425/149‑1 (MQ=255)
          cATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCtttt                                                                                                                                            <  3:92673/149‑1 (MQ=255)
           aTTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACTATGTGAATATAATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTc                                                                                                                                           <  6:237670/149‑1 (MQ=255)
               aCAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCGGAATAAAACAATTCACAGTCTAAATCTTTTCGCAc                                                                                                                                       >  5:122316/1‑149 (MQ=255)
                cAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCAc                                                                                                                                       <  2:372815/148‑1 (MQ=255)
                 aaCCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACtt                                                                                                                                     <  8:17497/149‑1 (MQ=255)
                 aaCCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACtt                                                                                                                                     >  6:172811/1‑149 (MQ=255)
                   ccTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGa                                                                                                                                   >  6:147524/1‑149 (MQ=255)
                   ccTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGa                                                                                                                                   <  5:241852/149‑1 (MQ=255)
                     tttttAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATc                                                                                                                                 >  5:15716/1‑149 (MQ=255)
                       tttAATTCCTTTTATTAACAATGTTTTATCGTTTTCTAACACGATTTGAATATTATATTTGTCTATATAGTAAATATAATGTGATACTTTGTGACGTTTTAGTTCAGAATAAAACAATTCACATTCTAAATCTTTTCGCACTTGATCg                                                                                                                                <  5:215611/148‑1 (MQ=255)
                       tttAAGTACTTTTATTAACAAGGCTTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATc                                                                                                                                 <  3:1025/147‑1 (MQ=255)
                         tAAGTCCTTTTATTAAGACGGGGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGaa                                                                                                                              <  4:121991/148‑1 (MQ=255)
                         tAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGaa                                                                                                                              <  3:156812/148‑1 (MQ=255)
                          aaGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGAt                                                                                                                                  <  5:37084/143‑1 (MQ=255)
                          aaGTCCTTTTATTAACACGGTGGTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAAAGTGAGACGTTGTGCCGTTTTAGTTCTGATGAAAACAATTCACAGTCTAAATCTTTTCGAACTAGATCGAata                                                                                                                            >  6:338990/1‑149 (MQ=255)
                             tCCTTTTATTATCACGGTGTTATCGTTTTCTTACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATAttt                                                                                                                         <  7:348087/149‑1 (MQ=255)
                             tCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAGTCTTTTCGCACTTGATCGAATAttt                                                                                                                         <  4:142806/149‑1 (MQ=255)
                               cTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTTAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTc                                                                                                                        <  6:44328/148‑1 (MQ=255)
                                ttttATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACTTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATATTTTCGCACTTGATCGAATATTTCtt                                                                                                                      >  7:373314/1‑149 (MQ=255)
                                ttttATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAATTCTTTTCGCACTTGATCGAAAATTTCtt                                                                                                                      >  4:87881/1‑149 (MQ=255)
                                ttttATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCtt                                                                                                                      >  5:83658/1‑149 (MQ=255)
                                ttttATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCtt                                                                                                                      >  6:304738/1‑149 (MQ=255)
                                 tttATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTCTTCGCACTTGATCGAATAtt                                                                                                                          <  8:127717/144‑1 (MQ=255)
                                       cacacGGTGTCATCGTTTTCAAACACGATGTCAATATTATCTGTGGCTCGATAGTAAATATAATGTGAGACGTTGTGACGTATTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAAt                                                                                                               <  2:324157/148‑1 (MQ=255)
                                       aaCACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGAt                                                                                                                                  <  3:317020/130‑1 (MQ=255)
                                        acacGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCATAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATg                                                                                                              <  3:6670/149‑1 (MQ=255)
                                        acacGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATg                                                                                                              <  1:243085/149‑1 (MQ=255)
                                        acacGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATg                                                                                                              <  7:154750/149‑1 (MQ=255)
                                         cacGGTGTTATAGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATgg                                                                                                             <  4:194167/149‑1 (MQ=255)
                                            ggTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCttt                                                                                                                                             >  7:13000/1‑114 (MQ=255)
                                             gtgtTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAAc                                                                                                         <  8:235422/149‑1 (MQ=255)
                                                ttATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAAGCt                                                                                                       >  4:255545/1‑148 (MQ=255)
                                                  aTCGTTTTCTAACACGATGTGAATATTATCTGTGTCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGa                                                                                                     <  8:207070/148‑1 (MQ=255)
                                                      ttttCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCAt                                                                                                <  8:151501/149‑1 (MQ=255)
                                                      ttttCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTAAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCAt                                                                                                >  3:269609/1‑149 (MQ=255)
                                                        ttCTAACACGATGTTAATATTATCTGTTGCTAGATATTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCAtt                                                                                               <  3:295713/148‑1 (MQ=255)
                                                        ttCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCAtt                                                                                               <  1:347889/148‑1 (MQ=255)
                                                         tCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTg                                                                                              <  1:129780/148‑1 (MQ=255)
                                                           tAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTATAAATCTTTTCGCACTTGATAGAATATTTCTTTAAAAATGGCAACCTGAGCCATTgtt                                                                                            >  3:345408/1‑146 (MQ=255)
                                                             acacGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCAGAGTCTAAATCTTTTCGCACCTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTa                                                                                           <  5:187830/147‑1 (MQ=255)
                                                             acacGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGt                                                                                            <  7:269175/146‑1 (MQ=255)
                                                             acacGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGt                                                                                            <  1:316801/146‑1 (MQ=255)
                                                             acacGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTaa                                                                                          >  4:312179/1‑148 (MQ=255)
                                                               acgaTGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAAc                                                                                       <  2:173245/146‑1 (MQ=255)
                                                               acgaTGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTATTCGCACTTGATCGAATATTTc                                                                                                                        <  7:97460/113‑1 (MQ=255)
                                                                   tgtgAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTc                                                                                   <  4:337383/149‑1 (MQ=255)
                                                                     tgAATATTATCTGTGGATAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCa                                                                                 <  3:43035/149‑1 (MQ=255)
                                                                       aaTATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCAt                                                                                <  2:153257/148‑1 (MQ=255)
                                                                         tattatCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGTTAAAACCTTCCAtgt                                                                              >  5:283029/1‑148 (MQ=255)
                                                                         tattatCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCAtgtg                                                                             <  7:103288/149‑1 (MQ=255)
                                                                         tattatCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCAtgtg                                                                             <  7:322952/149‑1 (MQ=255)
                                                                          attatCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGa                                                                            >  1:7586/1‑149 (MQ=255)
                                                                           ttatCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGAt                                                                           <  4:286877/149‑1 (MQ=255)
                                                                              tCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATAc                                                                         >  3:187375/1‑148 (MQ=255)
                                                                                tgtgGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAg                                                                      >  1:101147/1‑149 (MQ=255)
                                                                                 gtgGCTAGATAGTAATTATAATGTGAGACGTTGTGACGTTTTAGTTAAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACATGAGCAATTGGTAAGACCTTCCATGTGATTCGTgg                                                                     >  6:272217/1‑149 (MQ=255)
                                                                                  tgGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGa                                                                       <  2:102736/146‑1 (MQ=255)
                                                                                    ggtAGATAGTAAATATAATGTGAGTCGTTGTGACGTTTTAGTTCAGAATAAAACAATTCAGATTCTAAATGTTTTCGCACTTGATCTAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGgc                                                                   <  6:111505/146‑1 (MQ=255)
                                                                                      tAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGcgcg                                                                >  1:48131/1‑149 (MQ=255)
                                                                                          tAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAg                                                             <  1:184826/148‑1 (MQ=255)
                                                                                            gTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGAAACCTGATCCATTGGTAAAACCATCCATGTGATACGAGGCCGCGTAGtat                                                          >  7:239603/1‑147 (MQ=255)
                                                                                             tAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTg                                                         >  6:367710/1‑149 (MQ=255)
                                                                                               aaTATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGc                                                        <  8:346510/148‑1 (MQ=255)
                                                                                               aaTATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGc                                                        <  2:39952/148‑1 (MQ=255)
                                                                                               aaTATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGc                                                        <  2:200417/148‑1 (MQ=255)
                                                                                                atatAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTc                                                                                   >  8:32264/1‑120 (MQ=255)
                                                                                                 tataATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAACCCTTCCATGTGATACTAGGGCGCGAAGTTTGCAt                                                      >  8:370000/1‑148 (MQ=255)
                                                                                                        tgAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCAtgt                                                                              >  5:50297/1‑117 (MQ=255)
                                                                                                        tgAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGtt                                                           <  7:75545/136‑1 (MQ=255)
                                                                                                        tgAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGttt                                              >  4:368383/1‑149 (MQ=255)
                                                                                                        tgAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGtt                                               <  8:39686/148‑1 (MQ=255)
                                                                                                         gagaCGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGtttt                                             >  5:313606/1‑149 (MQ=255)
                                                                                                                 gtgACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCgtt                                      <  5:343734/148‑3 (MQ=255)
                                                                                                                  tgACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTa                                           <  8:71086/142‑1 (MQ=255)
                                                                                                                    aCGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGATGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCaa                                  >  1:218887/1‑149 (MQ=255)
                                                                                                                      gTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATc                                >  1:381944/1‑149 (MQ=255)
                                                                                                                         ttAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTgg                             >  6:72110/1‑149 (MQ=255)
                                                                                                                          tAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTgg                             <  1:158063/148‑1 (MQ=255)
                                                                                                                           aGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGt                            <  1:268530/148‑1 (MQ=255)
                                                                                                                             ttCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTg                         <  3:238059/149‑1 (MQ=255)
                                                                                                                                aGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAAc                                                                                                         >  2:352333/1‑66 (MQ=255)
                                                                                                                                aGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGAcc                      >  7:42090/1‑149 (MQ=255)
                                                                                                                                aGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGAc                       <  4:166321/148‑1 (MQ=255)
                                                                                                                                 gAATAAAACAATTCACAGTCAAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACct                     >  1:294019/1‑149 (MQ=255)
                                                                                                                                  aaTAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTc                                                                                   >  8:266484/1‑86 (MQ=255)
                                                                                                                                  aaTAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATc                                >  8:292357/1‑137 (MQ=255)
                                                                                                                                   aTAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACctcc                   >  4:221174/1‑148 (MQ=255)
                                                                                                                                   aTAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACctcc                   <  2:22619/149‑2 (MQ=255)
                                                                                                                                      aaaCAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGAc                       >  3:212891/1‑142 (MQ=255)
                                                                                                                                       aaCAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCg                                        <  2:12941/124‑1 (MQ=255)
                                                                                                                                       aaCAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTtgt               >  7:186529/1‑149 (MQ=255)
                                                                                                                                       aaCAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTtgt               <  4:58481/149‑1 (MQ=255)
                                                                                                                                       aaCAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTtgt               >  7:293107/1‑149 (MQ=255)
                                                                                                                                        aCAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGt                                                            >  6:203942/1‑103 (MQ=255)
                                                                                                                                              cacaGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTttgt        >  7:211732/1‑149 (MQ=255)
                                                                                                                                               acaGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTTAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTttgtt       <  3:132073/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGA  >  NC_000913/579392‑579688

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