Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,202,959 T→A 8.5% intergenic (‑26/‑65) cohE ← / → croE e14 prophage; repressor protein phage e14/e14 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,202,9590TA 8.5% 436.4 / 23.4 130intergenic (‑26/‑65)cohE/croEe14 prophage; repressor protein phage e14/e14 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (60/59);  new base A (2/9);  total (62/68)
Fisher's exact test for biased strand distribution p-value = 5.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 140 total aligned reads displayed.

TCCCATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACG  >  NC_000913/1202811‑1203106
                                                                                                                                                    |                                                                                                                                                   
tcCCATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGtt                                                                                                                                                     >  3:368050/1‑149 (MQ=255)
 cccATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGtt                                                                                                                                                     >  1:571032/1‑148 (MQ=255)
  ccATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTa                                                                                                                                                    <  2:888830/148‑1 (MQ=255)
   cATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAg                                                                                                                                                  >  3:406992/1‑149 (MQ=255)
   cATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAg                                                                                                                                                  <  3:213207/149‑1 (MQ=255)
   cATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCCGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTaa                                                                                                                                                   <  2:775046/148‑1 (MQ=255)
    aTAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAgg                                                                                                                                                 <  4:368050/149‑1 (MQ=255)
    aTAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAgg                                                                                                                                                 <  4:526278/149‑1 (MQ=255)
    aTAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAg                                                                                                                                                  <  2:672388/148‑1 (MQ=255)
       aaCCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCCGTTTTCTGGCGTTAGACCGCCGGATGTCATGGAATGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTAc                                                                                                                                              >  3:447849/1‑149 (MQ=255)
         ccTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTa                                                                                                                                               >  2:3771/1‑146 (MQ=255)
         ccTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTa                                                                                                                                               <  1:3771/146‑1 (MQ=255)
         ccTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACaa                                                                                                                                            >  4:176722/1‑149 (MQ=255)
           tgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCTTAACGGTACAAgg                                                                                                                                          <  3:493003/149‑2 (MQ=255)
           tgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCTTAACGGTACAAgg                                                                                                                                          <  4:186393/149‑2 (MQ=255)
           tgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAg                                                                                                                                           >  2:815643/1‑148 (MQ=255)
            gctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCTTAACGGTACAAgg                                                                                                                                          <  1:290946/148‑2 (MQ=255)
              tgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCTTAACGGTACAAgggtt                                                                                                                                       <  3:926124/149‑5 (MQ=255)
              gggTCTGGGTAGGCTCGCGATCAATCATAGTGACCAAGGAATAGTTGCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCTTAACGGTACAAgggt                                                                                                                                        <  2:735240/147‑4 (MQ=255)
               ggTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATc                                                                                                                                       >  1:500622/1‑148 (MQ=255)
                gTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCtt                                                                                                                                     <  3:945566/149‑1 (MQ=255)
                  cTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTg                                                                                                                                    >  1:35398/1‑148 (MQ=255)
                   tGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGaa                                                                                                                                  <  3:573720/149‑1 (MQ=255)
                   tGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGaa                                                                                                                                  <  4:326140/149‑1 (MQ=255)
                   tGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGa                                                                                                                                   <  2:736434/148‑1 (MQ=255)
                      gTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAgg                                                                                                                                >  2:603519/1‑148 (MQ=255)
                        tGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTc                                                                                                                             <  4:663086/149‑1 (MQ=255)
                            tCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCAttt                                                                                                                         >  3:38935/1‑149 (MQ=255)
                             cgcgATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCAttt                                                                                                                         >  2:168617/1‑148 (MQ=255)
                             cgcgATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCAttt                                                                                                                         >  1:588046/1‑148 (MQ=255)
                               cgATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCaa                                                                                                                      >  3:841200/1‑149 (MQ=255)
                                gATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAAt                                                                                                                     >  3:928043/1‑149 (MQ=255)
                                 atcaatcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAAt                                                                                                                     >  2:907649/1‑148 (MQ=255)
                                  tcaatcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCa                                                                                                                   >  4:887171/1‑149 (MQ=255)
                                   caatcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCa                                                                                                                   >  2:331628/1‑148 (MQ=255)
                                      tcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGt                                                                                                                >  1:270496/1‑148 (MQ=255)
                                       caTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTaa                                                                                                              <  4:185708/149‑1 (MQ=255)
                                       caTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTaa                                                                                                              <  4:677112/149‑1 (MQ=255)
                                       caTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTa                                                                                                               <  2:588046/148‑1 (MQ=255)
                                                aaaGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAGCCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTacac                                                                                                     >  4:888793/1‑149 (MQ=255)
                                                aaaGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTacac                                                                                                     >  4:454476/1‑149 (MQ=255)
                                                 aaGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTacac                                                                                                     >  1:314742/1‑148 (MQ=255)
                                                    gaagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGa                                                                                                 <  4:480997/149‑1 (MQ=255)
                                                     aagacTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAAGTAAAACACTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAg                                                                                                >  4:20124/1‑149 (MQ=255)
                                                     aagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGGTGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGa                                                                                                 <  2:353334/148‑1 (MQ=255)
                                                     aagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACa                                                                                                                                             >  2:261758/1‑104 (MQ=255)
                                                     aagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACa                                                                                                                                             <  1:261758/104‑1 (MQ=255)
                                                     aagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGa                                                                                                 >  2:271427/1‑148 (MQ=255)
                                                      agaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAg                                                                                                >  1:401238/1‑148 (MQ=255)
                                                       gaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGt                                                                                              >  4:544843/1‑149 (MQ=255)
                                                            cccGCCAACACTATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACata                                                                                         >  3:762888/1‑149 (MQ=255)
                                                              cGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATAtt                                                                                       <  3:544843/149‑1 (MQ=255)
                                                                ccAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGt                                                                                     <  4:776614/149‑1 (MQ=255)
                                                                  aaCACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTAt                                                                                   <  3:309563/149‑1 (MQ=255)
                                                                       cATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAg                                                                              >  4:682037/1‑149 (MQ=255)
                                                                        aTCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAg                                                                              >  1:236746/1‑148 (MQ=255)
                                                                        aTCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGATCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTAt                                                                                   >  3:695676/1‑143 (MQ=255)
                                                                            ctcAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAGCCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTAt                                                                         <  3:888793/149‑1 (MQ=255)
                                                                             tcAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTAt                                                                         <  1:175390/148‑1 (MQ=255)
                                                                              cAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTaa                                                                                                                                                   <  1:143802/73‑1 (MQ=255)
                                                                                       ggCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCTTAACGGTACAAgg                                                                                                                                          <  1:849607/73‑2 (MQ=38)
                                                                                       ggCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCTTAACGGTACAAgg                                                                                                                                          >  2:849607/1‑72 (MQ=38)
                                                                                         cGTTAGACCGCTGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTaa                                                            <  3:619172/149‑1 (MQ=255)
                                                                                           ttAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTggg                                                                     <  1:738010/138‑1 (MQ=255)
                                                                                             aGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCaa                                                        <  3:700882/149‑1 (MQ=255)
                                                                                              gACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCaaa                                                       <  3:176722/149‑1 (MQ=255)
                                                                                              gACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGGATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCaa                                                        <  2:509716/148‑1 (MQ=255)
                                                                                                   ccgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAca                                                  <  4:447849/149‑1 (MQ=255)
                                                                                                     ggATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAcag                                                 <  1:815578/148‑1 (MQ=255)
                                                                                                         gTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACATCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAggg                                            <  3:97804/149‑1 (MQ=255)
                                                                                                          tCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGc                                           <  3:488618/149‑1 (MQ=255)
                                                                                                           cATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGc                                           >  1:904425/1‑148 (MQ=255)
                                                                                                             tGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGa                                         >  2:159300/1‑148 (MQ=255)
                                                                                                              ggATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAg                                        <  2:485700/148‑1 (MQ=255)
                                                                                                                 tttttttCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTgg                                    <  3:831659/146‑1 (MQ=255)
                                                                                                                 ttGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAata                                                                                                            <  2:623284/77‑1 (MQ=255)
                                                                                                                            aCGAAATTAATACCCTTGTACCGTTAAGGTAGAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGTGGTACATATGGTATGAAAGCGTATTGGGACGCTTTAACCAAAAAACAGCAAGGCGAGTTGGCCGGAAAAGtt                         >  3:589054/1‑149 (MQ=255)
                                                                                                                            aCGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGt                          <  2:264643/148‑1 (MQ=255)
                                                                                                                            aCGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTAGCCGGAAAAGtt                         >  4:911252/1‑149 (MQ=255)
                                                                                                                              gAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGCTgg                       >  4:7391/1‑149 (MQ=255)
                                                                                                                                 aTTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTc                    >  3:106538/1‑149 (MQ=255)
                                                                                                                                   tAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCa                   >  1:808197/1‑148 (MQ=255)
                                                                                                                                         ccTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATGTTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTg            <  3:296573/149‑1 (MQ=255)
                                                                                                                                         ccTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTg            <  4:928043/149‑1 (MQ=255)
                                                                                                                                          cTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTg            >  2:125203/1‑148 (MQ=255)
                                                                                                                                          cTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTg            >  2:398942/1‑148 (MQ=255)
                                                                                                                                           ttGTACCTTAACGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTGCATATTGTATGAAAGCGTATTGGGTCTCTTTAACCAAAGAACAGCAGGGCGAGTTGGcc                                  >  3:890772/1‑125 (MQ=255)
                                                                                                                                           ttGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCa                                                                                                                   >  2:554095/1‑44 (MQ=255)
                                                                                                                                           ttGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCa                                                                                                                   <  1:554095/44‑1 (MQ=255)
                                                                                                                                           ttGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGc          <  3:436099/149‑1 (MQ=255)
                                                                                                                                           ttGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGCTGGCTCAACACCTGGc          <  3:7391/149‑1 (MQ=255)
                                                                                                                                            tGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGc          >  2:672538/1‑148 (MQ=255)
                                                                                                                                              tACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTa        <  1:672538/148‑1 (MQ=255)
                                                                                                                                                  gTTAAGGTACAAGTATCTTGAAGGTTCATTTTAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTAc   >  3:821117/1‑149 (MQ=255)
                                                                                                                                                   ttAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATAtt                                                                                       <  3:516737/64‑1 (MQ=255)
                                                                                                                                                   ttAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTAtt                                                                        <  3:741423/79‑1 (MQ=255)
                                                                                                                                                   ttAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTAtt                                                                        >  4:741423/1‑79 (MQ=255)
                                                                                                                                                   ttAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAacac               <  2:680903/136‑1 (MQ=255)
                                                                                                                                                   ttAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAacac               >  1:680903/1‑136 (MQ=255)
                                                                                                                                                    tAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACg  >  1:461918/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCCCATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACG  >  NC_000913/1202811‑1203106

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