Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,075,720 C→A 5.8% P314P (CCC→CCA)  yeeR → CP4‑44 prophage; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,075,7200CA 5.8% 217.8 / 2.7 69P314P (CCC→CCA) yeeRCP4‑44 prophage; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base C (25/40);  new base A (2/2);  total (27/42)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

ACATGGCAGGGTTTAACGAACAGTTGAAAGAGAACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTGCTCCTCTCTTTGTCGCGAAACA  >  NC_000913/2075572‑2075867
                                                                                                                                                    |                                                                                                                                                   
aCATGGCAGGGTTTAACGAACAGTTGAAAGAGAACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATccc                                                                                                                                                     <  2:27887/149‑1 (MQ=255)
 cATGGCAGGGTTTAACGTACTGTTGAAAGAGAACCTGTCATTCACACCTGATGAACTGATAACGCTCTTCGGGAACCGCCTGAATTTATGACCGGAAGTGAATGAGGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCt                                                                                                                                                    <  2:42393/149‑1 (MQ=255)
   tGGCAGGGTTTAACGAACAGTTGAAAGAGACCCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGGTGCCTGGACGGCGGTGACTGGTATCCCttt                                                                                                                                                  <  1:93820/149‑1 (MQ=255)
     gCAGGGTTTAACGAACAGTTGAAAGAGCACCTGTCATTCACACCTGATGAACTGCAACCGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTGTAGATCGTTGCCTGGACGGCGGTGACAGGTATCCCtttt                                                                                                                                                 >  6:264873/1‑148 (MQ=255)
         ggTTTAACGAACAGTTGAAAGAGAACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGaa                                                                                                                                            <  7:253888/149‑1 (MQ=255)
         ggTTTAACGAACAGTTGAAAGAGAACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGaa                                                                                                                                            <  3:42003/149‑1 (MQ=255)
           tttAACGAACAGTTGAAAGAGAACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAgg                                                                                                                                          <  1:306458/149‑1 (MQ=255)
               aCGAACAGTTGAAAGAGAACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCg                                                                                                                                       <  2:71628/148‑1 (MQ=255)
                    cAGTTGAAAGAGAACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGttt                                                                                                                                 <  8:317262/149‑1 (MQ=255)
                    cAGTTGAAAGAGAACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGttt                                                                                                                                 <  3:170644/149‑1 (MQ=255)
                     aGTTGAAAGAGCACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGtttt                                                                                                                                <  2:289803/149‑1 (MQ=255)
                     aGTTGAAAGAGAACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGACTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGtttt                                                                                                                                <  3:229227/149‑1 (MQ=255)
                       ttGAAAGAGAACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTg                                                                                                                               >  3:240088/1‑148 (MQ=255)
                       ttGAAAGAGAACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTg                                                                                                                               >  7:377299/1‑148 (MQ=255)
                       ttGAAAGAGAACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTg                                                                                                                               >  8:217181/1‑148 (MQ=255)
                         gAAAGAGAACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTa                                                                                                                             >  3:137075/1‑148 (MQ=255)
                             gagaACCTGTCATTCACACCTGATGAACTGTAAACGAGCTTCCGGAACCGCCTGAATTTATCACCGGAATCGAATGACGATTTTGTCGATCTTTGCCTGGACGGCGGTGACTGGTATCCCTTTTATGAAGGTCGCCGTTTTGTATCAtt                                                                                                                        <  1:171656/149‑1 (MQ=255)
                               gaACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAAGTTATCACCGGCAGCGAATGACGATGTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTcc                                                                                                                      <  5:264873/149‑1 (MQ=255)
                                  ccTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATc                                                                                                                   <  7:217181/149‑1 (MQ=255)
                                     gTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACg                                                                                                                >  1:321045/1‑149 (MQ=255)
                                       cATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGt                                                                                                               <  1:12033/148‑1 (MQ=255)
                                         ttCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGa                                                                                                            >  8:185941/1‑149 (MQ=255)
                                            acacCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGAt                                                                                                           >  2:149623/1‑147 (MQ=255)
                                              acCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGc                                                                                                       >  5:226716/1‑149 (MQ=255)
                                               ccTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCt                                                                                                      <  1:149623/149‑1 (MQ=255)
                                                   atgaACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCTTTTTGTATCATTCCATCACGTGGATGAGCTTCg                                                                                                   >  2:109910/1‑148 (MQ=255)
                                                   atgaACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTc                                                                                                    <  2:5696/147‑1 (MQ=255)
                                                   atgaACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTc                                                                                                    <  7:104073/147‑1 (MQ=255)
                                                    tgaACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGt                                                                                                  >  7:293649/1‑148 (MQ=255)
                                                        cTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTGCGTATCtg                                                                                             >  5:151505/1‑149 (MQ=255)
                                                         tGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCAGGTGGATGAGCTTCGTATCtgt                                                                                            >  6:237484/1‑149 (MQ=255)
                                                         tGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCtgt                                                                                            >  7:153687/1‑149 (MQ=255)
                                                         tGAAAACGATCTTCCGGAACCGACAGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGACGGTCGCCGTTTTGTATCCTTCCATCACGGGGATGAGCTTCGTATCtgt                                                                                            <  4:7890/149‑1 (MQ=255)
                                                           aaaaCGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCtgt                                                                                            <  3:313917/147‑1 (MQ=255)
                                                           aaaaCGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGc                                                                                          >  7:200837/1‑149 (MQ=255)
                                                           aaaaCGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGc                                                                                          <  1:109341/149‑1 (MQ=255)
                                                           aaaaCGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGACGAGCTTCGTATCtgt                                                                                            <  3:202892/147‑1 (MQ=255)
                                                              aCGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATc                                                                                                                   >  6:104314/1‑121 (MQ=255)
                                                              aCGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATc                                                                                                                   <  5:104314/121‑1 (MQ=255)
                                                              aCGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGcctc                                                                                       >  3:188778/1‑149 (MQ=255)
                                                              aCGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGATGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGAGGATGAGCTTCGTATCTGTGcct                                                                                        >  3:338232/1‑148 (MQ=255)
                                                                gCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGcctcct                                                                                     >  6:9045/1‑149 (MQ=255)
                                                                    ttCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGc                                                                                   <  1:151211/147‑1 (MQ=255)
                                                                      ccGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCggg                                                                                <  5:9045/148‑1 (MQ=255)
                                                                       cGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCt                                                                              <  3:255850/149‑1 (MQ=255)
                                                                         gAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGa                                                                            <  8:377299/149‑1 (MQ=255)
                                                                              gCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGaa                                                                       <  8:153687/149‑1 (MQ=255)
                                                                                cTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGt                                                                     <  4:137075/149‑1 (MQ=255)
                                                                                  gAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTAc                                                                   <  8:200837/149‑1 (MQ=255)
                                                                                   aaTTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTAc                                                                   <  7:344263/148‑1 (MQ=255)
                                                                                       tATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTAcatcat                                                              >  8:56020/1‑149 (MQ=255)
                                                                                       tATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTAca                                                                  <  7:18250/145‑1 (MQ=255)
                                                                                       tATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTAca                                                                  >  8:18250/1‑145 (MQ=255)
                                                                                                    gAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATc                                                 <  3:165330/149‑1 (MQ=255)
                                                                                                     aaTGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCa                                                <  6:226716/149‑1 (MQ=255)
                                                                                                     aaTGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCATTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCa                                                >  8:267421/1‑149 (MQ=255)
                                                                                                     aaTGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCATTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCa                                                >  8:272049/1‑149 (MQ=255)
                                                                                                       tGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATa                                              <  1:168375/149‑1 (MQ=255)
                                                                                                           gATTTTCTCGATCGTTGCCTGGACGGCGGTGACTTGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGcc                                          <  8:293649/149‑1 (MQ=255)
                                                                                                           gATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGcc                                          >  8:214070/1‑149 (MQ=255)
                                                                                                           gATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGcc                                          >  8:256028/1‑149 (MQ=255)
                                                                                                             ttttCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCg                                         >  7:191748/1‑148 (MQ=255)
                                                                                                                     aTCGTTGCCTGGTCGGCGGTGACTTGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTTGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTTACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGTATGCGGaa                                <  4:338232/149‑1 (MQ=255)
                                                                                                                        gTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATg                              <  2:321045/148‑1 (MQ=255)
                                                                                                                           gCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTa                           <  6:151505/148‑1 (MQ=255)
                                                                                                                             cTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGCCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATg                              <  2:97048/143‑1 (MQ=255)
                                                                                                                             cTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGACGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATg                              >  1:97048/1‑143 (MQ=255)
                                                                                                                             cTGGACGGCGGTGACTGGTATCCATTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTg                        <  7:272049/149‑1 (MQ=255)
                                                                                                                             cTGGACGGCGGTGACTGGTATCCATTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTg                        <  7:267421/149‑1 (MQ=255)
                                                                                                                                           cTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTGCTCCTCTCTTTGTc          <  4:188778/149‑1 (MQ=255)
                                                                                                                                                   ccTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTGCTCCTCTCTTTGTCGCgaaaca  <  4:240088/149‑1 (MQ=255)
                                                                                                                                                   ccTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTGCTCCTCTCTTTGTCGCgaaa    >  4:144161/1‑147 (MQ=255)
                                                                                                                                                    cTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTGCTCCTCTCTTTGTCGCgaaa    >  3:329881/1‑146 (MQ=255)
                                                                                                                                                    cTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTGCTCCTCTCTTTGTCGCgaaa    >  8:7806/1‑146 (MQ=255)
                                                                                                                                                    cTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGATTCGTATCTGTGCCTACTGCGGGGTGTAAGAAGTACATCATGCCCCGGAAAATAATAAGCAGGATCCGGAATGGTACTGCTCCTGTCTTTGTCGCgatt    >  7:327199/1‑144 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACATGGCAGGGTTTAACGAACAGTTGAAAGAGAACCTGTCATTCACACCTGATGAACTGAAAACGCTCTTCCGGAACCGCCTGAATTTATCACCGGAAGCGAATGACGATTTTCTCGATCGTTGCCTGGACGGCGGTGACTGGTATCCCTTTTCAGAAGGTCGCCGTTTTGTATCATTCCATCACGTGGATGAGCTTCGTATCTGTGCCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTGCTCCTCTCTTTGTCGCGAAACA  >  NC_000913/2075572‑2075867

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: ‑