Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,776,500 T→C 5.4% F39L (TTC→CTC)  yfjY → CP4‑57 prophage; putative DNA repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,776,5000TC 5.4% 224.7 / 3.3 77F39L (TTC→CTC) yfjYCP4‑57 prophage; putative DNA repair protein
Reads supporting (aligned to +/- strand):  ref base T (38/33);  new base C (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

GCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCTGTTTGCCGGTTCTATTAGCAGTACCGAGG  >  NC_000913/2776352‑2776647
                                                                                                                                                    |                                                                                                                                                   
gcgcCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCAtt                                                                                                                                                    <  1:1313276/150‑1 (MQ=255)
 cgcCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACATCGCACGTTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTc                                                                                                                                                   <  1:1329144/150‑1 (MQ=255)
 cgcCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTc                                                                                                                                                   <  2:1070944/150‑1 (MQ=255)
   ccATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCAc                                                                                                                                                 <  1:299208/150‑1 (MQ=255)
   ccATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCAc                                                                                                                                                 <  1:986054/150‑1 (MQ=255)
     aTTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCt                                                                                                                                               <  2:265722/150‑1 (MQ=255)
      ttATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTc                                                                                                                                              >  1:454865/1‑150 (MQ=255)
         tGCGCATTATTGAATAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTCACCGCACAACGCACGCTAAAACGAGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCAc                                                                                                                                           >  1:1330936/1‑150 (MQ=255)
          gcgcATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACt                                                                                                                                          <  1:257348/150‑1 (MQ=255)
                 attGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACCCCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTTACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCCCTCAGgctg                                                                                                                                   >  2:1089285/1‑150 (MQ=255)
                    gtttAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCACTCACCTCCACTCAGGCTGCTc                                                                                                                                >  1:876465/3‑150 (MQ=255)
                    gtttAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCACTCACCTCCACTCAGGCTGCTc                                                                                                                                >  1:880981/3‑150 (MQ=255)
                    gATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGCCTGCTc                                                                                                                                >  1:1172655/1‑150 (MQ=255)
                    gATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCGCCACTCAGGCTGCTc                                                                                                                                >  2:1137276/1‑150 (MQ=255)
                       tAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTg                                                                                                                             <  1:31747/150‑1 (MQ=255)
                            aTACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCCGGCTGCTCGTGACTgg                                                                                                                        >  2:332054/1‑150 (MQ=255)
                             tACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGgc                                                                                                                       <  1:1231997/150‑1 (MQ=255)
                                   tgatgGAACAGTCACTCATCCCACAGCCACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCCTTCACCTCCACTCAGGCTGCTCGTGACTGgctgacgc                                                                                                                 >  1:875974/1‑150 (MQ=255)
                                     atgGAACAGTCACTCATCCCACAGACCCCGGTACTTCCACTGCCCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGCCACCTGCGCGAAACAGGCGTGGCCTACGCCACAACTCAGGCTGCTCGTGACTGgctggtgggg                                                                                                               >  2:65464/1‑144 (MQ=255)
                                        gAACAGTCACTCATCCCACAGCCACCGGTACTTCCACTGACCGCACACCGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGaaa                                                                                                            >  2:1275026/1‑150 (MQ=255)
                                         aaCAGTCACTCATCCCACAGACACCGGTACTTCCACTGCCCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGaaaa                                                                                                           >  2:739758/1‑150 (MQ=255)
                                          aCAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAAt                                                                                                          <  1:169035/150‑1 (MQ=255)
                                           aaGTCACTCATCCCACAGACACCGGTACTTCCACTGACCTCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATg                                                                                                         <  2:204250/149‑1 (MQ=255)
                                            aGTCACACATCCCACAGCCACCGGTACTTCCACTGACCGCACACCGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCCCTCACGCTGCTCGTGACTGGCTGCAGCTGAAAATgg                                                                                                        >  2:1274507/1‑150 (MQ=255)
                                                aCTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCggg                                                                                                    <  2:17883/150‑1 (MQ=255)
                                                   cATCCCACAGACACCGATACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAGACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGCCTGGCTGAAGCTGAAAATGGCGGGGcg                                                                                                 >  2:1011467/1‑149 (MQ=255)
                                                       ccACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAg                                                                                             >  1:1234071/1‑150 (MQ=255)
                                                       ccACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAg                                                                                             >  1:276578/1‑150 (MQ=255)
                                                          caGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGt                                                                                          >  2:1222209/1‑150 (MQ=255)
                                                          caGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGt                                                                                          <  2:1159578/150‑1 (MQ=255)
                                                           aGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTg                                                                                         >  1:124/1‑150 (MQ=255)
                                                           aGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCATAACTCTGCTTGCCCGACACCTGCGCGAAACAGGCGTGGCATTCACCTACACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTg                                                                                         >  2:206732/1‑150 (MQ=255)
                                                                ccGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTgaagaa                                                                                    >  1:921295/1‑150 (MQ=255)
                                                                   gtccTTCCACTGACCGCACACCGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCAATCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAAttt                                                                                 <  2:838353/147‑1 (MQ=255)
                                                                   gTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAAttt                                                                                 >  1:763357/1‑150 (MQ=255)
                                                                     aCTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTCACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATCGCGGGGCTGGAGCGGGAAGAATTTat                                                                               <  2:58474/150‑1 (MQ=255)
                                                                            ctttccGCAAAACGGACTGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCTAAACAGGCGTGGCATTGACCTCCAATCAGGCTGGTCGTGAATGGCTGAAGCGTTAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTg                                                                        <  2:1172655/146‑1 (MQ=255)
                                                                              gACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGCCACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTa                                                                      >  1:165335/1‑150 (MQ=255)
                                                                                ccGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGATGAATTTATGATGCTGTATc                                                                    >  2:1080776/1‑150 (MQ=255)
                                                                                  gCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTg                                                                  >  1:765455/1‑150 (MQ=255)
                                                                                  gCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTg                                                                  >  1:1259463/1‑150 (MQ=255)
                                                                                    acaACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGaa                                                                <  2:631402/150‑1 (MQ=255)
                                                                                            cGGTAAAACGCGCCTTAACGCTGCTTGACAGACACCTCCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCACGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGATGTATCTGAACCAGCAGa                                                        >  1:727200/1‑150 (MQ=255)
                                                                                             ggTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGaa                                                       >  1:976461/1‑150 (MQ=255)
                                                                                             ggTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGaa                                                       >  1:606965/1‑150 (MQ=255)
                                                                                              gTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAAc                                                      >  1:1024844/1‑150 (MQ=255)
                                                                                               tAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAAcc                                                     <  2:492769/150‑1 (MQ=255)
                                                                                               tAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAAcc                                                     <  2:478484/150‑1 (MQ=255)
                                                                                                 aaaCGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCGGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAg                                                   <  2:848711/150‑1 (MQ=255)
                                                                                                    cgcgCCTTATCGGTGCTTGACCGACACCTGCGCGAAACTGGCGTGGCAATCACCTCCTCTGTGGCTGCTCGTGACTGGCTGAAGCTGAATTTGGCGGGGCTGGAGCGTGAAGTATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTg                                                <  1:1244465/150‑1 (MQ=255)
                                                                                                    cgcgCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTg                                                <  1:848150/150‑1 (MQ=255)
                                                                                                      cgcCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGAt                                              >  2:625595/1‑150 (MQ=255)
                                                                                                       gcCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCACTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGAtt                                             <  2:880981/150‑1 (MQ=255)
                                                                                                       gcCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCACTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGAtt                                             <  2:876465/150‑1 (MQ=255)
                                                                                                        ccTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTg                                            <  2:13402/150‑1 (MQ=255)
                                                                                                        ccTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTg                                            <  2:763357/150‑1 (MQ=255)
                                                                                                        ccTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTg                                            <  1:771771/150‑1 (MQ=255)
                                                                                                        ccTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAAACAGCAGAACCAGGTGAttt                                            >  2:122678/1‑149 (MQ=255)
                                                                                                         cTTAACGCTGCTTGACCGAAACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTGAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGc                                           <  1:168172/150‑1 (MQ=255)
                                                                                                                  gctTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAAATTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAAc                                  >  2:316364/1‑150 (MQ=255)
                                                                                                                   ctTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAAcc                                 >  2:790351/1‑150 (MQ=255)
                                                                                                                     tGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCt                               >  2:19526/1‑150 (MQ=255)
                                                                                                                       aCCGCCCCCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCGGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGACCCAGCAGCACCAGTTGATTGCCCACGACACCCTgg                             >  1:1064218/1‑149 (MQ=255)
                                                                                                                       aCCGACACCTGCGCGAACCAGGCGTGGCATTCACCACCACTCAGGCTGCTCGTGACTGGCTGAAGGTGAAAATGGCGGTGCTGGAGCGAGAAGAATTTATGATGCTGTATATGATCAAGCAGAACCAGTTGATGGCCCACGAAACCATGt                             >  2:685724/1‑150 (MQ=255)
                                                                                                                          gACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCTGTTTg                          <  2:456760/150‑1 (MQ=255)
                                                                                                                          gACACCTGCGCGAAACAGGCGGGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCTGTTTg                          <  1:605901/150‑1 (MQ=255)
                                                                                                                              ccTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCTGTTTGCCgg                      <  2:454865/150‑1 (MQ=255)
                                                                                                                                 gcgcGAAACAGGCGGGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCTGTTTGCCGGTTc                   <  1:885727/150‑1 (MQ=255)
                                                                                                                                        aCAGGCGTTGCATTCACCCCCCCTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCTGTTTGCCGGTTCTATTAGc            <  1:34511/150‑1 (MQ=255)
                                                                                                                                        aCAGGCGTGGCATTCACCTGCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCTGTTTGCCGGTTCTATTAGc            <  2:1139120/150‑1 (MQ=255)
                                                                                                                                          cggCGTGGCATTCACCTCCACTGCGGCTGCTCGTCACTGCCTGAAGCAGAAAATGGCGGGGCTGGAGCGTGAAGAATTAATGATGCTGTATCTGAACCAGCAGTACCAGTTGATTGCCCTCGAATCCCTGTTTGCCGGTTCTATTAGCAg          <  2:231647/149‑1 (MQ=255)
                                                                                                                                             cGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCTGTTTGCCGGTTCTATTAGCAGTAc       >  2:501966/1‑150 (MQ=255)
                                                                                                                                               tGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGAGTGAAGAATTTATGATGCTGTATCTGAACCATCAGAACCAGTTGATTGCCCACGAAACCCTGTTTGCCGGTTCTATTAGCAGTACCg     <  1:548811/150‑1 (MQ=255)
                                                                                                                                                  cATTCACCTCCACTCTGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCTGTTTGCCGGTTCTATTAGCAGTACCGAgg  <  2:276578/150‑1 (MQ=255)
                                                                                                                                                  cATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCTGTTTGCCGGTTCTATTAGCAGTACCGAgg  <  2:1234071/150‑1 (MQ=255)
                                                                                                                                                  cATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCTGTTTGCCGGTTCTATTAGCAGTACCGAgg  >  1:332359/1‑150 (MQ=255)
                                                                                                                                                  cATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCTGTTTGCCGGTTCTATTAGCAGTACCGAgg  >  1:1351972/1‑150 (MQ=255)
                                                                                                                                                  cATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCTGTTTGCCGGTTCTATTAGCAGTAACGAgg  >  2:811340/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGGAGCGTGAAGAATTTATGATGCTGTATCTGAACCAGCAGAACCAGTTGATTGCCCACGAAACCCTGTTTGCCGGTTCTATTAGCAGTACCGAGG  >  NC_000913/2776352‑2776647

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