Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 774,358 T→C 6.9% I322T (ATC→ACC)  PP_0665 → glyceraldehyde‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947774,3580TC 6.9% 168.2 / 2.1 58I322T (ATC→ACC) PP_0665glyceraldehyde‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (25/29);  new base C (2/2);  total (27/31)
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.12e-01

GGCGGCCGGGCAGAACTGCATCGGCACCCAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACAGTGCTCGAAAGCGTTGACCCCGTCAGCCGGCTGT  >  NC_002947/774212‑774501
                                                                                                                                                  |                                                                                                                                               
ggcggcCGGGCAGAACTGCATCGGCACCCAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATc                                                                                                                                                <  3:242240/148‑1 (MQ=255)
 gcggcCGGGCAGAACTGCATCGGCACCCAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCAc                                                                                                                                              >  6:124891/1‑149 (MQ=255)
      cGGGGAGAACTGCATCGGCACCCTGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAag                                                                                                                                         <  7:88974/149‑1 (MQ=255)
      cGGGCAGAACTGCATCGGCACCCAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAag                                                                                                                                         >  6:61471/1‑149 (MQ=255)
         gCAGAACTGCATCGGCACCCAGCGCATTCTCGTAGATGCCCCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAagcag                                                                                                                                      <  1:67808/149‑1 (MQ=255)
           aGAACTGCATCGGCACCCAGCGCATTCTCGCAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCt                                                                                                                                    >  4:207352/1‑149 (MQ=255)
                tGCATCGGCACCCAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGCTCTGGGGCGATCCGGGCCGGCGCGAACCCGATATGGGGCCGATCATCACCGAAGCAGCTGAACg                                                                                                                               >  8:141404/1‑149 (MQ=255)
                  cATCGGCACCCAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGc                                                                                                                             <  1:260903/149‑1 (MQ=255)
                   aTCGGCACCCAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAag                                                                                                                                         <  4:8397/136‑1 (MQ=255)
                   aTCGGCACCCAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAag                                                                                                                                         >  3:8397/1‑136 (MQ=255)
                    tCGGCACCCAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACgg                                                                                                                              >  1:82742/1‑146 (MQ=255)
                       gCACCCAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGACCACCGAAGCAGCTGCACgg                                                                                                                              >  1:258027/1‑143 (MQ=255)
                          cccAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCa                                                                                                                                 >  4:103329/1‑137 (MQ=255)
                          cccAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCa                                                                                                                                 <  3:103329/137‑1 (MQ=255)
                            cAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAg                                                                                                                    >  1:321844/1‑148 (MQ=255)
                            cAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATAAGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGc                                                                                                                   >  5:259151/1‑149 (MQ=255)
                             aGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCa                                                                                                                  >  4:246705/1‑149 (MQ=255)
                             aGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCa                                                                                                                  >  4:321215/1‑149 (MQ=255)
                             aGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCa                                                                                                                  >  6:109041/1‑149 (MQ=255)
                              gcgcATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCAAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGACCACCGAAGCAGCTGCACGGCAGATCGAAGCAc                                                                                                                 >  2:86141/1‑149 (MQ=255)
                                  aTTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCAcgc                                                                                                               <  2:239673/147‑1 (MQ=255)
                                  aTTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCAcgc                                                                                                               >  1:239673/1‑147 (MQ=255)
                                  aTTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGt                                                                                                             <  5:124891/149‑1 (MQ=255)
                                           gATGCCTAGTTCTACGAGGCGTTTCTAGAGCGCTTCGTGGCATTGACCCAGGCTATTGTCTTGGGCGATCCGGGGCGGCGCGAAACGGATATGGGGGCGATGATCACCGAAGCAGCTGCTCGGCAGATCGAAGCACGCGTGAACCAgg                                                                                                     <  5:55835/148‑1 (MQ=255)
                                              cgaggcaTCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGa                                                                                                  <  3:4386/142‑1 (MQ=255)
                                               ccTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGAtt                                                                                                >  1:34808/1‑149 (MQ=255)
                                               ccTCGATCTACGAGGCGTTTCGAGAGCGCTNCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGAtt                                                                                                >  1:29641/1‑149 (MQ=255)
                                                cTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTg                                                                                               <  5:61471/149‑1 (MQ=255)
                                                   gATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGa                                                                                              <  3:246705/147‑1 (MQ=255)
                                                        aCGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGt                                                                                       >  4:215619/1‑149 (MQ=255)
                                                         cGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAggg                                                                                           <  2:187615/144‑1 (MQ=255)
                                                         cGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGt                                                                                       <  3:321215/148‑1 (MQ=255)
                                                                  tCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTg                                                                              <  2:321844/148‑1 (MQ=255)
                                                                          gcTTCGTGGTATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGc                                                                     <  3:113988/149‑1 (MQ=255)
                                                                          gcTTCGTGGAATTGACCCAGGCGATGGTCTGGGGCGAACCGGGACGGCGCGAAACCGATATGGGGCAGATGATCACCGAAGCAGCTGAACGGAAGATAGAAGCACGCGTGAACCAGGCGATGGAGGGCGGTGCCCGCCTGCTAggcggc                                                                     >  4:243273/1‑149 (MQ=255)
                                                                                  gCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGc                                                              <  2:29641/148‑1 (MQ=255)
                                                                                  gCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGc                                                              <  2:34808/148‑1 (MQ=255)
                                                                                     ttGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCa                                                          >  5:19853/1‑149 (MQ=255)
                                                                                      tGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCAc                                                         >  3:310367/1‑149 (MQ=255)
                                                                                        aCCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACgg                                                       >  8:153273/1‑149 (MQ=255)
                                                                                          ccAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCg                                                                       >  7:90925/1‑131 (MQ=255)
                                                                                          ccAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCg                                                                       <  8:90925/131‑1 (MQ=255)
                                                                                          ccAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACgggg                                                     >  1:115999/1‑149 (MQ=255)
                                                                                                gATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCt                                               <  3:215619/149‑1 (MQ=255)
                                                                                                    gTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATgc                                           <  5:109041/149‑1 (MQ=255)
                                                                                                       tGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGcc                                        >  6:268228/1‑149 (MQ=255)
                                                                                                         gggCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGa                                      <  6:19853/149‑1 (MQ=255)
                                                                                                         gggCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGa                                      <  5:50981/149‑1 (MQ=255)
                                                                                                           gctaTCCGGGCCGGCGCGAAACCGATATGGGGACGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTTAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACa                                    <  3:243273/146‑1 (MQ=255)
                                                                                                               tCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACAGTGc                                >  1:173576/1‑149 (MQ=255)
                                                                                                                     ccGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACAGTGCTCGaaa                          <  3:207352/149‑1 (MQ=255)
                                                                                                                      cGGCGCGAAACCGATATGGGGCCGATGACCACCGAAGGAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACAGTGCTCGAAAg                         <  2:258027/149‑1 (MQ=255)
                                                                                                                            gAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACAGTGCTCGAAAGCGTTGa                   <  2:173576/149‑1 (MQ=255)
                                                                                                                            gAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACACTGCTCGAAAGCGTTGa                   >  1:339011/1‑149 (MQ=255)
                                                                                                                                ccGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACAGTGCTCGAAAGCGTTGAcccc               <  7:153273/149‑1 (MQ=255)
                                                                                                                                  gATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACAGTGCTCGAAAGCGTTGACCCCGt             <  7:141404/149‑1 (MQ=255)
                                                                                                                                          gCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACAGTGCTCGAAAGCGTTGACCCCGTCAGCCGGc     <  2:82742/149‑1 (MQ=255)
                                                                                                                                           ccGATGACCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTTAACTAGGCGATTGAGGGCTGTGCCCGCCTTCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACAGTGCTCGAAAGCGTTGACCCCGTCAGCCGGCt    <  1:86141/149‑1 (MQ=255)
                                                                                                                                             gatgatCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACAGTGCTCGAAAGCGTTGACCCCGTCAGCCGGCtgt  >  4:291985/1‑149 (MQ=255)
                                                                                                                                             gatgatCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACAGTGCTCGAAAGCGTTGACCCCGTCAGCCGGCtgt  <  2:115999/149‑1 (MQ=255)
                                                                                                                                             gatgatCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACAGTGCTCGAAAGCGTTGACCCCGTCAGCCGGCtgt  >  4:298298/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GGCGGCCGGGCAGAACTGCATCGGCACCCAGCGCATTCTCGTAGATGCCTCGATCTACGAGGCGTTTCGAGAGCGCTTCGTGGCATTGACCCAGGCGATGGTCTGGGGCGATCCGGGCCGGCGCGAAACCGATATGGGGCCGATGATCACCGAAGCAGCTGCACGGCAGATCGAAGCACGCGTGAACCAGGCGATTGAGGGCGGTGCCCGCCTGCTCTGCGGCCACCAGCGTCACGGGGCGGTCTATGCGCCGACAGTGCTCGAAAGCGTTGACCCCGTCAGCCGGCTGT  >  NC_002947/774212‑774501

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