Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 1,927,949 A→C R297R (CGT→CGG ptrB ← protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi1,927,9490AC95.3% 285.0 / ‑1.3 86R297R (CGT→CGGptrBprotease II
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base C (42/40);  total (44/42)
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.99e-01

TGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAGGC  >  NC_000913_3_bme_pgi/1927804‑1928091
                                                                                                                                                 |                                                                                                                                              
tGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGTATACGGGTCCGGt                                                                                                                                            >  4:637748/1‑150 (MQ=255)
tGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGt                                                                                                                                            <  2:874527/150‑1 (MQ=255)
  cGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGtat                                                                                                                                          <  2:983492/150‑1 (MQ=255)
     ttAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAg                                                                                                                                       >  1:573860/1‑150 (MQ=255)
          ttGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCaaa                                                                                                                                  >  4:524167/1‑150 (MQ=255)
                cAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGttttt                                                                                                                            <  3:736451/150‑1 (MQ=255)
                             ccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTg                                                                                                               >  3:682046/1‑150 (MQ=255)
                              cgcgCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTgg                                                                                                              <  2:1210462/150‑1 (MQ=255)
                                cgcTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGaa                                                                                                            <  1:1059096/150‑1 (MQ=255)
                                cgcTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGaa                                                                                                            <  2:573860/150‑1 (MQ=255)
                                  cTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGGACg                                                                                                          >  1:490910/1‑150 (MQ=255)
                                       gCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGAt                                                                                                     <  3:90156/150‑1 (MQ=255)
                                       gCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGAt                                                                                                     >  2:510406/1‑150 (MQ=255)
                                           ttCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaaa                                                                                                 >  4:104915/1‑150 (MQ=255)
                                           ttCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGACGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTCGCTCATCCCGCATACGGGTCCGGTATAAGCCAAAGTTCTTGCCGTGTGCGTCGGAACGAAGATaaaa                                                                                                 >  4:694428/1‑150 (MQ=255)
                                           ttCAACCACCAGCCAGTCGGTAAACAGCGTAAACACTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaaa                                                                                                 >  4:105104/1‑150 (MQ=255)
                                                 caccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGt                                                                                                                       >  4:160008/1‑122 (MQ=255)
                                                 caccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGt                                                                                                                       <  3:160008/122‑1 (MQ=255)
                                                 caccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATg                                                                                           >  2:156711/1‑150 (MQ=255)
                                                   ccagccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCt                                                                                         <  1:156711/150‑1 (MQ=255)
                                                             ggTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGAt                                                                                            >  1:15085/1‑137 (MQ=255)
                                                             ggTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGAt                                                                                            <  2:15085/137‑1 (MQ=255)
                                                                 aaCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAgg                                                                           <  1:468789/150‑1 (MQ=255)
                                                                   cAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCt                                                                         <  4:682046/150‑1 (MQ=255)
                                                                    aGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTg                                                                                                                           >  4:275075/1‑99 (MQ=255)
                                                                    aGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTg                                                                                                                           <  3:275075/99‑1 (MQ=255)
                                                                      cGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGt                                                                                                                                            >  1:1160520/1‑80 (MQ=255)
                                                                      cGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGt                                                                                                                                            <  2:1160520/80‑1 (MQ=255)
                                                                      cGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGt                                                                                   >  2:305447/1‑137 (MQ=255)
                                                                      cGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGt                                                                                   <  1:305447/137‑1 (MQ=255)
                                                                        tAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCa                                                                                                        >  3:355750/1‑114 (MQ=255)
                                                                        tAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCa                                                                                                        <  4:355750/114‑1 (MQ=255)
                                                                           aCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATaa                                                                                                                                        >  1:43380/1‑79 (MQ=255)
                                                                           aCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATaa                                                                                                                                        <  2:43435/79‑1 (MQ=255)
                                                                           aCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATaa                                                                                                                                        <  2:43380/79‑1 (MQ=255)
                                                                           aCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATaa                                                                                                                                        >  1:43435/1‑79 (MQ=255)
                                                                               ttCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCt                                                                         <  2:618190/138‑1 (MQ=255)
                                                                               ttCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCt                                                                         >  1:618190/1‑138 (MQ=255)
                                                                               ttCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCt                                                                         >  2:185062/1‑138 (MQ=255)
                                                                               ttCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATc                                                             >  1:162989/1‑150 (MQ=255)
                                                                               ttCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCt                                                                         <  1:185062/138‑1 (MQ=255)
                                                                                   aGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGc                                                                                                         <  2:168588/102‑1 (MQ=255)
                                                                                   aGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGc                                                                                                         >  1:168588/1‑102 (MQ=255)
                                                                                     cATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGttttt                                                                                                                            >  1:450398/1‑81 (MQ=255)
                                                                                     cATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGttttt                                                                                                                            <  2:450398/81‑1 (MQ=255)
                                                                                     cATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcg                                                       <  3:524167/150‑1 (MQ=255)
                                                                                       tgatgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggc                                                     >  2:1021616/1‑150 (MQ=255)
                                                                                       tgatgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggc                                                     >  1:1066557/1‑150 (MQ=255)
                                                                                        gatgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGAt                                                                                                     <  2:168922/101‑1 (MQ=255)
                                                                                        gatgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGAt                                                                                                     >  1:168922/1‑101 (MQ=255)
                                                                                        gatgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcg                                                    <  1:425145/150‑1 (MQ=255)
                                                                                        gatgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcg                                                    >  2:1057441/1‑150 (MQ=255)
                                                                                         atgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCttt                                                          >  3:141640/1‑143 (MQ=255)
                                                                                         atgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCttt                                                          <  4:141640/143‑1 (MQ=255)
                                                                                         atgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcgc                                                   >  2:131887/1‑150 (MQ=255)
                                                                                           gTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcgcgg                                                 <  1:131887/150‑1 (MQ=255)
                                                                                               tCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGa                                             >  3:99675/1‑150 (MQ=255)
                                                                                                   gTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGaaaaa                                         <  1:510406/150‑1 (MQ=255)
                                                                                                    tGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaac                                        >  2:1039837/1‑150 (MQ=255)
                                                                                                                ctctTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGc                            >  4:49723/1‑150 (MQ=255)
                                                                                                                   ttCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGc                                                                                                                          <  2:473526/53‑1 (MQ=255)
                                                                                                                   ttCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGc                                                                                                                          >  1:473526/1‑53 (MQ=255)
                                                                                                                     cccACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgc                                                        >  2:566785/1‑117 (MQ=255)
                                                                                                                     cccACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgc                                                        <  1:566785/117‑1 (MQ=255)
                                                                                                                      ccACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcgcggc                                                >  3:696827/1‑124 (MQ=255)
                                                                                                                      ccACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcgcggc                                                <  4:696827/124‑1 (MQ=255)
                                                                                                                       cACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCa                                                                                                        <  2:363197/67‑1 (MQ=255)
                                                                                                                       cACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCa                                                                                                        >  1:363197/1‑67 (MQ=255)
                                                                                                                             tGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTcc               <  2:1066557/150‑1 (MQ=255)
                                                                                                                              gCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaac                                        >  4:73046/1‑124 (MQ=255)
                                                                                                                              gCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaac                                        <  3:73046/124‑1 (MQ=255)
                                                                                                                               cTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCcgc             <  3:105104/150‑1 (MQ=255)
                                                                                                                               cTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCcgc             <  3:104915/150‑1 (MQ=255)
                                                                                                                                tcatcaCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCcgcg            >  2:390862/1‑150 (MQ=255)
                                                                                                                                 catcaCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGt           >  4:246767/1‑150 (MQ=255)
                                                                                                                                 catcaCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGt           <  1:1021616/150‑1 (MQ=255)
                                                                                                                                    caCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCca        <  1:83346/150‑1 (MQ=255)
                                                                                                                                     aCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCcag       >  2:752704/1‑150 (MQ=255)
                                                                                                                                     aCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCcag       <  1:390862/150‑1 (MQ=255)
                                                                                                                                      cGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCcagc      >  4:538296/1‑150 (MQ=255)
                                                                                                                                       gCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCcagca     >  1:301117/1‑150 (MQ=255)
                                                                                                                                         aTACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGCTCAAGGCTGTATTCGTGa                                                               >  1:1068386/1‑90 (MQ=255)
                                                                                                                                         aTACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGa                                                               <  2:1068386/90‑1 (MQ=255)
                                                                                                                                          tACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAGGc  <  1:1057441/150‑1 (MQ=255)
                                                                                                                                             gggTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGcc                     <  3:529426/128‑1 (MQ=255)
                                                                                                                                             gggTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGcc                     >  4:529426/1‑128 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAGGC  >  NC_000913_3_bme_pgi/1927804‑1928091

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