Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,653,199 T→G 8.4% L362* (TTA→TGA)  PP_3218 → NtaA/SnaA/SoxA family monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,653,1990TG 8.4% 219.3 / 10.4 71L362* (TTA→TGA) PP_3218NtaA/SnaA/SoxA family monooxygenase
Reads supporting (aligned to +/- strand):  ref base T (36/29);  new base G (3/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TACGACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCGTTGCT  >  NC_002947/3653060‑3653346
                                                                                                                                            |                                                                                                                                                   
tACGACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCATTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAgg                                                                                                                                            >  8:163386/1‑149 (MQ=255)
tACGACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTAACTGACGGCCaaa                                                                                                                                            <  5:157353/149‑3 (MQ=255)
   gACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGc                                                                                                                                         >  1:9880/1‑149 (MQ=255)
   gACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGc                                                                                                                                         >  4:387509/1‑149 (MQ=255)
    aCCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGcc                                                                                                                                        <  3:387509/149‑1 (MQ=255)
              cccGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGtt                                                                                                                              >  1:258371/1‑149 (MQ=255)
                  ttCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGc                                                                                                                          <  2:258371/149‑1 (MQ=255)
                     ccGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCg                                                                                                                        >  7:57132/1‑148 (MQ=255)
                       gCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGaaa                                                                                                                     >  4:192255/1‑149 (MQ=255)
                           gagagCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGt                                                                                                                               <  6:241761/135‑1 (MQ=255)
                           gagagCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTAAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGt                                                                                                                               >  5:241761/1‑135 (MQ=255)
                                  tGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTggcggc                                                                                                           >  4:84611/1‑148 (MQ=255)
                                   ggAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTggc                                                                                                              >  8:118395/1‑144 (MQ=255)
                                   ggAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTggc                                                                                                              <  7:118395/144‑1 (MQ=255)
                                         tGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGGCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATg                                                                                                    <  4:11714/148‑1 (MQ=255)
                                         tGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATg                                                                                                    <  3:192255/148‑1 (MQ=255)
                                         tGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTGGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATg                                                                                                    <  1:299124/148‑1 (MQ=255)
                                            aCGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAg                                                                                                                                             <  1:314742/104‑1 (MQ=255)
                                            aCGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAg                                                                                                                                             >  2:314742/1‑104 (MQ=255)
                                            aCGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGcc                                                                                                 >  6:153211/1‑148 (MQ=255)
                                                          tCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTg                                                                                                                >  7:227673/1‑119 (MQ=255)
                                                          tCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTg                                                                                                                <  8:227673/119‑1 (MQ=255)
                                                            aCGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGcgcg                                                                                <  8:57132/149‑1 (MQ=255)
                                                             cGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCgg                                                                               >  1:246542/1‑149 (MQ=255)
                                                              gAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCggg                                                                              >  8:167922/1‑149 (MQ=255)
                                                                     aaGCGTATTGCCGAGCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTgg                                                                                                                                >  7:351308/1‑92 (MQ=255)
                                                                     aaGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTgg                                                                                                                                <  8:351308/92‑1 (MQ=255)
                                                                     aaGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGa                                                                       >  3:103542/1‑149 (MQ=255)
                                                                     aaGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGa                                                                       <  3:84611/149‑1 (MQ=255)
                                                                      aGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGAt                                                                      <  7:251439/149‑1 (MQ=255)
                                                                      aGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGAt                                                                      <  4:149413/149‑1 (MQ=255)
                                                                            ttGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGtt                                                                 <  7:167922/148‑1 (MQ=255)
                                                                               ccGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGATCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTc                                                                                        >  1:54572/1‑123 (MQ=255)
                                                                               ccGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTc                                                                                        <  2:54572/122‑1 (MQ=255)
                                                                               ccGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCaa                                                              >  2:85785/1‑148 (MQ=255)
                                                                                 gACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGTGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAatat                                                           >  1:23479/1‑149 (MQ=255)
                                                                                   ccAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCt                                                         <  2:246542/149‑1 (MQ=255)
                                                                                              aCAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCAt                                              >  2:106064/1‑149 (MQ=255)
                                                                                               cAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATg                                                                                                    >  2:316339/1‑94 (MQ=255)
                                                                                               cAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATg                                                                                                    <  1:316339/94‑1 (MQ=255)
                                                                                                 aGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCt                                           <  5:153211/149‑1 (MQ=255)
                                                                                                          gTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAat                                                             >  2:325910/1‑122 (MQ=255)
                                                                                                          gTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAat                                                             <  1:325910/122‑1 (MQ=255)
                                                                                                           tCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTAc                                                       <  4:367748/127‑1 (MQ=255)
                                                                                                           tCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTAc                                                       >  3:367748/1‑127 (MQ=255)
                                                                                                            cTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCcg                                >  4:177783/1‑149 (MQ=255)
                                                                                                            cTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCcg                                >  5:150304/1‑149 (MQ=255)
                                                                                                               cGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCcggcg                             >  2:21532/1‑149 (MQ=255)
                                                                                                               cGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCcggcg                             >  2:340603/1‑149 (MQ=255)
                                                                                                                     ggCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTa                        >  7:82640/1‑149 (MQ=255)
                                                                                                                      gCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTAc                       <  1:106064/149‑1 (MQ=255)
                                                                                                                        tGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTa                                                        <  8:350396/114‑1 (MQ=255)
                                                                                                                        tGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTa                                                        >  7:350396/1‑114 (MQ=255)
                                                                                                                        tGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTAcc                      >  4:7738/1‑148 (MQ=255)
                                                                                                                        tGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTAcc                      >  1:295418/1‑148 (MQ=255)
                                                                                                                          ttGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCg                     >  6:286080/1‑147 (MQ=255)
                                                                                                                            gAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCgaag                  <  3:177783/148‑1 (MQ=255)
                                                                                                                            gAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCgaag                  <  6:150304/148‑1 (MQ=255)
                                                                                                                             aGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGTGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCgaagaa                <  2:23479/149‑1 (MQ=255)
                                                                                                                              gTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCgaagaa                <  6:59746/148‑1 (MQ=255)
                                                                                                                              gTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCgaag                  >  5:59746/1‑146 (MQ=255)
                                                                                                                              ccctGGCCGTAACTGACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGaagaag               >  8:187996/4‑149 (MQ=255)
                                                                                                                                ctGGCCGTAACTGACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCcggcgg                            <  3:244551/133‑1 (MQ=255)
                                                                                                                                ctGGCCGTAACTGACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCcggcgg                            >  4:244551/2‑134 (MQ=255)
                                                                                                                                ctGGCCGTAACTGACGGCCAGGGGCCTTCCTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTg             >  2:232168/2‑149 (MQ=255)
                                                                                                                                 tGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACAGTGGCGGCGAAGATGGCCGAGTGGTTAGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGTTTTCCCGAAGAAGTgg            >  6:386102/1‑149 (MQ=255)
                                                                                                                                   gCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTa          <  1:340603/149‑1 (MQ=255)
                                                                                                                                   gCCGTAACTGACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTa          <  7:187996/149‑1 (MQ=255)
                                                                                                                                      gTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCg       >  4:18092/1‑149 (MQ=255)
                                                                                                                                         aGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGc                                                                                                  <  7:232490/55‑1 (MQ=255)
                                                                                                                                         aGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGc                                                                                                  >  8:232490/1‑55 (MQ=255)
                                                                                                                                           ttACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCGTtgct  >  2:208360/1‑149 (MQ=255)
                                                                                                                                           ttACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCGTtgct  >  6:74181/1‑149 (MQ=255)
                                                                                                                                           ttACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCGTtgct  <  2:9880/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
TACGACCTGGATGCCCCGTTCCCGCAAGAGAGCCTGGAGTTTGCACGGAGCAGCTTTTTCACGAATGCCAAGCGTATTGCCGACCAGGCGCGGGACAAGGCCTTCAGTCTGCGTCAGGCTGTTGAGTTTGGCCGTCAGTTACGGCCAGGGGCCTTCGTTGGTTCGGCCGAAACCGTGGCGGCGAAGATGGCCGAGTGGTTCGAAGCGCGGGCGCTGGATGGGTTCAATATCTACATTGGTCATCCTGGCCAGTTCCGGCGGTTTACCGAAGAAGTGGTACCGTTGCT  >  NC_002947/3653060‑3653346

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