Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,429,160 C→G 8.7% L38V (CTG→GTG)  stfR → Rac prophage; putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,1600CG 8.7% 299.8 / 18.3 92L38V (CTG→GTG) stfRRac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base C (40/44);  new base G (2/6);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAA  >  NC_000913/1429032‑1429308
                                                                                                                                |                                                                                                                                                    
ttttttGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGAt                                                                                                                                  <  4:64991/149‑1 (MQ=255)
  ttttGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATg                                                                                                                                 >  7:324970/1‑148 (MQ=255)
      gtgGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGcc                                                                                                                            >  5:139984/1‑149 (MQ=255)
       tgGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGcc                                                                                                                            <  3:322371/148‑1 (MQ=255)
       tgGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCg                                                                                                                           <  2:176766/149‑1 (MQ=255)
         gggTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTAAAGAACTGAACAATCCAGCTGAAAGCAAAAAGTAACAGCACCACGGTGGTGGTGACCACGCTGCCCTCAGAAAATCCGGATGAAGCCggg                                                                                                                         >  5:207159/1‑149 (MQ=255)
          ggTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGc                                                                                                                        <  5:114439/149‑1 (MQ=255)
            tGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGt                                                                                                                      <  5:132791/149‑1 (MQ=255)
              aaTATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGtt                                                                                                                     <  3:378080/148‑1 (MQ=255)
               atatGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTAc                                                                                                                   <  3:242169/149‑1 (MQ=255)
                  tGGCAGTAAAGATTTCAGGTGTTCTGAAAGACGTCACAGGAAAACCGGTACTGAACTGGACAATCCAGCTGAAACGAAAACGTAACAGCACCACGGTGGTTGTGAACACGCTGTCCTCAGAAATTCCGGATGAAGCCGGGCGTTACAg                                                                                                                 <  3:287775/148‑1 (MQ=255)
                    gCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCa                                                                                                               <  1:357244/148‑1 (MQ=255)
                     cAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATg                                                                                                             <  4:50403/149‑1 (MQ=255)
                      aGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCAAAATCAAGCTGAAAGCAAAACGTAACAGCACCAAGGTGGTGGTGAAAAAGCTGGCCTCATAAAATCCTTATGAAGACGGTCGTTACAGCATg                                                                                                             >  1:403860/1‑148 (MQ=255)
                       gTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGa                                                                                                           >  4:211108/1‑149 (MQ=255)
                       gTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGa                                                                                                           >  4:131771/1‑149 (MQ=255)
                        tAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGAc                                                                                                          >  2:77119/1‑149 (MQ=255)
                          aaGTTTTCAGGTGTACTGAAAGACGTCACAGGAAAACCGGTACATAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGATCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGt                                                                                                        <  6:207159/149‑1 (MQ=255)
                           aGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTTCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTTGCCTCAGAAAATCCGGATTAAGCCGGGCGTTACAGCATGGACGtt                                                                                                       <  7:242040/149‑1 (MQ=255)
                           aGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCgg                                                                                                                          >  2:309154/1‑130 (MQ=255)
                           aGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCgg                                                                                                                          <  1:309154/130‑1 (MQ=255)
                           aGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGtt                                                                                                       <  7:36967/149‑1 (MQ=255)
                           aGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGtt                                                                                                       <  3:211108/149‑1 (MQ=255)
                           aGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGtt                                                                                                       <  6:139984/149‑1 (MQ=255)
                           aGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGtt                                                                                                       <  1:196759/149‑1 (MQ=255)
                                tCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCTTTACAGCATGGACGTTGAGt                                                                                                   <  5:42562/148‑1 (MQ=255)
                                tCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACAGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTa                                                                                                  >  7:245071/1‑149 (MQ=255)
                                  aGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACg                                                                                                <  5:328771/149‑1 (MQ=255)
                                     tgtACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGt                                                                                              >  3:311530/1‑148 (MQ=255)
                                     agtccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTc                                                                                             <  5:42843/145‑1 (MQ=255)
                                      gtccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGGGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCa                                                                                            <  8:342172/146‑1 (MQ=255)
                                       tACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCa                                                                                            <  3:131771/148‑1 (MQ=255)
                                       tACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAACCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGTGGGTTAGAGCATGGACGTTGAGTACGGTCAg                                                                                           >  4:264369/1‑149 (MQ=255)
                                        aCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCCCGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGt                                                                                              >  5:301044/1‑145 (MQ=255)
                                        aCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGt                                                                                              <  6:301044/145‑1 (MQ=255)
                                         cTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTa                                                                                         <  3:108215/149‑1 (MQ=255)
                                          tGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTa                                                                                         >  6:204354/1‑148 (MQ=255)
                                          tGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTa                                                                                         >  4:76680/1‑148 (MQ=255)
                                          tGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTa                                                                                         >  4:87618/1‑148 (MQ=255)
                                          tGAAAGACGGCACAGGAAAACAGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTAc                                                                                        >  6:287203/1‑149 (MQ=255)
                                            aaaGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACa                                                                                       <  8:324970/148‑1 (MQ=255)
                                             aaGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGtg                                                                                                       >  7:378086/1‑130 (MQ=255)
                                             aaGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGtg                                                                                                       <  8:378086/131‑2 (MQ=255)
                                              aGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtg                                                                                    <  4:17006/149‑3 (MQ=255)
                                                     acaGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCt                                                                             <  3:360390/149‑1 (MQ=255)
                                                     acaGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTT‑TTCTg                                                                            >  2:187373/1‑149 (MQ=255)
                                                     acaGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGCTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGGTATTCt                                                                             >  4:296502/1‑149 (MQ=255)
                                                       aGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTAt                                                                                <  6:342721/144‑1 (MQ=255)
                                                       aGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTAt                                                                                >  5:342721/1‑144 (MQ=255)
                                                        ggAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACGACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGtt                                                                          <  3:87618/149‑1 (MQ=255)
                                                        ggAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGtt                                                                          >  1:105654/1‑149 (MQ=255)
                                                               cGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTtggtgg                                                                     >  6:39030/1‑147 (MQ=255)
                                                               cGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGaa                                                                   <  4:234807/149‑1 (MQ=255)
                                                                     aGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTc                                                             <  4:311530/149‑1 (MQ=255)
                                                                       aaCTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTccc                                                           >  4:181589/1‑149 (MQ=255)
                                                                        aCTGAACAATCCAGCTGAAAGCAAAACGTAACATCACCACGGTGGTGTTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACTTTGAGTACGGTCAGTACAGAGTTATTCTGTTGGTTGAAGTATTccc                                                           <  2:403860/148‑1 (MQ=255)
                                                                               aaTCCAGCTGAAAGCAAAACGTAACAGCAGCACGGTGGTTGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACTGTCATTAAAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTcac                                                   <  2:36532/149‑1 (MQ=255)
                                                                               aaTCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACa                                                                                                                  >  8:230629/1‑86 (MQ=255)
                                                                               aaTCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACa                                                                                                                  <  7:230629/86‑1 (MQ=255)
                                                                               aaTCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTca                                                    >  7:29058/1‑148 (MQ=255)
                                                                                   cAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGc                                               >  5:324201/1‑149 (MQ=255)
                                                                                   cAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGc                                               >  1:92829/1‑149 (MQ=255)
                                                                                    aGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGcc                                              <  3:76680/149‑1 (MQ=255)
                                                                                      cTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGGCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCgg                                            <  5:287203/149‑1 (MQ=255)
                                                                                         aaaGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGAc                                         >  2:19276/1‑149 (MQ=255)
                                                                                         aaaGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAAACCGGATGCAGCCGGGCGTTACAGCATGGACGTAGAGTACGGTCAGTACAGCGTTAATCTGTGGGTGGAAGGATTCCCGCCGTCACATGCCgggcc                                         >  2:330687/1‑147 (MQ=255)
                                                                                          aaGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGAcc                                        <  6:289743/149‑1 (MQ=255)
                                                                                          aaGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGAc                                         >  6:380930/1‑148 (MQ=255)
                                                                                           agccagaCGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCa                                       >  8:232231/7‑149 (MQ=255)
                                                                                             cAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCAt                                      >  8:262502/1‑148 (MQ=255)
                                                                                              aaaaCGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTAc                                                                                                 >  3:126322/1‑88 (MQ=255)
                                                                                              aaaaCGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTAc                                                                                                 <  4:126322/88‑1 (MQ=255)
                                                                                              aaaaCGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTAAGTACAGAGTTATTCTGTTGGTGGAAGGATTCCCGCCGTACCATGCCGGGACCATTa                                    >  6:249657/1‑149 (MQ=255)
                                                                                               aaaCGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTAc                                   >  8:406268/1‑149 (MQ=255)
                                                                                                aaCGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATACGGATGACGCCGGGCGTGACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTAGTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTAcc                                  >  5:5955/1‑149 (MQ=255)
                                                                                                   gTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTT‑TTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCgtgt                              <  1:187373/149‑1 (MQ=255)
                                                                                                    tAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCgtgt                              <  2:105654/149‑1 (MQ=255)
                                                                                                    tAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCgtgt                              <  8:29058/149‑1 (MQ=255)
                                                                                                           accacGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATg                           >  6:155590/1‑145 (MQ=255)
                                                                                                              acGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATtc                    <  7:406268/149‑1 (MQ=255)
                                                                                                               cGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTtg                                                                         <  7:272473/95‑1 (MQ=255)
                                                                                                               cGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTtg                                                                         >  8:272473/1‑95 (MQ=255)
                                                                                                               cGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATtct                   <  2:92829/149‑1 (MQ=255)
                                                                                                                ggtggtggtGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATtctc                  <  6:324201/149‑1 (MQ=255)
                                                                                                                   ggtggtgAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAAc               >  6:150792/1‑149 (MQ=255)
                                                                                                                    gtggtgAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAAcc              <  5:39030/149‑1 (MQ=255)
                                                                                                                        tgAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCgg           >  5:364277/1‑148 (MQ=255)
                                                                                                                        tgAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGt          <  3:181589/149‑1 (MQ=255)
                                                                                                                            cacGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGc      <  7:232231/149‑1 (MQ=255)
                                                                                                                             acGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGACAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCt     >  2:182976/1‑149 (MQ=255)
                                                                                                                               gCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGa   <  8:245071/149‑1 (MQ=255)
                                                                                                                               gCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGGTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTg    >  3:261433/1‑148 (MQ=255)
                                                                                                                                cTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATTGACGTTGAGTACGGGCAGTACAGCGTTATTCTGTGTGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGAGTCTCAACCCGGTACGCTGaa  <  5:249657/149‑1 (MQ=255)
                                                                                                                                |                                                                                                                                                    
TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAA  >  NC_000913/1429032‑1429308

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