Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,431,928 T→C S960S (AGT→AGC 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,431,9280TC100.0% 252.8 / NA 68S960S (AGT→AGCstfRRac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/35);  total (33/35)

GTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTACGGCACTAAATCCACGAATAACACCGGGGCACATACAC  >  NC_000913/1431781‑1432073
                                                                                                                                                   |                                                                                                                                                 
gTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGt                                                                                                                                                <  2:366654/151‑1 (MQ=255)
 tGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTc                                                                                                                                               >  2:1081667/1‑151 (MQ=255)
  ggCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTc                                                                                                                                               <  2:909222/150‑1 (MQ=255)
   gCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGAGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCgg                                                                                                                                             <  1:1109019/151‑1 (MQ=255)
    ccATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCggg                                                                                                                                            <  1:859081/151‑1 (MQ=255)
     cATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGc                                                                                                                                           >  1:922523/1‑151 (MQ=255)
      aTCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGcc                                                                                                                                          >  1:1116154/1‑151 (MQ=255)
        cAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGt                                                                                                                                        <  2:616137/151‑1 (MQ=255)
          gATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTAt                                                                                                                                      <  1:883452/151‑1 (MQ=255)
             aCCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGt                                                                                                                                   <  2:770049/151‑1 (MQ=255)
               cGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGtct                                                                                                                                 >  2:304743/1‑151 (MQ=255)
                   ccGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAgg                                                                                                                             <  2:210905/151‑1 (MQ=255)
                       cTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAAc                                                                                                                          >  1:253862/1‑150 (MQ=255)
                         ggTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACa                                                                                                                         <  2:71378/149‑1 (MQ=255)
                         ggTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAgg                                                                                                                       <  2:424251/151‑1 (MQ=255)
                         ggTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAgg                                                                                                                       >  1:774833/1‑151 (MQ=255)
                            tATGCCCTGATGCAGGGGCAGGCTTTGGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGTTCCCTGATATGCGTGGCGGGACGATTAAGGGCAAACCTGCCAGCGGGCGGGCCGTATTGTCTCAGGAACAGGAc                                                                                                                     <  1:259612/150‑1 (MQ=255)
                             aTGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACgg                                                                                                                   >  2:998855/1‑151 (MQ=255)
                                    gATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTaa                                                                                                             >  2:171749/1‑150 (MQ=255)
                                       gCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATc                                                                                                          <  2:722974/150‑1 (MQ=255)
                                       gCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATc                                                                                                          <  1:397136/150‑1 (MQ=255)
                                       gCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCg                                                                                                         >  1:254921/1‑151 (MQ=255)
                                        cAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGc                                                                                                        >  1:645466/1‑151 (MQ=255)
                                         aGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGc                                                                                                        >  2:197504/1‑150 (MQ=255)
                                            ggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATAc                                                                                                    <  1:291406/151‑1 (MQ=255)
                                             gCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATAcc                                                                                                   >  1:1123296/1‑151 (MQ=255)
                                               aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATAccc                                                                                                  >  1:430660/1‑150 (MQ=255)
                                               aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATAccc                                                                                                  >  2:819726/1‑150 (MQ=255)
                                               aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCca                                                                                                 <  1:491803/151‑1 (MQ=255)
                                                 gCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCcaca                                                                                               <  1:825001/151‑1 (MQ=255)
                                                  cTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCcac                                                                                                >  2:753720/1‑149 (MQ=255)
                                                   ttttGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagc                                                                                             >  1:129237/1‑151 (MQ=255)
                                                       gACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcgcc                                                                                          >  2:117037/1‑150 (MQ=255)
                                                          aaaTCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcgccagcg                                                                                      >  1:803534/1‑151 (MQ=255)
                                                                gCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCa                                                                                    <  2:662315/147‑1 (MQ=255)
                                                                gCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATcc                                                                                 >  2:387267/1‑150 (MQ=255)
                                                                gCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCa                                                                                <  1:384818/151‑1 (MQ=255)
                                                                     cccGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTAc                                                                            <  1:349978/150‑1 (MQ=255)
                                                                      ccGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAg                                                                          <  1:462928/151‑1 (MQ=255)
                                                                         aaaCTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGAtt                                                                       <  2:774844/151‑1 (MQ=255)
                                                                          aaCTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGAttt                                                                      <  1:255229/151‑1 (MQ=255)
                                                                            cTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTgg                                                                    >  2:272142/1‑151 (MQ=255)
                                                                              tGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTgggg                                                                  >  2:266505/1‑151 (MQ=255)
                                                                               gCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGa                                                                 <  2:1158187/151‑1 (MQ=255)
                                                                                     gCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGaaaa                                                           <  2:8041/151‑1 (MQ=255)
                                                                                        tATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAAcc                                                         <  2:75119/150‑1 (MQ=255)
                                                                                            cGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACAtc                                                    >  1:330387/1‑151 (MQ=255)
                                                                                               cAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACAtcgt                                                  <  2:608147/150‑1 (MQ=255)
                                                                                               cAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACAt                                                     <  2:623311/147‑1 (MQ=255)
                                                                                                 ggCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATcgtcgt                                               >  1:599325/1‑151 (MQ=255)
                                                                                                  gCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGtt                                              <  1:628581/151‑1 (MQ=255)
                                                                                                    gTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTg                                            >  1:399569/1‑151 (MQ=255)
                                                                                                    gTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTg                                            >  1:1092518/1‑151 (MQ=255)
                                                                                                      gATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGAt                                          <  2:403825/151‑1 (MQ=255)
                                                                                                       aTCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGAtt                                         <  1:773465/151‑1 (MQ=255)
                                                                                                         cccTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTa                                        <  1:541509/150‑1 (MQ=255)
                                                                                                         cccTGATATGCGTGGCTGGACGAATAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTAc                                       <  2:703606/151‑1 (MQ=255)
                                                                                                             gATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTACGGc                                    >  1:441211/1‑150 (MQ=255)
                                                                                                               tatGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTACGGCACt                                 >  2:723384/1‑151 (MQ=255)
                                                                                                                     tGGCTGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTACGGCACTAAATcc                           >  2:295748/1‑151 (MQ=255)
                                                                                                                         tGGACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTACGGCACTAAATCCACg                        >  2:326083/1‑150 (MQ=255)
                                                                                                                          ggACGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTACGGCACTAAATCCACGa                       >  2:214955/1‑150 (MQ=255)
                                                                                                                             cGATTAAGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTACGGCACTAAATCCACGAATaa                   <  1:809291/151‑1 (MQ=255)
                                                                                                                                  aaGGGCAAACCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTACGGCACTAAATCCACGAATAACAcc               <  2:54718/150‑1 (MQ=255)
                                                                                                                                         aaCCTGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTACGGCACTAAATCCACGAATAACACCGGGGcac        >  1:1177550/1‑150 (MQ=255)
                                                                                                                                             tGCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTACGGCACTAAATCCACGAATAACACCGGGGCACATac    >  1:1009945/1‑150 (MQ=255)
                                                                                                                                              gCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTACGGCACTAAATCCACGAATAACACCGGGGCACATacac  >  2:166473/1‑151 (MQ=255)
                                                                                                                                              gCCAGCGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTACGGCACTAAATCCACGAATAACACCGGGGCACATaca   <  1:295450/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAACCACATCGTCGTTTGATTACGGCACTAAATCCACGAATAACACCGGGGCACATACAC  >  NC_000913/1431781‑1432073

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