Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,240,291 C→A 6.6% A15S (GCT→TCT)  PP_3717 ← LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,240,2910CA 6.6% 251.0 / 4.4 76A15S (GCT→TCT) PP_3717LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (39/32);  new base A (2/3);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGTGATGATCAACG  >  NC_002947/4240142‑4240439
                                                                                                                                                        |                                                                                                                                                    
gCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGc                                                                                                                                                      >  2:1780815/1‑150 (MQ=255)
  cGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCt                                                                                                                                                    >  2:561940/1‑150 (MQ=255)
      ggCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCTGCAGTCTAGCCTGCCa                                                                                                                                                <  1:1246021/150‑1 (MQ=255)
       gCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCTGGCAGGCTAGACTGccgg                                                                                                                                               <  2:276536/150‑3 (MQ=255)
       gCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGAAGCGCGCCCATTGCGAGGCTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAg                                                                                                                                               <  1:237581/150‑1 (MQ=255)
         aGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGc                                                                                                                                             <  1:1780815/150‑1 (MQ=255)
           cGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGAGGTAAGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAgatcg                                                                                                                                           >  1:990964/1‑146 (MQ=255)
               cacaTACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGccc                                                                                                                                       <  1:1178958/150‑1 (MQ=255)
                acaTACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCcg                                                                                                                                      >  2:888152/1‑150 (MQ=255)
                   tcagccAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCgcgc                                                                                                                                   >  1:1542936/5‑150 (MQ=255)
                      tCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGcaca                                                                                                                                <  1:1580471/150‑1 (MQ=255)
                         aGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGcc                                                                                                                             <  2:1542936/150‑1 (MQ=255)
                              cTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTg                                                                                                                        >  1:1329815/1‑150 (MQ=255)
                              cTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTg                                                                                                                        >  1:922548/1‑150 (MQ=255)
                               tCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTgc                                                                                                                       >  1:1516383/1‑150 (MQ=255)
                                   gggTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCa                                                                                                                   >  2:1396099/1‑150 (MQ=255)
                                   cccgcttaccctcacaatGGGCGCGCTGCAGCAGTGGCACGTGGC‑‑‑GTTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCa                                                                                                                   <  1:258868‑M1/132‑1 (MQ=255)
                                          aCCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCAt                                                                                                            >  2:1054698/1‑150 (MQ=255)
                                                 tcacaatgggcgcgctgcagcagtggcacgtggcgttgcagcgcgcccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGc                                                                                                        <  2:258868‑M2/102‑1 (MQ=255)
                                                         GCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCTGGCAGGCTAGACTGCCGggcatga                                                                                                                                         <  1:1941586‑M1/106‑4 (MQ=255)
                                                         GCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCTGGCAGGCTAGACTGCCGggcatga                                                                                                                                         >  2:1941586‑M1/1‑103 (MQ=255)
                                                               CAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAgc                                                                                       >  2:1516383‑M1/1‑150 (MQ=255)
                                                                  cagcagtggcacgtggcgttgcagcgcgcccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAgc                                                                                       >  1:1625383‑M2/32‑150 (MQ=255)
                                                                  c‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGcc                                                                                    <  2:1329815/150‑1 (MQ=255)
                                                                       ccACGTGCCACTGTTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGt                                                                                  >  1:15095/1‑150 (MQ=255)
                                                                       ccACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCCCACAAACTTTTGCGCCAGGTTCATCGGCCTTGGCACAGGCCAAGCGCCGt                                                                                  <  1:1901249/150‑1 (MQ=255)
                                                                         aCGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTg                                                                                                    <  1:1561968/130‑1 (MQ=255)
                                                                         aCGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTg                                                                                                    >  2:1561968/1‑130 (MQ=255)
                                                                             gCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTc                                                                            >  2:761727/1‑150 (MQ=255)
                                                                                aCTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCAcg                                                                         >  1:672530/1‑150 (MQ=255)
                                                                                aCTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCAcg                                                                         >  1:1115718/1‑150 (MQ=255)
                                                                                    ctgcAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCAcgcgcg                                                                     <  1:761727/150‑1 (MQ=255)
                                                                                      gcagcgcgcccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGt                                                                   >  2:476643‑M2/12‑150 (MQ=255)
                                                                                           gcgcccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCaa                                                              <  1:1054698‑M2/144‑1 (MQ=255)
                                                                                           gcgcccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCaa                                                              <  1:1396099‑M2/144‑1 (MQ=255)
                                                                                               ccATTGTGAGGGTAAGCGGGGTATCGGGAGGGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAatgat                                                          >  1:1724341‑M2/3‑150 (MQ=255)
                                                                                               ccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAat                                                             >  1:1704604‑M2/3‑147 (MQ=255)
                                                                                               ccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAat                                                             <  2:1704604‑M2/145‑1 (MQ=255)
                                                                                                cATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGTGCCAGGTTCATCGGCCTTg                                                                                                    <  1:1034515‑M2/106‑1 (MQ=255)
                                                                                                cATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGTGCCAGGTTCATCGGCCTTg                                                                                                    >  2:1034515‑M2/2‑107 (MQ=255)
                                                                                                 aTTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGc                                                        >  1:1211834/1‑150 (MQ=255)
                                                                                                  ttGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCa                                                       >  1:1145190/1‑150 (MQ=255)
                                                                                                   tgtgAGGGTAAGCGGGGTATCGGGAGGGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCAc                                                      <  2:1724341/150‑1 (MQ=255)
                                                                                                   tgtgAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCAc                                                      >  1:135203/1‑150 (MQ=255)
                                                                                                    gtgAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACt                                                     >  2:757227/1‑150 (MQ=255)
                                                                                                     tgAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTg                                                    <  2:1625383/150‑1 (MQ=255)
                                                                                                          gTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCa                                               >  2:656102/1‑150 (MQ=255)
                                                                                                          gTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCa                                               >  1:349212/1‑150 (MQ=255)
                                                                                                           tAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAg                                              <  2:187644/150‑1 (MQ=255)
                                                                                                                 gggTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGtctc                                        <  1:888152/150‑1 (MQ=255)
                                                                                                                        gggAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCTACGAtc                                 <  1:476643/150‑1 (MQ=255)
                                                                                                                          gagaGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGAtctc                               >  1:327096/1‑150 (MQ=255)
                                                                                                                             agGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTg                                                    <  2:868836/126‑1 (MQ=255)
                                                                                                                             agGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTg                                                    >  1:868836/1‑126 (MQ=255)
                                                                                                                               gCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGAtc                                 >  1:868945/1‑143 (MQ=255)
                                                                                                                               gCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGAtc                                 <  2:868945/143‑1 (MQ=255)
                                                                                                                                   gggTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTg                      >  1:742301/1‑150 (MQ=255)
                                                                                                                                      tCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCaa                                                              <  2:796906/107‑1 (MQ=255)
                                                                                                                                      tCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCaa                                                              >  1:796906/1‑107 (MQ=255)
                                                                                                                                          gCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATgg               <  1:561940/150‑1 (MQ=255)
                                                                                                                                          gCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATgg               >  2:184384/1‑150 (MQ=255)
                                                                                                                                           cccGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATggg              >  2:1010639/1‑150 (MQ=255)
                                                                                                                                            ctggCAGGCTAGACTGCCGGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCAcgcgcg                                                                     <  2:1918634/92‑1 (MQ=255)
                                                                                                                                            ctggCAGGCTAGACTGCCGGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCAcgcgcg                                                                     >  1:1918634/3‑94 (MQ=255)
                                                                                                                                            ccGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGt             >  2:938945/1‑150 (MQ=255)
                                                                                                                                              ggCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAatga                                                           >  2:1897874/1‑102 (MQ=255)
                                                                                                                                              ggCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAatga                                                           <  1:1897874/102‑1 (MQ=255)
                                                                                                                                               gCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGtgat          <  2:135203/150‑1 (MQ=255)
                                                                                                                                                   tCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGTGATGATc      >  2:325545/1‑150 (MQ=255)
                                                                                                                                                   tCTAGCCTGCCAGGCACGCCCGCGCACAGCACTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACt                                                     >  1:126134/1‑103 (MQ=255)
                                                                                                                                                   tCTAGCCTGCCAGGCACGCCCGCGCACAGCACTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACt                                                     <  2:126134/103‑1 (MQ=255)
                                                                                                                                                    cTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAatgat                                                          >  2:691949/1‑97 (MQ=255)
                                                                                                                                                    cTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAatgat                                                          <  1:691949/97‑1 (MQ=255)
                                                                                                                                                      aGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGTGATGATCAAc   <  2:349212/150‑1 (MQ=255)
                                                                                                                                                       gCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGtc                                          >  1:966586/1‑110 (MQ=255)
                                                                                                                                                       gCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGtc                                          <  2:966586/110‑1 (MQ=255)
                                                                                                                                                       gCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGTGATGATCAACg  <  1:184384/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
GCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGTGATGATCAACG  >  NC_002947/4240142‑4240439

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