Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,056,547 G→T 6.9% intergenic (‑208/+78) glnL ← / ← glnA sensory histidine kinase NtrB/glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,056,5470GT 6.9% 181.0 / 3.7 68intergenic (‑208/+78)glnL/glnAsensory histidine kinase NtrB/glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base T (2/3);  total (33/35)
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.83e-01

AGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTCATCTTCTT  >  NC_000913/4056407‑4056689
                                                                                                                                            |                                                                                                                                              
aGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCa                                                                                                                                              >  1:914481/1‑143 (MQ=255)
aGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCa                                                                                                                                              >  1:503581/1‑143 (MQ=255)
   aTCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACATTGAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATcc                                                                                                                                           >  1:397954/1‑143 (MQ=255)
      agagAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGAAAAAAAAATTTAATTGCTTTTGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGc                                                                                                                                        <  1:686065/143‑1 (MQ=255)
          aaTTGCCGGAGAAAAAAGCCCACGCAGAGATGGGCTACAGATGGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGa                                                                                                                                    <  1:79268/143‑1 (MQ=255)
          aaTTGACGGCGAAAAAATCCGATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGa                                                                                                                                    <  2:532131/143‑1 (MQ=255)
          aaTTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATCAGACGCATTTGCGTCGCATCCGGCGTGa                                                                                                                                    >  2:422705/1‑143 (MQ=255)
           cgttcctgagaAAAAAGCCCATGAAAAGATGGGCTACAGCTAGCTGACCCCGCTTGCGTTGGAGAACGACTGACCTGCTCCGGCCGACCGCATCCTCCGCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGAt                                                                                                                                   <  2:152454/136‑1 (MQ=255)
             tGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATGCGACGCAAATGCGTCTTATCCGGcctacacg                                                                                                                                 <  1:42411/143‑8 (MQ=37)
             tGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATGCGACGCAAATGCGTCTTATCCGGcctacacg                                                                                                                                 <  2:848941/143‑8 (MQ=37)
               acgcggaACAAAGCCCATGGATGGACGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCGCCTGCTCCGGCCTACCACAGCGTCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGc                                                                                                                               <  2:352216/138‑1 (MQ=255)
                 ggAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCtt                                                                                                                             <  1:498703/143‑1 (MQ=255)
                  gatgagaaaGCCGATGGAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACAGCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTg                                                                                                                            <  2:297797/137‑1 (MQ=255)
                   agaAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGc                                                                                                                           >  1:412560/1‑143 (MQ=255)
                     aaaaaaGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCa                                                                                                                         <  1:656122/143‑1 (MQ=255)
                     aaaaaaGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCa                                                                                                                         >  1:236621/1‑143 (MQ=255)
                     aaaaaaGCCCATGCAGAGAGGGCCTACAGATAGCTGCCAAACTTCACGTTGGAGAGCGCCTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCa                                                                                                                         <  2:300902/143‑1 (MQ=255)
                             ccATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTCGCCCGGATAAGACGCATCTGCGCCGCATCCCCCGAGATGTACCTGCCACCGAcccc                                                                                                                 >  2:611398/1‑143 (MQ=255)
                              cATGCAGAGATGGGCGACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCCGCTCCGGCCTACCACATCATAACCGCGTAGGCCGGATAAGACGCCCTTGCGGCGCATCCGGCGTGATGCGCCGGCCACCGAccccc                                                                                                                >  2:794796/1‑142 (MQ=255)
                                   agtgaTGCGCTACATCTCGCTGCCAACCTTCACGGTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGGATGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCt                                                                                                           <  2:756220/140‑1 (MQ=255)
                                         gggCTACGGATAGCGGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCATgagaga                                                                                                     <  2:382240/143‑6 (MQ=255)
                                         gggCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCaga                                                                                                     <  1:382099/143‑1 (MQ=255)
                                         gggCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCaga                                                                                                     <  2:270645/143‑1 (MQ=255)
                                                   tAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGCTCCCAGGCGTGCCAGATACAGGCGaa                                                                                           >  1:842462/1‑143 (MQ=255)
                                                     gCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCGCATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGaaaa                                                                                         <  2:832110/143‑1 (MQ=255)
                                                     gCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGaaaa                                                                                         <  1:278382/143‑1 (MQ=255)
                                                       tGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGGGACAGGCGAagagc                                                                                       >  1:105134/1‑142 (MQ=255)
                                                       tGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGAGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGt                                                                                       >  2:90613/1‑143 (MQ=255)
                                                         acagagctCACGGTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGttt                                                                                     <  1:835302/136‑1 (MQ=255)
                                                           aaaCTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATGTGCGTCCCATCCGGCGTGACGTGCTTGCCCCCGGTCCCAGGCCTGGCACAGACAGGCGAAAAGCTTcg                                                                                   >  2:51352/1‑142 (MQ=255)
                                                            agccgcACGTCGGCGAGCGGCGCACCTGCTCCGGCCTACCAAATCATCACCGTGTAGGCCGGATAAGACGCATTCCCGGCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCa                                                                                  <  1:901632/138‑1 (MQ=255)
                                                            aaCTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCa                                                                                  >  1:296721/1‑143 (MQ=255)
                                                             aCTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCAGCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCAc                                                                                 <  1:11256/143‑1 (MQ=255)
                                                              cTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAACAGTCCCCACg                                                                                >  2:405742/1‑143 (MQ=255)
                                                                tCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGc                                                                              <  1:163600/143‑1 (MQ=255)
                                                                  ccGTTGGAGAGCGACGCACCTGCGCCGGCCTACCACATCAGCACCGTGTAGGCTCGTTATTACACATTTGCGTCTCATCCGGCGTGATGTGCTTGCCACCGATCCCCTGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCaa                                                                            <  2:108394/142‑1 (MQ=255)
                                                                   cGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCCGCGAGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAACAGTTTCCACGGCAAc                                                                           >  2:765636/1‑143 (MQ=255)
                                                                     ttGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTCTCCACGGCaaatg                                                                         >  1:745307/1‑140 (MQ=255)
                                                                       ggAGAGCGACGCACCTGCTCCGGCCTACCACACCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTaaa                                                                       <  2:437860/143‑1 (MQ=255)
                                                                         agagCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAGAAGGTTCCACGGCAACTaaacc                                                                     >  2:722066/1‑141 (MQ=255)
                                                                             cgccTCGCCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACtt                                                                 <  1:241862/140‑1 (MQ=255)
                                                                              gACTCACCTGCTCCGGCCTACCACGTCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATGCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCTGTGTCAAACCAAAACTTTCCACGGCAACTAAAACACTTa                                                                <  2:924857/143‑1 (MQ=255)
                                                                              gACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGCTCCCAGGCCTGCCAGAGACAGGCGCCACGTTTCCACGGCCACTAAAACACtaa                                                                >  1:451731/1‑141 (MQ=255)
                                                                                 tCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGAc                                                             >  1:531607/1‑143 (MQ=255)
                                                                                 tCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGAGAAAACTTTCCACGGCAACTAAAACACTAAGAc                                                             >  2:580988/1‑143 (MQ=255)
                                                                                  cACCTGCTCCAGCCAATCACGTCTGCGCCGTGTAGGCCGGAGAAGACGCAGTTGCGTCGCTTCCGGCGTGATGCGCTTGCCACCGCGCCCAGGCCTGCCAGAGACAGGCGAAGAGTTTCCACGGCAACTAAAACACTTAGACg                                                            <  2:861913/143‑1 (MQ=255)
                                                                                   aCCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCCTCAGGCGTGATGTGCTTGCCACCGATCGCAGGCCTGCCCGACACAGGGGGAACGATGCCACGGCACCTGACACACTTAAACGc                                                           >  2:504330/1‑143 (MQ=255)
                                                                                           cGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACCGGCGACAAGTTTCCCCGGCAACTACAACACTTAGACGCCGTAGTAc                                                   >  2:924006/1‑143 (MQ=255)
                                                                                           cGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTAc                                                   >  1:906249/1‑143 (MQ=255)
                                                                                             gCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACCAAAACACTTAGACGCCGTAGTACAg                                                 >  2:49046/1‑143 (MQ=255)
                                                                                                       tcatcaCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGGGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGGTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAActct                                       >  1:829680/1‑143 (MQ=255)
                                                                                                          tcaCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTCAACCACTTCGCCGCCGTCGTCCAGCTCAACCTCTAcc                                    >  1:64163/1‑143 (MQ=255)
                                                                                                               gtgtAGGCCGGATCAGACGCATTTGCGCCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATg                               <  2:84053/143‑1 (MQ=255)
                                                                                                                   aGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGa                           <  1:316858/143‑1 (MQ=255)
                                                                                                                    ggCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAACGTTTCCACGGCAACTAAAACACTTAGACGCTGTCGTACAGCTCAAAGTCTACCGGATGCGGac                          >  2:713853/1‑142 (MQ=255)
                                                                                                                    ggCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAg                          <  1:191197/143‑1 (MQ=255)
                                                                                                                    cgCCGGATGCGACGCAAAGGCGTCTTATCCGGCGTGATGTGATTGCCAACGATCCCAGGCCTGAGAGAGGGGGGCGAAAGGTGTCCACGGCAACTAAAACACTTAGACGCTGTAGTAGAGCTGAAACTCTACCGGACGCGGAg                          >  1:132266/2‑143 (MQ=255)
                                                                                                                     gCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGt                         >  1:277198/1‑143 (MQ=255)
                                                                                                                     gCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGt                         <  2:800049/143‑1 (MQ=255)
                                                                                                                      ccGGATGCGACGCAAATGCGTCTTATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCAACCGGATGCGGAGTc                        >  2:224905/1‑143 (MQ=255)
                                                                                                                      ccGGATAAGTCGCAGTTGCGTCGCATACGGCGTGATGTTCTTGCCATTGATCTCAGGCCTGTCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTc                        <  1:719765/143‑1 (MQ=255)
                                                                                                                      ccGGATAAGACGCGTTTGCGTCGCATCCGGCGTGAGGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTc                        <  2:741028/143‑1 (MQ=255)
                                                                                                                              gACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATacgcac                >  1:230409/1‑143 (MQ=255)
                                                                                                                               aCGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAACCACCTAGACGCTGTAGCACAGCTCAAACCCTAACGGATCCGGAGACATacgcacg               >  2:384261/1‑143 (MQ=255)
                                                                                                                                 gCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCAGCAGACCCAGTCCGGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCAcgcg             <  2:914481/143‑1 (MQ=255)
                                                                                                                                 gCATTTGCGTCGCATCCGGCGTGATGTGCTCGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGTCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCAcgcg             <  2:503581/143‑1 (MQ=255)
                                                                                                                                 gCATTTGCGTCGCATCCGGCGTGAAGTGCTTGCCACCGATCCCAGGCCTGCCACAGACAGGCGAACAGTTTCCACGGCAACTAAAACACTCAGACGCTGAAGCACAGCCCAAACTCTACCGGATGCGGAGTCAAACGCCcgcg             >  2:443015/1‑143 (MQ=255)
                                                                                                                                       gCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTCAtc       <  1:230847/143‑1 (MQ=255)
                                                                                                                                            gCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTCATcttctt  <  2:236621/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                              
AGCATCAGAGAATTGACGGAGAAAAAAGCCCATGCAGAGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACTAAAACACTTAGACGCTGTAGTACAGCTCAAACTCTACCGGATGCGGAGTCATACGCACGCGGTCATCTTCTT  >  NC_000913/4056407‑4056689

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