Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,435,965 G→A 6.0% Q278Q (CAG→CAA)  gpJ → baseplate/tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,435,9650GA 6.0% 235.1 / 2.4 78Q278Q (CAG→CAA) gpJbaseplate/tail fiber protein
Reads supporting (aligned to +/- strand):  ref base G (40/33);  new base A (2/3);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAAGGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGCTGACGTCGAAGGTGGTGACATGAGC  >  NC_002947/3435817‑3436114
                                                                                                                                                    |                                                                                                                                                     
ggCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAAGGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGt                                                                                                                                                      >  2:871115/1‑150 (MQ=255)
 gCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAAGGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGtt                                                                                                                                                     <  1:1919644/150‑1 (MQ=255)
  ctccGCGTGGCCGCCGGCGTACCGACTGAGTTCAAGGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTg                                                                                                                                                    >  2:1451176/3‑150 (MQ=255)
         tGGCCGCCGGCGTACCGACTGAGTTAACGGTGGAGGGGGTGGTGGCGCTCGGGGCTTGAGGATGCCGTGACGTGGTAAAGGCCACGGCGCGGGCTGGCCTCAACCCGGCCACTCCAGACGACAGCCTCCAGCAGGTGCTAGTGCCCCTgg                                                                                                                                             <  2:909890/150‑2 (MQ=255)
          ggCCGCCGGCGTACCGACTGAGTTCAAGGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTc                                                                                                                                            <  1:1451176/150‑1 (MQ=255)
            ccgccgGCGTACCGACTGAGTTCAAGGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCGACAGGGGCAACTGCCCCTgcagg                                                                                                                                          <  2:232428/150‑5 (MQ=255)
                         gACTGAGTTCAAGGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTgca                                                                                                                                            <  2:1725339/135‑3 (MQ=255)
                         gACTGAGTTCAAGGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTgca                                                                                                                                            >  1:1725339/1‑133 (MQ=255)
                         gACTGAGTTCAAGGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCcgc                                                                                                                             <  2:871175/150‑1 (MQ=255)
                             gAGTTCAAGGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTa                                                                                                                         <  1:1724694/150‑1 (MQ=255)
                              aGTTCAAGGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTac                                                                                                                        >  1:1664167/1‑150 (MQ=255)
                                ttCAAGGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTacac                                                                                                                      >  1:553650/1‑150 (MQ=255)
                                        ggAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACggg                                                                                                              >  2:1065722/1‑150 (MQ=255)
                                          aGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGa                                                                                                            >  2:681724/1‑150 (MQ=255)
                                              cGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCGGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGGGTTACACGTAACGGGGATTCg                                                                                                        >  2:453546/1‑150 (MQ=255)
                                              cGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCg                                                                                                        >  1:1195672/1‑150 (MQ=255)
                                               gAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGt                                                                                                       >  1:1210303/1‑150 (MQ=255)
                                                    gCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACg                                                                                                                <  2:1861687/136‑1 (MQ=255)
                                                    gCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACg                                                                                                                >  1:1861687/1‑136 (MQ=255)
                                                     cAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGccc                                                                                                                                                 >  2:1565277/1‑102 (MQ=255)
                                                     cAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGccc                                                                                                                                                 <  1:1565277/102‑1 (MQ=255)
                                                       gATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGAc                                                                                               <  2:553650/150‑1 (MQ=255)
                                                       gATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGAc                                                                                               <  2:1664167/150‑1 (MQ=255)
                                                       gATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGAc                                                                                               <  1:180607/150‑1 (MQ=255)
                                                        aTCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGAcc                                                                                              <  1:1065722/150‑1 (MQ=255)
                                                             ggAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGaa                                                                                         >  2:134947/1‑150 (MQ=255)
                                                              gAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAg                                                                                        >  2:1328158/1‑150 (MQ=255)
                                                              gAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAg                                                                                        <  1:134947/150‑1 (MQ=255)
                                                                 gggCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCAt                                                                                     >  1:1266395/1‑150 (MQ=255)
                                                                   gCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGc                                                                                       >  2:1253830/1‑146 (MQ=255)
                                                                   gCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGc                                                                                       <  1:1253830/146‑1 (MQ=255)
                                                                    ccAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACgg                                                                                                               >  1:631558/1‑121 (MQ=255)
                                                                    ccAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACgg                                                                                                               <  2:631558/121‑1 (MQ=255)
                                                                     cAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAg                                                                                 >  2:1060321/1‑150 (MQ=255)
                                                                      aGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGc                                                                                <  1:681724/150‑1 (MQ=255)
                                                                       gACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCAACCAGCg                                                                               >  1:581107/1‑150 (MQ=255)
                                                                              aCGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCgggg                                                                        <  2:1266395/150‑1 (MQ=255)
                                                                                      aaGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCATCGCAGGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGgtgtgt                                                                <  1:1087596/150‑1 (MQ=255)
                                                                                           cACGGCGCGGGCTGGCCTGAATGTGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATaa                                                           >  2:1052476/1‑150 (MQ=255)
                                                                                                 gcgGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTacac                                                                                                                      >  1:1739766/1‑85 (MQ=255)
                                                                                                 gcgGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTacac                                                                                                                      <  2:1739766/85‑1 (MQ=255)
                                                                                                 gcgGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCa                                                     <  1:1060321/150‑1 (MQ=255)
                                                                                                 gcgGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCGGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCa                                                     <  2:146262/150‑1 (MQ=255)
                                                                                                       tGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGc                                                                                                    <  2:44755/97‑1 (MQ=255)
                                                                                                       tGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGc                                                                                                    >  1:44755/1‑97 (MQ=255)
                                                                                                       tGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGAAGTTGCCGCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAGGCATCCAGCGGCCGGGGTGGCGTGTGATAAGCGGCATTAccc                                               >  2:815938/1‑150 (MQ=255)
                                                                                                         gCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAg                                                                                                                                                       <  2:456682/44‑1 (MQ=255)
                                                                                                         gCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAg                                                                                                                                                       >  1:456682/1‑44 (MQ=255)
                                                                                                           cTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGt                                           >  2:1540483/1‑150 (MQ=255)
                                                                                                           cTGAATGCGGACATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACAAGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGt                                           >  1:621304/1‑150 (MQ=255)
                                                                                                            tGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGtt                                          >  1:1639423/1‑150 (MQ=255)
                                                                                                              aaTGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGc                                        >  2:1856979/1‑150 (MQ=255)
                                                                                                               aTGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCa                                       >  1:1483170/1‑150 (MQ=255)
                                                                                                                   ggCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGAt                                                             <  2:976393/124‑1 (MQ=255)
                                                                                                                   ggCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGAt                                                             >  1:976393/1‑124 (MQ=255)
                                                                                                                      cATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCggg                                                                         >  1:1258197/1‑109 (MQ=255)
                                                                                                                      cATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCggg                                                                         <  2:1258197/109‑1 (MQ=255)
                                                                                                                      cATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGt                                >  1:1429463/1‑150 (MQ=255)
                                                                                                                       aTTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCgg                                                                          <  1:614111/107‑1 (MQ=255)
                                                                                                                       aTTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCgg                                                                          >  2:614111/1‑107 (MQ=255)
                                                                                                                        ttGCAGAAGCCCGTATCCTGCAGGGGCAGTTCCCCCTGTCCGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTCGAAGCATCCAGCTGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCATCTGCACCTCGGCTa                              <  2:671268/150‑1 (MQ=255)
                                                                                                                               aGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGCTGAc                       <  1:1052476/150‑1 (MQ=255)
                                                                                                                                gCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGCTGAc                       >  1:60816/1‑149 (MQ=255)
                                                                                                                                gCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGCTGAc                       <  2:60816/149‑1 (MQ=255)
                                                                                                                                gCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGCTGACg                      >  1:391760/1‑150 (MQ=255)
                                                                                                                                  ccGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGCTGACGTc                    >  2:857885/1‑150 (MQ=255)
                                                                                                                                     tCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGCTGACGTCGaa                 >  2:1345100/1‑150 (MQ=255)
                                                                                                                                        ccgacAGGGGCAACTGCCCCTGCCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGCTGACGTCGAAggt              >  2:735227/5‑150 (MQ=255)
                                                                                                                                        ccTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGCTGACGTCGAAggt              >  1:1272416/1‑150 (MQ=255)
                                                                                                                                         cgacAGGGGCAACTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGc                           >  2:1325221/4‑136 (MQ=255)
                                                                                                                                         cgacAGGGGCAACTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGc                           <  1:1325221/133‑1 (MQ=255)
                                                                                                                                           gCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGCTGACGTCGAAggtggt           <  2:1210303/150‑1 (MQ=255)
                                                                                                                                           gCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGCTGACGTCGAAggtggt           <  1:825226/150‑1 (MQ=255)
                                                                                                                                           gCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACACTGACGTCGAAggtggt           <  1:871115/150‑1 (MQ=255)
                                                                                                                                            cAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACgg                                                                                                               >  2:156968/1‑49 (MQ=255)
                                                                                                                                            cAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACgg                                                                                                               <  1:156968/49‑1 (MQ=255)
                                                                                                                                                ggCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGtt                                          >  1:127476/1‑114 (MQ=255)
                                                                                                                                                ggCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGtt                                          <  2:127476/114‑1 (MQ=255)
                                                                                                                                                 gCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCgg                                                       >  2:1455472/1‑100 (MQ=255)
                                                                                                                                                 gCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCgg                                                       <  1:1455472/100‑1 (MQ=255)
                                                                                                                                                    gTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGCTGACGTCGAAGGTGGTGACATGAGc  <  2:742396/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAAGGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTCGTAGCGTGACCTTGAAGCATCCAGCGGCCGGGGTGGTGTGTGATAAGCGGCATTACCCCAGTTGCACGTCGGTTACGCTGACGTCGAAGGTGGTGACATGAGC  >  NC_002947/3435817‑3436114

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