Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,411,873 G→A 6.5% A21V (GCT→GTT)  PP_1232 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,411,8730GA 6.5% 182.8 / 2.0 61A21V (GCT→GTT) PP_1232hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/31);  new base A (2/2);  total (28/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.64e-01

GTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTCACACGTTGACGACCAGTGGCGCTATCCTAGCCA  >  NC_002947/1411729‑1412012
                                                                                                                                                |                                                                                                                                           
gTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATg                                                                                                                                         >  4:10072/1‑149 (MQ=255)
 tttCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATgg                                                                                                                                        >  6:75581/1‑149 (MQ=255)
  ttCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGc                                                                                                                                       >  1:73657/1‑149 (MQ=255)
           gTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAgg                                                                                                                               <  7:256463/148‑1 (MQ=255)
            tCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAACATGGCCGGGAAGGGc                                                                                                                             >  3:66465/1‑149 (MQ=255)
            tCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAACATGGCCGGGAAGGGc                                                                                                                             >  3:61884/1‑149 (MQ=255)
              cTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGcc                                                                                                                            >  7:53120/1‑148 (MQ=255)
              cTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGcc                                                                                                                            <  6:87583/148‑1 (MQ=255)
                 tagttgtgGGTCATTGATCTTATTGAACGTGCCGCGGAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATGGCAGTAGTGGCATCACGCTGTGATGGGAGGTCGTGTGCTTGGGGGGTAAGGGGCATc                                                                                                                         <  6:283342/144‑1 (MQ=255)
                   gCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATcagg                                                                                                                      <  2:333899/149‑1 (MQ=255)
                   gCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATcagg                                                                                                                      <  3:88719/149‑1 (MQ=255)
                   gCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATcagg                                                                                                                      <  2:63941/149‑1 (MQ=255)
                     tgtgGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATcaggca                                                                                                                    <  2:336694/149‑1 (MQ=255)
                          gTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAgg                                                                                                                               >  7:137345/1‑133 (MQ=255)
                          gTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAgg                                                                                                                               <  8:137345/133‑1 (MQ=255)
                           tCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAg                                                                                                              >  7:133597/1‑149 (MQ=255)
                           tCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAg                                                                                                              >  5:112834/1‑149 (MQ=255)
                             aTTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATcaggcag                                                                                                                   <  7:221940/142‑1 (MQ=255)
                             aTTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATcaggcag                                                                                                                   >  8:221940/1‑142 (MQ=255)
                                 aGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTcag                                                                                                        <  8:53120/149‑1 (MQ=255)
                                  gCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTcagc                                                                                                       <  5:72660/149‑1 (MQ=255)
                                      aTTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAggg                                                                                                   >  1:111942/1‑149 (MQ=255)
                                          aCCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGG‑CAGCGTCAGCAGGGTAgg                                                                                               <  5:75581/148‑1 (MQ=255)
                                            cTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAACATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGgc                                                                                             <  4:61884/149‑1 (MQ=255)
                                            cTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAACATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGgc                                                                                             <  4:66465/149‑1 (MQ=255)
                                             tGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGATGGTGCTCCTGTTGCGGCGAGACGATAGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGCTAGGgcg                                                                                            >  2:252594/1‑149 (MQ=255)
                                               gCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAgg                                                                                                    >  1:240421/1‑139 (MQ=255)
                                               gCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAgg                                                                                                    <  2:240421/139‑1 (MQ=255)
                                                 cgcgCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAgg                                                                                               >  8:152248/1‑142 (MQ=255)
                                                 cgcgCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAgg                                                                                               <  7:152248/142‑1 (MQ=255)
                                                 cgcgCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGt                                                                                        <  5:106918/149‑1 (MQ=255)
                                                   cgcAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATcaggca                                                                                                                    >  2:83742/1‑119 (MQ=255)
                                                   cgcAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATcaggca                                                                                                                    <  1:83742/119‑1 (MQ=255)
                                                   cgcAGCAGGCGGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAg                                                                                      <  6:200630/149‑1 (MQ=255)
                                                                    ggCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTc                                                                     <  6:112834/149‑1 (MQ=255)
                                                                     gCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTcat                                                                                >  4:170825/1‑137 (MQ=255)
                                                                     gCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTcat                                                                                <  3:170825/137‑1 (MQ=255)
                                                                      cACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTcgc                                                                   <  3:10072/149‑1 (MQ=255)
                                                                            ccAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGTCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGa                                                             <  5:145912/149‑1 (MQ=255)
                                                                              aGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATc                                                           <  8:133597/149‑1 (MQ=255)
                                                                                      cTCCTGTTGCGGCTAGACTACCGCAGAACTGGCATCACCCCGTGTTGGCTGGTCGTCTGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCTGGGTAGGGCGCAGTAGATTCTTGCATGAAGCTCGCGTCGTATCaaacaaa                                                    <  2:190370/148‑5 (MQ=255)
                                                                                       tCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATAcc                                                  >  1:12617/1‑149 (MQ=255)
                                                                                       tCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATAcc                                                  >  5:207737/1‑149 (MQ=255)
                                                                                           gTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTcatgc                                                                              >  6:157340/1‑117 (MQ=255)
                                                                                                   gagaCGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAg                                          <  1:134665/145‑1 (MQ=255)
                                                                                                      aCGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCt                                    <  6:207737/148‑1 (MQ=255)
                                                                                                      aCGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCt                                    >  8:4837/1‑148 (MQ=255)
                                                                                                       cGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTca                                  <  2:111942/149‑1 (MQ=255)
                                                                                                        gATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTca                                  >  6:178352/1‑148 (MQ=255)
                                                                                                          tCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCGGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTCacac                               >  3:164664/1‑149 (MQ=255)
                                                                                                                 tcTGCCATCACTGGGTGATGGCAGGTAGTCAGCATGGTCGGGAAGTGCCATCAGGCAGGCCAGCGTCAGGAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGTGCTCACACGTTgac                         <  4:203950/147‑1 (MQ=255)
                                                                                                                       aTCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGAt                                                     >  3:105690/1‑114 (MQ=255)
                                                                                                                       aTCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGAt                                                     <  4:105690/114‑1 (MQ=255)
                                                                                                                        tCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTCACACGTTGACGACCAGTg                 >  3:42533/1‑149 (MQ=255)
                                                                                                                          aCCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTCACACGTTGACGACCAGTgg                <  7:223848/148‑1 (MQ=255)
                                                                                                                          aCCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGACATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTAACACGTTGACGACCAGTgtc               >  8:112589/1‑147 (MQ=255)
                                                                                                                            ccAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTCACACGTTGACGACCAGTGgcgc             >  1:136770/1‑149 (MQ=255)
                                                                                                                             cAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTCACACGTTGACGACCAGTGgcgc             >  4:61029/1‑148 (MQ=255)
                                                                                                                                   tGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCGAGCGTCAGCAGTGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAATATACCTTCTGTAGCGGGCTCACACGTTGACGACCAGTGGCGCTATCCTa      <  1:252594/149‑1 (MQ=255)
                                                                                                                                     gCAGGTAGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCATCGTCAGCAGGGTAGGGAGCAGTATATTAATTCATGAAGCTCGCGTCGGATCAAAGCTACCTACTGTAGCGGGCTCACACGTTGACGACCAGTGGCGCTATCCTagc    >  5:242447/1‑149 (MQ=255)
                                                                                                                                       aGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTCACACGTTGACGACCAGTGGCGCTATCCTagcca  <  2:73657/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                           
GTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTCACACGTTGACGACCAGTGGCGCTATCCTAGCCA  >  NC_002947/1411729‑1412012

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