Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,088,847 G→C 5.7% G338A (GGC→GCC)  aphA → acetylpolyamine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,088,8470GC 5.7% 213.2 / 4.8 70G338A (GGC→GCC) aphAacetylpolyamine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base G (28/38);  new base C (2/2);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01

ACCCCATCTCGCAGTTCAAGCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCTGCATGCACAAGCCGAGCAGCGCACCCTGCGTGTATATAACTGGTTCGATTACATCAC  >  NC_002947/6088700‑6088996
                                                                                                                                                   |                                                                                                                                                        
acCCCATCTCGCAGTTCAAGCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCt                                                                                                                                                        <  1:900939/150‑1 (MQ=255)
 ccccATCTCGCAGTTCAAGCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCtt                                                                                                                                                       >  1:793371/1‑150 (MQ=255)
  cccATCTCGCAGTTCAAGCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGCCTTc                                                                                                                                                      >  1:796670/1‑150 (MQ=255)
        tcGCAGTTCAAGCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAgcgc                                                                                                                                              <  2:251617/150‑1 (MQ=255)
        tcGCAGTTCAAGCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAgcgc                                                                                                                                              <  2:1525555/150‑1 (MQ=255)
        tcGCAGTTAAAGCTGGACAGCCCGGATCACCTGGAAATGGGCAAGCGCATCGTCCAACTCGGCAAGTCGACCCTGTTCGTGATGGAAGGGGGCTACGCTGTGGAAGAAATCGGCACCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAgcgc                                                                                                                                              <  1:636222/150‑1 (MQ=255)
         cGCAGTTCAAGCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAgcgcg                                                                                                                                             <  1:600370/150‑1 (MQ=255)
         cGCAGTTCAAGCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAgcgcg                                                                                                                                             <  2:578399/150‑1 (MQ=255)
         cGCAGTTCAAGCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGCCTTC‑‑CAgcaca                                                                                                                                             <  2:796670/150‑4 (MQ=255)
         cGCAGTTCAAGCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAAGCGGCATCAACGCGGTCAATGCGCTGGAAGCCTTC‑‑CAgcaca                                                                                                                                             <  1:986492/150‑4 (MQ=255)
                cAAGCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTc                                                                                                                                      <  1:94019/150‑1 (MQ=255)
                 aaGCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCa                                                                                                                                     <  1:1161081/150‑1 (MQ=255)
                   gCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAg                                                                                                                                    >  2:1643476/1‑149 (MQ=255)
                   gCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAg                                                                                                                                    <  1:1643476/149‑1 (MQ=255)
                   gCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCGGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAgg                                                                                                                                   >  2:743493/1‑150 (MQ=255)
                    cTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGa                                                                                                                                  <  2:1887874/150‑1 (MQ=255)
                       gACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGcc                                                                                                                              <  1:1316129/150‑1 (MQ=255)
                        aCAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGccc                                                                                                                             >  2:289958/1‑150 (MQ=255)
                            cccGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATa                                                                                                                         >  1:517663/1‑150 (MQ=255)
                            cccGGATGACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATa                                                                                                                         >  1:828162/1‑150 (MQ=255)
                               ggATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGATTGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATg                                                                                                                      >  1:989490/1‑150 (MQ=255)
                                   tACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTcc                                                                                                                  <  2:517663/150‑1 (MQ=255)
                                      cTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGTAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGt                                                                                                                                       <  2:1426857/127‑1 (MQ=255)
                                      cTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGt                                                                                                                                       >  1:1426857/1‑127 (MQ=255)
                                       tGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACt                                                                                                              >  1:1544506/1‑150 (MQ=255)
                                       tGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACt                                                                                                              <  2:1544506/150‑1 (MQ=255)
                                        ggAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTc                                                                                                             <  2:793371/150‑1 (MQ=255)
                                            acgggCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCCTCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAAc                                                                                                         <  1:565965/148‑1 (MQ=255)
                                            aTGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAAc                                                                                                         <  2:1653241/150‑1 (MQ=255)
                                                gCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAgcgc                                                                                                                                              <  1:637622/110‑1 (MQ=255)
                                                gCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAgcgc                                                                                                                                              >  2:637622/1‑110 (MQ=255)
                                                gCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCACGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCt                                                                                                     >  2:781056/1‑150 (MQ=255)
                                                        aTCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGc                                                                                                   >  1:1176307/1‑144 (MQ=255)
                                                        aTCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGc                                                                                                   <  2:1176307/144‑1 (MQ=255)
                                                        aTCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGccc                                                                                             >  2:1470052/1‑150 (MQ=255)
                                                             ccAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTc                                                                                        <  1:1470052/150‑1 (MQ=255)
                                                              cAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAATAAATCGGCAACCACGCGGACCAAGCTCGTATGCCGTCTTCTGCTTGACAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTCCTCGCCCCGTTCa                                                                                       >  1:443012/1‑150 (MQ=255)
                                                               aaCTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCt                                                                                                                                                        <  1:1408343/87‑1 (MQ=255)
                                                               aaCTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCGGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCt                                                                                                                                                        >  2:1408343/1‑87 (MQ=255)
                                                                       aaGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCa                                                                              >  1:886894/1‑150 (MQ=255)
                                                                       aaGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCa                                                                              <  2:828162/150‑1 (MQ=255)
                                                                               ccTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCCAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTc                                                                      <  2:1122971/150‑1 (MQ=255)
                                                                               ccTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTc                                                                      <  2:923574/150‑1 (MQ=255)
                                                                                   ttCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAg                                                                  <  1:289958/150‑1 (MQ=255)
                                                                                      gTGATTGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGgcg                                                               >  2:919987/1‑150 (MQ=255)
                                                                                         aTGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGCCGccc                                                            >  1:30061/1‑150 (MQ=255)
                                                                                          tGGAAGGCGGCTACGCTGTTGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCATGAGCCCGATAATGGTCCGACTCAANCGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGcccn                                                           <  1:96346/150‑2 (MQ=255)
                                                                                           ggAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCtg                                                          <  2:989490/150‑1 (MQ=255)
                                                                                               ggcggcTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCtgcatg                                                      >  2:1902237/1‑150 (MQ=255)
                                                                                                  ggcTACGCTGTGGAAGAAATCGACATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGccc                                                                                             >  1:294532/1‑108 (MQ=255)
                                                                                                  ggcTACGCTGTGGAAGAAATCGACATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGccc                                                                                             <  2:294532/108‑1 (MQ=255)
                                                                                                   gcTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGc                                                                                               >  2:1291551/1‑105 (MQ=255)
                                                                                                   gcTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGc                                                                                               <  1:1291551/105‑1 (MQ=255)
                                                                                                   gcTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCTGCATGcaca                                                  >  2:1794381/1‑150 (MQ=255)
                                                                                                     tACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCTGCATGCACAAg                                                >  2:1005877/1‑150 (MQ=255)
                                                                                                        gCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAg                                                                  >  1:1636625/1‑129 (MQ=255)
                                                                                                        gCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAg                                                                  <  2:1636625/129‑1 (MQ=255)
                                                                                                           gtgGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCAGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCTGCATGCACAAGCCGagc                                          <  1:919987/150‑1 (MQ=255)
                                                                                                            tgGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCTGCATGCACAAGCCGagca                                         <  1:743493/150‑1 (MQ=255)
                                                                                                                 gaaATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCACCGTCAGGATCTCTACAATTGTCCTGCTTAAACGTCGGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCTGCATGCACAAGCCGAGCAGCGCa                                    <  2:1188354/150‑1 (MQ=255)
                                                                                                                        gCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGt                                                                                                                                       >  2:625670/1‑45 (MQ=255)
                                                                                                                        gCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGt                                                                                                                                       <  1:625670/45‑1 (MQ=255)
                                                                                                                            cAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCCGAGGCGCCGAACACTGTCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCTGCATGCACAAGCCGAGCAGCGCACCCTGCGTGta                         <  2:1020570/150‑1 (MQ=255)
                                                                                                                             aaCGCGGTCAATGTGCTGGAAGGCAAA‑‑CAACGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCTGCATGCACAAGCCGAGCAGCGCACCCTGCGTGtat                        >  1:1020570/1‑150 (MQ=255)
                                                                                                                               cgctGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCATTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGcc                                                                                              <  1:174657/78‑1 (MQ=255)
                                                                                                                               cgcgGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGcc                                                                                              >  2:174657/1‑78 (MQ=255)
                                                                                                                               cgcgGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCTGCATGCACAAGCCGAGCAGCGCACCCTGCGTGtatat                      >  1:616665/1‑150 (MQ=255)
                                                                                                                                          tgCTGGAAGGTTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGTCCTGCATGCACAAGCCGAGCAGCGCACCCTGCGTGTATATAACTGGTTCGa           <  2:886894/150‑1 (MQ=255)
                                                                                                                                             tGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCTGCATGCACAAGCCGAGCAGCGCACCCTGCGTGTATATAACTGGTTCGATTa        <  1:1005877/150‑1 (MQ=255)
                                                                                                                                               gAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCAGCGCGTACAAGCCGAGCAGCGCACCCTGCGTGTATATAACTGGTTCGATTACa      >  1:1574051/1‑150 (MQ=255)
                                                                                                                                                  ggCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCTGCATGCACAAGc                                               >  1:950706/1‑106 (MQ=255)
                                                                                                                                                  ggCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCTGCATGCACAAGc                                               <  2:950706/106‑1 (MQ=255)
                                                                                                                                                   gCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCTGCATGCACAAGCCGAGCAGCGCACCCTGCGTGTATATAACTGGTTCGATTACATCAc  >  2:1750261/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
ACCCCATCTCGCAGTTCAAGCTGGACAGCCCGGATTACCTGGAAATGGGCAAGCGCATCGGCCAACTCGGCAAGCCGACCCTGTTCGTGATGGAAGGCGGCTACGCTGTGGAAGAAATCGGCATCAACGCGGTCAATGTGCTGGAAGGCTTC‑‑CAGCGCGCCCAGTCAGGAGCCCGATAATGGTCCGACTCAAACGTCTGCTCGCCCCGTTCATCGCTGCCACCCTGTTCGCAGGCGCCCTGCATGCACAAGCCGAGCAGCGCACCCTGCGTGTATATAACTGGTTCGATTACATCAC  >  NC_002947/6088700‑6088996

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: