Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,453,523 A→G 7.0% intergenic (‑34/+79) PP_3069 ← / ← ppiC‑I outer membrane autotransporter/peptidyl‑prolyl cis‑trans isomerase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,453,5230AG 7.0% 201.3 / 3.5 71intergenic (‑34/+79)PP_3069/ppiC‑Iouter membrane autotransporter/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base A (32/34);  new base G (3/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

ACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCA  >  NC_002947/3453376‑3453672
                                                                                                                                                   |                                                                                                                                                     
aCCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCCGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACTAGGGAAACCTGTGCACCACGCGTGGTGcgca                                                                                                                                                     <  1:914983/150‑1 (MQ=255)
 cccGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCTAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAg                                                                                                                                                    >  1:1642919/1‑150 (MQ=255)
  ccgggccgggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGGCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAgg                                                                                                                                                   >  2:1069980/1‑150 (MQ=255)
   cgggccgggTGAGAAGGTACATGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCAGCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGTGCACCACGCGTGGTGCGCAggt                                                                                                                                                  <  1:78153/150‑2 (MQ=255)
   cgggccgggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAggg                                                                                                                                                  <  2:624490/150‑1 (MQ=255)
    gggccgggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGcgca                                                                                                                                                     <  1:818210/146‑1 (MQ=255)
        cgggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTAcc                                                                                                                                             >  2:811980/1‑150 (MQ=255)
          ggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCgg                                                                                                                                           >  1:515873/1‑150 (MQ=255)
            tGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCgggg                                                                                                                                         <  2:982500/150‑1 (MQ=255)
                   gTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGtata                                                                                                                                  >  2:692351/1‑150 (MQ=255)
                   gTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGtata                                                                                                                                  >  2:614751/1‑150 (MQ=255)
                   gTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGtata                                                                                                                                  <  1:614751/150‑1 (MQ=255)
                   gTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGtata                                                                                                                                  <  1:346819/150‑1 (MQ=255)
                     aCAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATaga                                                                                                                                >  1:1086112/1‑150 (MQ=255)
                      cAAGCCGCGTGTGCAGCGGGCGAGCAAATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATngag                                                                                                                               <  2:1186360/150‑1 (MQ=255)
                        aGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAgagga                                                                                                                             <  1:692351/150‑1 (MQ=255)
                        aGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAgagga                                                                                                                             <  2:1536811/150‑1 (MQ=255)
                         gCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGta                                                                                                                                    <  1:216990/142‑1 (MQ=255)
                         gCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGta                                                                                                                                    >  2:216990/1‑142 (MQ=255)
                         gCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAgaggag                                                                                                                            >  2:1232861/1‑150 (MQ=255)
                                    aGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGAc                                                                                                                 <  2:1086112/150‑1 (MQ=255)
                                       gggCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACcac                                                                                                              <  2:1210359/150‑1 (MQ=255)
                                            aGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTg                                                                                                         <  1:279788/150‑1 (MQ=255)
                                             ggCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTgg                                                                                                        >  2:517946/1‑150 (MQ=255)
                                             ggCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTgg                                                                                                        <  2:515873/150‑1 (MQ=255)
                                                    cccAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGACACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAAt                                                                                                 >  1:230258/1‑150 (MQ=255)
                                                           gcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGt                                                                                          <  2:515050/150‑1 (MQ=255)
                                                           gcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGt                                                                                          <  2:230258/150‑1 (MQ=255)
                                                           gcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCACCAGGATATCCTCCAGTCACCAGGGAAACCTGCGCCACACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGt                                                                                          <  1:1069980/150‑1 (MQ=255)
                                                           gcGGGAGCGAACCAAGGGTCAGCATCAGGCGTAGGCGCATCAGGAAAGCGTTTATGCAACAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGCTGTATAGAGGAGTGCGCCACGCGCAAAACTGGCAAAAATGTCAGGt                                                                                          >  2:741328/1‑150 (MQ=255)
                                                               gAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCa                                                                                      >  1:298335/1‑150 (MQ=255)
                                                                   gAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAgg                                                                                  >  1:30028/1‑150 (MQ=255)
                                                                     aCCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAgg                                                                                  >  2:754522/1‑148 (MQ=255)
                                                                     aCCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAgg                                                                                  <  1:754522/148‑1 (MQ=255)
                                                                        cAGGGTCAGCATCAGGCGTAGACACATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCTACTATGTCAGGTGGCAAAGGGCTAg                                                                             >  2:950808/1‑150 (MQ=255)
                                                                         aGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAgc                                                                            >  1:297906/1‑150 (MQ=255)
                                                                             tCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGt                                                                        <  1:1232861/150‑1 (MQ=255)
                                                                               aGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTa                                                                      <  2:298335/150‑1 (MQ=255)
                                                                               aGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTa                                                                      <  1:950808/150‑1 (MQ=255)
                                                                                    cAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCAAGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTcgc                                                                 <  2:1566694/150‑1 (MQ=255)
                                                                                     aGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTcgcg                                                                >  1:1147015/1‑150 (MQ=255)
                                                                                      ggCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCCGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCgg                                                               <  2:603600/150‑1 (MQ=255)
                                                                                         gTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTAACAAATTTGTAGCGCGTTAGTCGCGGTAg                                                            >  2:441384/1‑150 (MQ=255)
                                                                                               gCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCgg                                                               <  2:502073/141‑1 (MQ=255)
                                                                                               gCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCgg                                                               >  1:502073/1‑141 (MQ=255)
                                                                                                gaTCAGGAAATCCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCGAGCGCGTTAGTCGCGGTAGAACGttt                                                     <  2:1176196/149‑1 (MQ=255)
                                                                                                cATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGATCACAACTGGCAAAAATGTCAGGTGGCAAAGGTCTAGCGCGTTAGTCGCGGTAGAACGttt                                                     <  1:663962/150‑1 (MQ=255)
                                                                                                cATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGttt                                                     <  2:763004/150‑1 (MQ=255)
                                                                                                cATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGttt                                                     <  2:297906/150‑1 (MQ=255)
                                                                                                       aaaaCCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAg                                              <  2:669896/150‑1 (MQ=255)
                                                                                                          aCCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtg                                           >  2:294978/1‑150 (MQ=255)
                                                                                                          aCCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtg                                           <  2:437217/150‑1 (MQ=255)
                                                                                                            cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTCGCGCGTTAGTCCCTGTAGAACGTTTGCACCAggtggt                                         >  1:479781/1‑150 (MQ=255)
                                                                                                            cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtggt                                         >  2:165847/1‑150 (MQ=255)
                                                                                                            cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtggt                                         >  2:124598/1‑150 (MQ=255)
                                                                                                             tCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTa                                        <  2:170596/150‑1 (MQ=255)
                                                                                                               cATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTcgc                                                                 <  1:1690502/123‑1 (MQ=255)
                                                                                                               cATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTcgc                                                                 >  2:1690502/1‑123 (MQ=255)
                                                                                                                       aGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTg                              >  1:1312379/1‑150 (MQ=255)
                                                                                                                        gggAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGc                             <  1:811980/150‑1 (MQ=255)
                                                                                                                        gggAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGc                             <  1:1894853/150‑1 (MQ=255)
                                                                                                                          gAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGctc                           <  1:1412200/150‑1 (MQ=255)
                                                                                                                          gAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGctc                           <  2:1147015/150‑1 (MQ=255)
                                                                                                                             aCCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTg                        >  2:1803758/1‑150 (MQ=255)
                                                                                                                              ccTGTGCACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGttt                                                     >  1:923387/1‑120 (MQ=255)
                                                                                                                              ccTGTGCACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGttt                                                     <  2:923387/120‑1 (MQ=255)
                                                                                                                               cTGTGCACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGAt                      >  2:91778/1‑150 (MQ=255)
                                                                                                                               cTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGAt                      >  2:210761/1‑150 (MQ=255)
                                                                                                                                tgtgcACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATc                     >  1:1161271/4‑150 (MQ=255)
                                                                                                                                     accacGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGccc                >  2:162652/1‑150 (MQ=255)
                                                                                                                                     accacGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGccc                >  1:93515/1‑150 (MQ=255)
                                                                                                                                     accacGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGccc                >  1:927450/1‑150 (MQ=255)
                                                                                                                                          gcgTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTa           >  2:858441/1‑150 (MQ=255)
                                                                                                                                                   acaGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAATGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCa  >  2:1092989/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCA  >  NC_002947/3453376‑3453672

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