Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,383,558 G→C 6.9% T114S (ACT→AGT)  PP_3858 ← tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,383,5580GC 6.9% 227.2 / 6.7 72T114S (ACT→AGT) PP_3858tail fiber protein
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base C (2/3);  total (35/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 = 8.45e-01

TTTCAGTCTCGACCTTGGTATAGGTCTCTGCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAGTGGTGGCCTTGAGCGCGAAGTTGGCGACCATCTTCTGCGCCGC  >  NC_002947/4383409‑4383706
                                                                                                                                                      |                                                                                                                                                    
tttCAGTCTCGACCTTGGTATAGGTCTCTGCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAg                                                                                                                                                      >  2:521425/1‑150 (MQ=255)
tttCAGTCTCGACCTTGGTATAGGTCTCTGCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAg                                                                                                                                                      >  1:364135/1‑150 (MQ=255)
    aGTCTCGACCTTGGTATAGGTCTCTGCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGAt                                                                                                                                                  <  2:612716/150‑1 (MQ=255)
      tctcGACCTTGGTATAGGTCTCTGCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGc                                                                                                                                                >  1:1026191/1‑150 (MQ=255)
      tctcGACCTTGGTATAGGTCTCTGCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCACTGATGc                                                                                                                                                <  1:1149267/150‑1 (MQ=255)
       ctcGACCTTGGTATAGGTCTCTGCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGcc                                                                                                                                               <  2:1026191/150‑1 (MQ=255)
           aCCTTGGTATAGGTCTCTGCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTaa                                                                                                                                           >  1:1006055/1‑150 (MQ=255)
            ccTTGGTATAGGTCTCTGCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAac                                                                                                                                          <  2:836737/150‑1 (MQ=255)
              ttGGTATAGGTCTCTGCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAacca                                                                                                                                        <  1:1286327/150‑1 (MQ=255)
                    taGGTCTCTGCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGc                                                                                                                                  >  1:708219/1‑150 (MQ=255)
                          tctGCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAac                                                                                                                                          <  1:1035047/136‑1 (MQ=255)
                          tctGCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAac                                                                                                                                          >  2:1035047/1‑136 (MQ=255)
                             gCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGcc                                                                                                                         >  2:1095913/1‑150 (MQ=255)
                               tcaaCTGTACGTCTCGCCTTTGGTGGAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTACTGTTCAGCAGCGTGTTGATTTCTGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCa                                                                                                                       >  2:799626/4‑150 (MQ=255)
                                  aCTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGt                                                                                                                                                     >  2:999779/1‑117 (MQ=255)
                                  aCTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGt                                                                                                                                                     <  1:999779/117‑1 (MQ=255)
                                   cTGTACGTCTCGCCTTTGGTGTAAGCATACCAGATGACTCTCGGTGCCTGGGGTTTGGCGTTGTTCGCCTTGGTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTg                                                                                                                   <  1:318671/150‑1 (MQ=255)
                                    tGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCACTGATGcctac                                                                                                                                            >  1:968666/1‑121 (MQ=255)
                                    tGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCACTGATGcctac                                                                                                                                            <  2:968666/124‑4 (MQ=255)
                                      tACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAg                                                                                                                                                      >  2:716584/1‑112 (MQ=255)
                                      tACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAg                                                                                                                                                      <  1:716584/112‑1 (MQ=255)
                                       aCGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGAATCAGCCAAGCCTGAACCACCCAGCGTGATCGCCCAGTTGGCCt                                                                                                               >  2:1182942/1‑150 (MQ=255)
                                        cGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCtt                                                                                                              <  1:3824/150‑1 (MQ=255)
                                           ctcGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTg                                                                                                           >  2:845026/1‑150 (MQ=255)
                                              gCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGa                                                                                                        >  1:81863/1‑150 (MQ=255)
                                              gCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGa                                                                                                        <  1:1019216/150‑1 (MQ=255)
                                                 tttGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGAt                                                                                                                                                  <  2:692133/105‑1 (MQ=255)
                                                 tttGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGAt                                                                                                                                                  >  1:692133/1‑105 (MQ=255)
                                                  ttGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGacaac                                                                                                    >  2:993522/1‑150 (MQ=255)
                                                      tgtAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCACTGATGccta                                                                                                                                             >  1:179446/1‑103 (MQ=255)
                                                      tgtAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCACTGATGccta                                                                                                                                             <  2:179446/105‑3 (MQ=255)
                                                       gtAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTa                                                                                               >  1:721042/1‑150 (MQ=255)
                                                           gCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTcc                                                                                           <  2:1241649/150‑1 (MQ=255)
                                                            cATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCg                                                                                          >  2:549755/1‑150 (MQ=255)
                                                              tCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTc                                                                                        <  2:1420247/150‑1 (MQ=255)
                                                              tCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTc                                                                                        >  1:1371363/1‑150 (MQ=255)
                                                                aCCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTg                                                                                      >  2:444763/1‑150 (MQ=255)
                                                                    aTGCCATAGCCGGCCAAGGTGTTGTCGATCTTCTCCTTGCTGTTCAGCAGCGTGTTGCTTTCGGTTTTTCATTGTGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCAcc                                                                                  <  1:1095913/150‑1 (MQ=255)
                                                                      gCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCtt                                                                                <  2:561013/150‑1 (MQ=255)
                                                                       ccATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTg                                                                               <  1:1182942/150‑1 (MQ=255)
                                                                         aTAGCAGGCCAAGGTGTTGTCTTTTTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGTTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCATTTGGCCTTCTGCGACTACAGGTAGTCCGACTGCAccatagt                                                                             <  2:81863/150‑6 (MQ=255)
                                                                             ccGGCCAAGTTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACCCCTGGTAGTCCGTCTGCACCTTGGTGTAg                                                                         <  1:549755/150‑1 (MQ=255)
                                                                             ccGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAg                                                                         >  1:1308504/1‑150 (MQ=255)
                                                                             ccGGCCAAGGTATTGGCGTTGTTCGCCTTGCTTTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTGCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAg                                                                         >  2:67167/1‑150 (MQ=255)
                                                                                  ccAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGTGTGATCGCCCAGTTGGCCTTCTGCGACAACAgg                                                                                                 <  1:979060/122‑1 (MQ=255)
                                                                                  ccAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGTGTGATCGCCCAGTTGGCCTTCTGCGACAACAgg                                                                                                 >  2:979060/1‑122 (MQ=255)
                                                                                          ttGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAgg                                                                                                 >  2:937079/1‑114 (MQ=255)
                                                                                          ttGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAgg                                                                                                 <  1:937079/114‑1 (MQ=255)
                                                                                            ggCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAAt                                                              >  1:899862/1‑147 (MQ=255)
                                                                                            ggCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAAt                                                              <  2:899862/147‑1 (MQ=255)
                                                                                            ggCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAACCATGcc                                                           <  2:708219/150‑1 (MQ=255)
                                                                                            ggCGTTGTGCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGcc                                                           <  2:603285/150‑1 (MQ=255)
                                                                                             gCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGt                                                                                                                                             >  1:676408/1‑67 (MQ=255)
                                                                                             gCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGt                                                                                                                                             <  2:676408/67‑1 (MQ=255)
                                                                                             gCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCg                                                          <  1:521425/150‑1 (MQ=255)
                                                                                               gttgttCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTa                                                        >  2:1000410/1‑150 (MQ=255)
                                                                                                         ttGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGTGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCaaa                                              >  2:601648/1‑150 (MQ=255)
                                                                                                         ttGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCaac                                              >  2:254391/1‑149 (MQ=255)
                                                                                                            cTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAgtg                                           >  1:742575/1‑150 (MQ=255)
                                                                                                             tGTTCAGCAGCGGGTTGATTGCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAgtgg                                          <  2:220387/150‑1 (MQ=255)
                                                                                                               ttCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAgtggtg                                        <  1:1000410/150‑1 (MQ=255)
                                                                                                                tCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAGtggtgg                                       <  2:1006055/150‑1 (MQ=255)
                                                                                                                    cagcGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTg                                                                               <  2:867038/106‑1 (MQ=255)
                                                                                                                    cagcGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTg                                                                               >  1:867038/1‑106 (MQ=255)
                                                                                                                      gcGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAGTGGTGGCCTTGa                                 <  1:444763/150‑1 (MQ=255)
                                                                                                                              aTTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTa                                                                          <  2:630126/101‑1 (MQ=255)
                                                                                                                              aTTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTa                                                                          >  1:630126/1‑101 (MQ=255)
                                                                                                                                 tCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAGTGGTGGCCTTGAGCGCGAAGTTg                      <  1:993522/150‑1 (MQ=255)
                                                                                                                                   ggTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAGGAGGTATTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAGTGGTGGCCTTGAGCGCGAAGTTGGc                    <  1:799626/150‑1 (MQ=255)
                                                                                                                                     tCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAGTGGTGGCCTTGAGCGCGAAGTTGGc                    >  2:1310789/1‑148 (MQ=255)
                                                                                                                                     tCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAGTGGTGGCCTTGAGCGCGAAGTTGGc                    <  1:1310789/148‑1 (MQ=255)
                                                                                                                                     tCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAGTGGTGGCCTTGAGCGCGAAGTTGGCGa                  >  2:50088/1‑150 (MQ=255)
                                                                                                                                        tGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAGTGGTGGCCTTGAGCGCGAAGTTGGCGACCa               <  1:845026/150‑1 (MQ=255)
                                                                                                                                             tAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAGTGGTGGCCTTGAGCGCGAAGTTGGCGACCAtcttc          >  2:1273294/1‑150 (MQ=255)
                                                                                                                                                  aTCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAGTGGTGGCCTTGAGCGCGAAGTTGGCGACCATCTTCTgcgc     >  1:335331/1‑150 (MQ=255)
                                                                                                                                                     aGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGACGAACCCCGCCAAAGTGGTGGCCTTGAGCGCGAAGTTGGCGACCATCTTCTGcgccgc  >  1:234719/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TTTCAGTCTCGACCTTGGTATAGGTCTCTGCTCGACTGTACGTCTCGCCTTTGGTGTAAGCATCACCGATGCCATAGCCGGCCAAGGTGTTGGCGTTGTTCGCCTTGCTGTTCAGCAGCGTGTTGATTTCGGTCTTGGTGTAGGCATCAGTGATGCCGTAACCACCCAGCGTGATCGCCCAGTTGGCCTTCTGCGACAACAGGTAGTCCGTCTGCACCTTGGTGTAGGCATCAAGAATGCCGTACCCCGCCAAAGTGGTGGCCTTGAGCGCGAAGTTGGCGACCATCTTCTGCGCCGC  >  NC_002947/4383409‑4383706

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: