Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,918,174 C→G 8.1% L191L (CTC→CTG)  PP_4328 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,918,1740CG 8.1% 187.8 / 5.5 61L191L (CTC→CTG) PP_4328hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base G (2/3);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

GCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACG  >  NC_002947/4918031‑4918320
                                                                                                                                               |                                                                                                                                                   
gCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCgg                                                                                                                                                >  2:175232/1‑149 (MQ=255)
gCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGGCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCgg                                                                                                                                                >  6:57801/1‑149 (MQ=255)
      tACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTGAGCggcacgaa                                                                                                                                          <  8:167395/149‑5 (MQ=255)
      tACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCggcagg                                                                                                                                            >  3:134677/1‑147 (MQ=255)
         cGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGa                                                                                                                                       <  2:287061/149‑1 (MQ=255)
               cTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACt                                                                                                                                 <  5:210918/149‑1 (MQ=255)
                tGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACt                                                                                                                                 >  5:414128/1‑148 (MQ=255)
                  gTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCa                                                                                                                               <  8:177598/148‑1 (MQ=255)
                  gTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAg                                                                                                                              >  1:212721/1‑149 (MQ=255)
                      gCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCAtt                                                                                                                          >  3:35764/1‑149 (MQ=255)
                       ccTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTGCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTg                                                                                                                         <  4:348430/149‑1 (MQ=255)
                        cTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGc                                                                                                                        <  3:236068/149‑1 (MQ=255)
                         tGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGAGGCAGGCAGGACCAACTCAGCATTGcc                                                                                                                       >  3:242106/1‑149 (MQ=255)
                           cTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGcc                                                                                                                       <  4:35764/147‑1 (MQ=255)
                            tCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCaa                                                                                                                    >  2:238042/1‑149 (MQ=255)
                              aaCAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAAc                                                                                                                   <  4:80458/148‑1 (MQ=255)
                                 aGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAgct                                                                                                               >  8:378539/1‑149 (MQ=255)
                                  gCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAgctg                                                                                                              <  4:134677/149‑1 (MQ=255)
                                        aggtctgTGCCACACTTAGCATCGAGAGCCCGCGCCCGCTGGTCTTTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCATCTCCGCATTGCCCAACAGCTGCTGAcc                                                                                                        <  8:169251/143‑1 (MQ=255)
                                         ggCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGAcc                                                                                                        >  1:322822/1‑148 (MQ=255)
                                          gCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGAcc                                                                                                        >  6:413322/1‑147 (MQ=255)
                                          gCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCa                                                                                                      <  5:329539/149‑1 (MQ=255)
                                            cGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAgctgc                                                                                                             >  2:392548/1‑140 (MQ=255)
                                            cGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAgctgc                                                                                                             <  1:392548/140‑1 (MQ=255)
                                                   acacTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTGAGCggcacgaa                                                                                                                                          >  8:27749/1‑100 (MQ=255)
                                                   acacTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTGAGCggcacgaa                                                                                                                                          <  7:27749/104‑5 (MQ=255)
                                                          gCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGc‑c                                                                                     >  5:185300/1‑149 (MQ=255)
                                                            aTCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTc                                                                                   <  3:404237/149‑1 (MQ=255)
                                                                    cccgcgcccgcTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCg                                                                           >  2:90520/1‑149 (MQ=255)
                                                                     ccgcgcccgcTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCgg                                                                          >  1:173726/1‑149 (MQ=255)
                                                                        cgcccgcTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAg                                                                                                     >  4:359453/1‑120 (MQ=255)
                                                                        cgcccgcTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAg                                                                                                     <  3:359453/120‑1 (MQ=255)
                                                                           ccgcTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCt                                                                    >  8:291856/1‑149 (MQ=255)
                                                                                 gTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGc                                                              >  4:321923/1‑149 (MQ=255)
                                                                                     ttGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCACCGCCCTg                                                          >  5:280077/1‑149 (MQ=255)
                                                                                       ggCAGCCTTCTCAGCGCCCTGGTTCAGTGTGATCAGTCATTGCTCTTCGTTCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCACAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCTGCCTGCTCAACGCCCTGc                                                         <  4:242106/148‑1 (MQ=255)
                                                                                       ggCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGc                                                         <  6:414128/148‑1 (MQ=255)
                                                                                           gCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAg                                                    <  3:321923/149‑1 (MQ=255)
                                                                                           gCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAg                                                    <  3:322716/149‑1 (MQ=255)
                                                                                            ccTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGc                                                   <  5:57801/149‑1 (MQ=255)
                                                                                            ccTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGc                                                   <  1:238042/149‑1 (MQ=255)
                                                                                                     cgcCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCcgc                                           >  7:321200/1‑148 (MQ=255)
                                                                                                        ccTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGAcc                                       <  2:212721/149‑1 (MQ=255)
                                                                                                        ccTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGAcc                                       >  1:295666/1‑149 (MQ=255)
                                                                                                            gTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAgg                                   <  6:185300/149‑1 (MQ=255)
                                                                                                             tGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGc                                  <  1:90520/149‑1 (MQ=255)
                                                                                                                      aTCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGc                         >  6:290574/1‑149 (MQ=255)
                                                                                                                      aTCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGc                         <  2:295666/149‑1 (MQ=255)
                                                                                                                      aTCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAACTGAACCGCCAAGCGCGCGGTACCCGGCCTGCTAAACGCCCATCAACCGCACACCCGCCACCCCGGCGTCGATgg                          >  1:330875/1‑149 (MQ=255)
                                                                                                                      aTCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCACAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCACGGAGTCGATGGc                         >  4:236612/1‑149 (MQ=255)
                                                                                                                          gTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGAcc                     >  4:239861/1‑149 (MQ=255)
                                                                                                                             aTTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTgc                  >  2:87442/1‑149 (MQ=255)
                                                                                                                                gcgcTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTCCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACGCGGCCTGCTCAACGCCCTGCAACAGCTCCCCCGCGACCCAGGCGTCGATGGCGACCTgcgcg               >  6:343791/1‑149 (MQ=255)
                                                                                                                                 cgcTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTgcgcgc              <  7:291856/149‑1 (MQ=255)
                                                                                                                                    ttCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCcg                                            <  1:172776/116‑1 (MQ=255)
                                                                                                                                    ttCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCcg                                            >  2:172776/1‑116 (MQ=255)
                                                                                                                                      cctGCCGCTGAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAgcg         >  3:249510/3‑149 (MQ=255)
                                                                                                                                          ccGCTGAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGcc       <  4:249510/147‑1 (MQ=255)
                                                                                                                                           cGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGa     >  2:286649/1‑148 (MQ=255)
                                                                                                                                           cGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGa     <  1:286649/148‑1 (MQ=255)
                                                                                                                                             cTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACg  >  3:399191/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGC‑CTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACG  >  NC_002947/4918031‑4918320

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