Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,201,086 A→C 6.4% intergenic (‑49/+66) PP_3690 ← / ← PP_3691 hypothetical protein/DNA helicase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,201,0860AC 6.4% 196.1 / 2.9 62intergenic (‑49/+66)PP_3690/PP_3691hypothetical protein/DNA helicase‑related protein
Reads supporting (aligned to +/- strand):  ref base A (27/31);  new base C (2/2);  total (29/33)
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.80e-01

GCATTTGTGTGCACCGCTGGGTCACCGTGATCCGACTAAAATGCAGAACCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTGGGCAAGCTTGCTGACCCCAGCTTCTTTGG  >  NC_002947/4200938‑4201231
                                                                                                                                                     |                                                                                                                                                   
gCATTTGTGTGCACCGCTGGGTCACCGTGATCCGACTAAAATGCAGAACCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTa                                                                                                                                                     <  4:42238/149‑1 (MQ=255)
 cATTTGTGTGCACCGCTGGGTCACCGTGATCCGACTAAAATGCAGAACCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAg                                                                                                                                                    <  3:203669/149‑1 (MQ=255)
 cATTTGTGTGCACCGCTGGGTCACCGTGATCCGACTAAAATGCAGAACCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGATTAATCTGAGCGATGCGCAGCTAgg                                                                                                                                                   <  3:25684/149‑1 (MQ=255)
       tgtgCACCGCTGGGTCACCGTGATCCGACTAAAATGCAGAACCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTcaac                                                                                                                                              <  3:257243/149‑1 (MQ=255)
            aCCGCTGGGTCACCGTGATCCGACTAAAATGCAGAACCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGa                                                                                                                                         <  2:243681/149‑1 (MQ=255)
             ccGCTGGGTCACGGTTATCCGACTTAAATTCAGAAGGATCTCACGTTTGATGAAGGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATATGAGCGATGCGCATCTATGTCAACAACGaa                                                                                                                                        <  3:201988/149‑1 (MQ=255)
             ccGCTGGGTCACCGTGATCCGACTAAAATGCAGAACCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGaa                                                                                                                                        <  5:133162/149‑1 (MQ=255)
             ccGCTGGGTCACCGTGATCCGACTAAAATGCAGAACCAACTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGaa                                                                                                                                        <  5:81290/149‑1 (MQ=255)
                cTGGGTCACCGTGATCCGACTAAAATGCAGAACCATCTCACCTTTGATGAAAGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGc                                                                                                                                     <  8:225788/149‑1 (MQ=255)
                             aTCCGACTAAAATGCAGAACCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTgg                                                                                                                        <  5:102071/149‑1 (MQ=255)
                                 gACTAAAATGCAGAACCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATg                                                                                                                     >  7:23196/1‑148 (MQ=255)
                                 gACTAAAATGCAGAACCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATg                                                                                                                     >  2:202283/1‑148 (MQ=255)
                                 gACTAAAATGCAGAACCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATAAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATg                                                                                                                     >  1:113438/1‑148 (MQ=255)
                                           cAGAACCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGgc                                                                                                           <  8:23196/148‑1 (MQ=255)
                                             gAACCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCATGGATGAAAAAAGGGCGCt                                                                                                        >  5:91468/1‑149 (MQ=255)
                                             gAACCATCTCACCTTTGATGAAAGCCTCAACGTACACACCATGAGCTTAGGTCATTCTGAGATCCTAGGATGGTGGGGGCATTAATCTGAGCGATGCGAAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGCCGCt                                                                                                        >  5:261976/1‑149 (MQ=255)
                                              aaCCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTg                                                                                                       >  2:154758/1‑149 (MQ=255)
                                                ccATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGt                                                                                                     <  1:154758/149‑1 (MQ=255)
                                                 cctctcACCGTTGAGTCACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTa                                                                                                    <  8:219023/147‑1 (MQ=255)
                                                     tcACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAgc                                                                                                >  6:277501/1‑149 (MQ=255)
                                                       aCCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTaaa                                                                                                  <  7:145993/145‑1 (MQ=255)
                                                         cTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGccc                                                                                           >  7:225951/1‑149 (MQ=255)
                                                              atgaACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGa                                                                                      >  5:181288/1‑149 (MQ=255)
                                                                gaACGCCTCAACGTTCACACCATGAGCTTTGGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAgg                                                                                                                                      >  1:160817/1‑101 (MQ=255)
                                                                gaACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAgg                                                                                                                                      <  2:160817/100‑1 (MQ=255)
                                                                  aCGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTAt                                                                                  >  3:116666/1‑149 (MQ=255)
                                                                   cGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATc                                                                                 <  8:225951/149‑1 (MQ=255)
                                                                    gCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGt                                                                                                     <  1:270148/129‑1 (MQ=255)
                                                                    gCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGt                                                                                                     >  2:270148/1‑129 (MQ=255)
                                                                         aaCGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTgg                                                                                                                        >  7:163163/1‑105 (MQ=255)
                                                                         aaCGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTgg                                                                                                                        <  8:163163/105‑1 (MQ=255)
                                                                         aaCGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTg                                                                                                       >  7:8601/1‑122 (MQ=255)
                                                                         aaCGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTg                                                                                                       <  8:8601/122‑1 (MQ=255)
                                                                               cacaCCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAAt                                                                      <  5:277501/148‑1 (MQ=255)
                                                                                acacCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGTGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTgg                                                                                                      >  6:212759/1‑116 (MQ=255)
                                                                                acacCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTgg                                                                                                      <  5:212759/116‑1 (MQ=255)
                                                                                  acCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGcc                                                                  >  3:109957/1‑149 (MQ=255)
                                                                                            ttGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCttt                                                        >  8:102604/1‑149 (MQ=255)
                                                                                                gTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCttttt                                                      >  4:107869/1‑148 (MQ=255)
                                                                                                 tCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTc                                                     >  4:220798/1‑148 (MQ=255)
                                                                                                             ccTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGc                                        >  6:44286/1‑149 (MQ=255)
                                                                                                              cTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGc                                        <  1:202283/149‑1 (MQ=255)
                                                                                                                 ggTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAg                                   >  1:146355/1‑149 (MQ=255)
                                                                                                                      gtggGGGGATTAATCTGAGCGATGCGCAGCTCGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGAc                                          >  4:217030/1‑138 (MQ=255)
                                                                                                                      gtggGGGGATTAATCTGAGCGATGCGCAGCTCGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGAc                                          >  4:213280/1‑138 (MQ=255)
                                                                                                                      gtggGGGGATTAATCTGAGCGATGCGCAGCTCGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGAc                                          <  3:217030/138‑1 (MQ=255)
                                                                                                                      gtggGGGGATTAATCTGAGCGATGCGCAGCTCGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGAc                                          <  3:213280/138‑1 (MQ=255)
                                                                                                                      gtggGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTgg                              <  6:181288/149‑1 (MQ=255)
                                                                                                                      gtggGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTg                               <  3:220798/148‑1 (MQ=255)
                                                                                                                           gggATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTGGGCAAg                         >  1:255396/1‑149 (MQ=255)
                                                                                                                                  aTCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTggg                             <  1:149689/138‑1 (MQ=255)
                                                                                                                                  aTCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTggg                             >  2:149689/1‑138 (MQ=255)
                                                                                                                                     tGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTGGGCAAGCTTGCTGAcc               >  6:75693/1‑149 (MQ=255)
                                                                                                                                     tGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTGGGCAAGCTTGCTGAcc               >  2:172898/1‑149 (MQ=255)
                                                                                                                                     tGAGCGATGCGAAGCAAGCTCAACAACGAAGGCAGAGCTATCCTGGAAGAAAAAAGCGCGCTGGTAAAGACCCCCTCGAGTAACACTCATAAAATCTCCGCATGTCTATTTCAATACGCTGACTCCAAGCCTGGCCAAGCTTGCTGAcg               >  6:185360/1‑148 (MQ=255)
                                                                                                                                       aGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTGGGCAAGCTTGCTGa                 <  4:109957/145‑1 (MQ=255)
                                                                                                                                                cAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTGGGCAAGCTTGCTGACCCCAGCTTCttt    <  2:255396/149‑1 (MQ=255)
                                                                                                                                                 aGCTAGGTCAACAACGAAGGCAGATCTCTCCTGGATGAAAAAAGTGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCTGATGGCTTTTTCAATACGCTGACGCCGAGCTTGGGCAAGCTTGCTGACCCCAGCTTCTTTg   <  5:185360/149‑1 (MQ=255)
                                                                                                                                                 aGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTGGGCAAGCTTGCTGACCCCAGCTTCTTTg   <  7:97211/149‑1 (MQ=255)
                                                                                                                                                 aGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTGGGCAAGCTTGCTGACCCCAGCTTCTTTg   >  8:97211/1‑149 (MQ=255)
                                                                                                                                                  gCTATGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTGGGCAAGCTTGCTGACCCCAGCTTCTTTgg  <  6:91468/149‑1 (MQ=255)
                                                                                                                                                  gCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAAGCGTGGGCAAGCTTGCTGACCCCAGCTTCTTTg   <  1:172898/149‑1 (MQ=255)
                                                                                                                                                  gCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTGGGCAAGCTTGCTGACCCCAGCTTCTTTgg  <  6:207960/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCATTTGTGTGCACCGCTGGGTCACCGTGATCCGACTAAAATGCAGAACCATCTCACCTTTGATGAACGCCTCAACGTTCACACCATGAGCTTTGGTCATTCTGAGTTCCTAGGTTGGTGGGGGGATTAATCTGAGCGATGCGCAGCTAGGTCAACAACGAAGGCAGAGCTCTCCTGGATGAAAAAAGGGCGCTGGTAAAGCGCCCCTCGAGTATCACTCATCAAATCGCCGGATGGCTTTTTCAATACGCTGACGCCCAGCGTGGGCAAGCTTGCTGACCCCAGCTTCTTTGG  >  NC_002947/4200938‑4201231

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