Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,792,273 G→A 6.5% P372S (CCT→TCT)  gltD ← glutamate synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,792,2730GA 6.5% 175.6 / 2.6 61P372S (CCT→TCT) gltDglutamate synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (25/32);  new base A (2/2);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

ACCTTGGCTGTCCAGCTGAATGCCATGCTGTTCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCTTAACCTCTTTGCGCGAACCCGGCATGTTGGCTTCG  >  NC_002947/5792127‑5792418
                                                                                                                                                  |                                                                                                                                                 
aCCTTGGCTGTCCAGCTGAATGCCATGCTGTTCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTc                                                                                                                                                <  2:1370596/150‑1 (MQ=255)
   ttGGCTGTCCAGCTGAATGCCATGCTGTTCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGACTCGcc                                                                                                                                             >  2:433189/1‑150 (MQ=255)
          tCCAGCTGAATGCCATGCTGTTCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGa                                                                                                                                      <  2:388991/150‑1 (MQ=255)
               cTGAATGCCATGCTGTTCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGTCTTCAGGATCTCTTCGGCGCCCGGGACGGGCACAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGtctc                                                                                                                                 >  1:1386454/1‑150 (MQ=255)
               cTGAATGCCATGCTGTTCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGtctc                                                                                                                                 <  1:221237/150‑1 (MQ=255)
                tGAATGCCATGCTGTTCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGACTCGCCTGGGCGAgcctga                                                                                                                                <  1:433189/150‑6 (MQ=255)
                  aaTGCCATGCTGTTCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGAc                                                                                                                              <  1:504762/150‑1 (MQ=255)
                       cATGCTGTTCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACtt                                                                                                                         <  1:944340/150‑1 (MQ=255)
                       cATGCTGTTCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACtt                                                                                                                         >  1:1222183/1‑150 (MQ=255)
                       cATGCTGTTCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACtt                                                                                                                         >  2:906284/1‑150 (MQ=255)
                            tGTTCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACg                                                                                                                    >  2:990288/1‑150 (MQ=255)
                              ttCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGAc                                                                                                                     <  1:926337/147‑1 (MQ=255)
                              ttCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGAc                                                                                                                     >  2:926337/1‑147 (MQ=255)
                                 aaaCCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTc                                                                                                               <  2:531149/150‑1 (MQ=255)
                                   aCCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGa                                                                                                             >  1:920279/1‑150 (MQ=255)
                                    ccACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGAc                                                                                                            >  1:1210300/1‑150 (MQ=255)
                                     cACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGAcc                                                                                                           <  2:1210300/150‑1 (MQ=255)
                                             cGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGtc                                                                                                                                   <  2:344958/118‑1 (MQ=255)
                                             cGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGtc                                                                                                                                   >  1:344958/1‑118 (MQ=255)
                                             cGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTc                                                                                                   <  2:89556/150‑1 (MQ=255)
                                                       cGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCAcaatg                                                                                         <  2:423107/150‑1 (MQ=255)
                                                        ggAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCAcaatgg                                                                                        <  2:1407497/150‑1 (MQ=255)
                                                                aaGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGc                                                                                <  2:1222183/150‑1 (MQ=255)
                                                                   tcGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGa                                                                             <  1:1014037/149‑1 (MQ=255)
                                                                     gATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTAGCCAACCTTGGCACCGGGCTGata                                                                           >  2:1014037/1‑148 (MQ=255)
                                                                      aTTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACgg                                                                          <  1:990288/150‑1 (MQ=255)
                                                                      aTTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACgg                                                                          >  1:1060348/1‑150 (MQ=255)
                                                                          ccACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGAGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGt                                                                      >  2:346424/1‑150 (MQ=255)
                                                                           cACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGa                                                                                                                               <  1:1140217/92‑1 (MQ=255)
                                                                           cACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGa                                                                                                                               >  2:1140217/1‑92 (MQ=255)
                                                                           cACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTa                                                                     <  2:688910/150‑1 (MQ=255)
                                                                           cACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGACTCGCCTTGGCGAgc                                                                                                                                    <  2:1417450/87‑2 (MQ=255)
                                                                           cACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGACTCGCCTGGGCGAgc                                                                                                                                    >  1:1417450/1‑86 (MQ=255)
                                                                             cGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTaga                                                                   <  2:1193263/150‑1 (MQ=255)
                                                                              ggCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTagag                                                                  <  2:1060348/150‑1 (MQ=255)
                                                                              ggCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTagag                                                                  <  1:1420459/150‑1 (MQ=255)
                                                                               gCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAgg                                                                 <  2:920279/150‑1 (MQ=255)
                                                                                     gCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTc                                                           <  1:1250664/150‑1 (MQ=255)
                                                                                         gCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGTGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGc                                                       <  1:346424/150‑1 (MQ=255)
                                                                                           aGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGccc                                                     <  1:906284/150‑1 (MQ=255)
                                                                                                    tCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTg                                            <  2:2708/150‑1 (MQ=255)
                                                                                                      gtaGCCCGGCATTGGCTATAATCTGCGGCATCCACGGGCATCACGCTCGCCCAGGCGAGTCTCGACATATTTTACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGGTGTAGAGGAACTTCACGCCCTCTTCCTTGGc                                          <  1:1380785/148‑1 (MQ=255)
                                                                                                          cccGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTc                                      >  1:457601/1‑150 (MQ=255)
                                                                                                           ccGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCt                                     >  1:468415/1‑150 (MQ=255)
                                                                                                            cGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCtt                                    >  2:1039181/1‑150 (MQ=255)
                                                                                                              ggATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCTTaa                                  <  2:674305/150‑1 (MQ=255)
                                                                                                                    gCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCTTAACCTCtt                            >  1:323341/1‑150 (MQ=255)
                                                                                                                      tCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCAcaatgg                                                                                        >  2:343044/1‑88 (MQ=255)
                                                                                                                      tCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCAcaatgg                                                                                        <  1:343044/88‑1 (MQ=255)
                                                                                                                      tCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCTTAACCTCTTTg                          >  1:702237/1‑150 (MQ=255)
                                                                                                                                 gcggcCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCTTAACCTCTTTGCGCGAACCCgg               >  1:159388/1‑150 (MQ=255)
                                                                                                                                         ggggATCAGGAACGCCCAAGCGACTCTCGACCACCTTGACGCCGTCGACCTTGTCCTCGCCCACAATAACACTCGGCTGACGCTTGAACAGGAAAAAAACACGCTCTTCCTTAGCGTTCGTAACCTCTTTGCGCCAACCCGGCATGTTgg       >  2:87981/1‑150 (MQ=255)
                                                                                                                                           gCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCTTAACCTCTTTGCGCGAACCCGGCATGTTGGCt     >  2:277177/1‑150 (MQ=255)
                                                                                                                                             aTCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCTGTCGTTACCGCTGGGCGCGAACCCGGCATGTTGGCTTc   >  2:233743/1‑150 (MQ=255)
                                                                                                                                             aTCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCTTAACCTCTTTGCGCGAACCCGGCATGTTGGCTTc   <  1:277177/150‑1 (MQ=255)
                                                                                                                                             aTCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCTTAACCTCTTTGCGCGAACCCGGCATGTTGGCTTc   >  1:208607/1‑150 (MQ=255)
                                                                                                                                             aTCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCTTAACCTCTTTGCGCGAACCCGGCATGTTGGCTTc   <  2:580808/150‑1 (MQ=255)
                                                                                                                                             aTCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCTTAACCTCTTTGCGCGAACCCGGCATGTTGGCTTc   <  2:208607/150‑1 (MQ=255)
                                                                                                                                              tCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCTTAACCTCTTTGCGCGAACCCGGCATGTTGGCTTCg  >  2:780025/1‑150 (MQ=255)
                                                                                                                                               cAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCTTAACCTCTTTGCGCGAACCCGGCATGTTGGCTTc   <  2:345905/148‑1 (MQ=255)
                                                                                                                                               cAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCTTAACCTCTTTGCGCGAACCCGGCATGTTGGCTTc   >  1:345905/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACCTTGGCTGTCCAGCTGAATGCCATGCTGTTCAAACCACGGTGCCGGGCTCGGGCGGAAGCCGAAGGCGATTACCACGGCGTCTGCCGGCAGGATCTCTTCGGAGCCCGGGATTGGCTCAGGGCTGCGGCGGCCACGGGCATCAGGCTCGCCCAGGCGAGTCTCGACCACTTTGACGCCTTCGACCTTGTCCTCGCCCACAATGGCAATGGGCTGACGGTTGTAGAGGAACTTCACGCCCTCTTCCTTGGCGTTCTTAACCTCTTTGCGCGAACCCGGCATGTTGGCTTCG  >  NC_002947/5792127‑5792418

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