Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,382,866 T→C 10.1% intergenic (‑33/+229) PP_3857 ← / ← PP_3858 hypothetical protein/tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,382,8660TC10.1% 185.8 / 9.5 70intergenic (‑33/+229)PP_3857/PP_3858hypothetical protein/tail fiber protein
Reads supporting (aligned to +/- strand):  ref base T (35/27);  new base C (5/2);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GTCCGTATCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCC  >  NC_002947/4382717‑4383013
                                                                                                                                                     |                                                                                                                                                    
gTCCGTATCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGt                                                                                                                                                      >  2:136385/1‑150 (MQ=255)
  ccGTATCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCCCCGTgg                                                                                                                                                    >  1:208119/1‑150 (MQ=255)
    gTATCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCg                                                                                                                                                  >  1:1002787/1‑150 (MQ=255)
         gACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGttt                                                                                                                                            >  2:559508/1‑150 (MQ=255)
                      aCATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGAtt                                                                                                                               <  1:559508/150‑1 (MQ=255)
                        aTAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCa                                                                                                                                      >  1:396114/1‑141 (MQ=255)
                        aTAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCa                                                                                                                                      <  2:396114/141‑1 (MQ=255)
                             cATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAg                                                                                                                        <  1:136385/150‑1 (MQ=255)
                               tgctgcAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTc                                                                                                                                       >  1:428376/1‑133 (MQ=255)
                               tgctgcAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTc                                                                                                                                       <  2:428376/133‑1 (MQ=255)
                                          cGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGAttt                                                                                                                              <  1:401486/131‑1 (MQ=255)
                                          cGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGAttt                                                                                                                              >  2:401486/1‑131 (MQ=255)
                                                  tctTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGctct                                                                                                   <  1:326880/150‑1 (MQ=255)
                                                         ttgcattgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGa                                                                                            >  1:1251751/1‑150 (MQ=255)
                                                           gcattgcCGCTTCGTACTACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGAtt                                                                                          >  1:685291/1‑150 (MQ=255)
                                                            cattgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAACGCTATCGCCGGGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTg                                                                                         >  1:79465/1‑150 (MQ=255)
                                                             attgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                          <  1:860345/100‑1 (MQ=255)
                                                             attgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                          >  2:860345/1‑100 (MQ=255)
                                                                 ccGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGa                                                                                            <  2:1061711/142‑1 (MQ=255)
                                                                 ccGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGa                                                                                            >  1:1061711/1‑142 (MQ=255)
                                                                          acaacaGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGcatca                                                                           >  2:1016713/1‑150 (MQ=255)
                                                                           caacaGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGcatcat                                                                          >  1:1099967/1‑150 (MQ=255)
                                                                                  ccTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCggcggc                                                                   <  2:1251751/150‑1 (MQ=255)
                                                                                   cTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCg                                                                                                                   >  1:474183/1‑101 (MQ=255)
                                                                                   cTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCg                                                                                                                   <  2:474183/101‑1 (MQ=255)
                                                                                     gAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCaaa                                                                                                                                    <  2:122828/82‑1 (MQ=255)
                                                                                     gAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCaaa                                                                                                                                    >  1:122828/1‑82 (MQ=255)
                                                                                     gAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAg                                                                <  2:308252/150‑1 (MQ=255)
                                                                                           gCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCa                                                                                                                         <  1:48759/87‑1 (MQ=255)
                                                                                           gCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCa                                                                                                                         >  2:48759/1‑87 (MQ=255)
                                                                                           gCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTAcgc                                                                                                             >  1:1031883/1‑99 (MQ=255)
                                                                                           gCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTAcgc                                                                                                             <  2:1031883/99‑1 (MQ=255)
                                                                                           gCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGttt                                                          >  2:201733/1‑150 (MQ=255)
                                                                                            cTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTa                                                         <  2:271533/150‑1 (MQ=255)
                                                                                            cTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTa                                                         <  1:201733/150‑1 (MQ=255)
                                                                                               ccAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCa                                                      >  1:51278/1‑150 (MQ=255)
                                                                                                 aaaGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGa                                                                                            >  2:1137370/1‑110 (MQ=255)
                                                                                                 aaaGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGa                                                                                            <  1:1137370/110‑1 (MQ=255)
                                                                                                 aaaGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCAtt                                                    >  2:610853/1‑150 (MQ=255)
                                                                                                            atCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTAc                                                                                                 <  2:273383/94‑1 (MQ=255)
                                                                                                            atCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTAc                                                                                                 >  1:273383/1‑94 (MQ=255)
                                                                                                                   aTTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTAGGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatccat                                  >  1:1202397/1‑150 (MQ=255)
                                                                                                                    ttCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatccatc                                 >  2:1209263/1‑150 (MQ=255)
                                                                                                                       aCTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAAt                              >  2:1001237/1‑150 (MQ=255)
                                                                                                                           ttCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCtt                          <  1:1016713/150‑1 (MQ=255)
                                                                                                                              tacAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCAt                       >  2:234016/3‑150 (MQ=255)
                                                                                                                                cAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGa                     >  1:630654/1‑150 (MQ=255)
                                                                                                                                 aGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGa                     >  1:1004452/1‑149 (MQ=255)
                                                                                                                                 aGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGa                     <  2:1004452/149‑1 (MQ=255)
                                                                                                                                   gCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTc                         <  1:1076602/143‑1 (MQ=255)
                                                                                                                                   gCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTc                         >  2:1076602/1‑143 (MQ=255)
                                                                                                                                   gCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGc                  <  2:1099967/150‑1 (MQ=255)
                                                                                                                                     aTCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCa                                                      >  2:334739/1‑112 (MQ=255)
                                                                                                                                     aTCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCa                                                      <  1:334739/112‑1 (MQ=255)
                                                                                                                                      tCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAGAACGTCATCGATTAATGCTTCATGACGTTCa               >  2:917100/1‑150 (MQ=255)
                                                                                                                                      tCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCa               <  2:208119/150‑1 (MQ=255)
                                                                                                                                      tCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCa               >  1:1274515/1‑150 (MQ=255)
                                                                                                                                       cAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACTCTCAt              <  1:1110700/150‑1 (MQ=255)
                                                                                                                                        cccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCAtt             >  1:617621/4‑150 (MQ=255)
                                                                                                                                        aaaaaaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCAtt             >  1:143291/1‑150 (MQ=255)
                                                                                                                                         ccaaaCCCGCCACGGCGGGTTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGGTTCATGACGACCAtt             >  2:123920/3‑150 (MQ=255)
                                                                                                                                         ccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTAc                                                        <  1:1092131/104‑1 (MQ=255)
                                                                                                                                         ccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTAc                                                        >  2:1092131/3‑106 (MQ=255)
                                                                                                                                         ccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGa                     <  2:255006/139‑1 (MQ=255)
                                                                                                                                         ccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGa                     >  1:255006/3‑141 (MQ=255)
                                                                                                                                         ccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTc            >  1:26831/3‑150 (MQ=255)
                                                                                                                                         aaaaaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTc            >  1:600281/1‑150 (MQ=255)
                                                                                                                                           aaccccGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGt          <  1:917100/147‑1 (MQ=255)
                                                                                                                                                cgccgTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCg     >  1:1085017/1‑150 (MQ=255)
                                                                                                                                                 gccgTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCgg    >  2:14992/1‑150 (MQ=255)
                                                                                                                                                   cgTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGCTCCCCACTGGCCAAGTGGTTTAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGcc  <  2:685291/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTCCGTATCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCC  >  NC_002947/4382717‑4383013

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