Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 901,646 A→C F202C (TTT→TGT)  artQ ← arginine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913901,6460AC98.4% 229.2 / NA 61F202C (TTT→TGT) artQarginine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (31/29);  total (32/29)
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

AGCCCTTTCATCAGTTCGGGTAAATACTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTAAATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAGCCACTGGTTACCGAGGCCAGGCAGCGCATGACGCCACATCTGCGGCAT  >  NC_000913/901500‑901791
                                                                                                                                                  |                                                                                                                                                 
aGCCCTTTCATCAGTGCGGGTAAATCCTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGGGCAGGTCAATGCGTTTGAGAGTGTACAGACTGAGCGGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATg                                                                                                                                                <  2:205371/150‑1 (MQ=255)
aGCCCTTCCATCAGTTCGGGTAAATTCTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCTTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATg                                                                                                                                                <  2:427852/150‑1 (MQ=255)
  cccTTTCATCAGTTCGGGTAAATACTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCCCCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGtt                                                                                                                                             <  2:277855/151‑1 (MQ=255)
   ccTTTCATCAGTTCGGGTAAATACTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTc                                                                                                                                            <  1:295301/151‑1 (MQ=255)
     tttCATCAGTTCGGGTAAATACTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGGGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGAACATGGTTcc                                                                                                                                           >  2:427372/1‑150 (MQ=255)
             gTTCGGGTAAATCCTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGa                                                                                                                                  <  2:857495/151‑1 (MQ=255)
             gTTCGGGTAAATACTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGa                                                                                                                                  <  1:679721/151‑1 (MQ=255)
                cGGGTAAATACTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTa                                                                                                                               <  1:861379/151‑1 (MQ=255)
                cGGGTAAATACTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTa                                                                                                                               <  2:282694/151‑1 (MQ=255)
                 gggTAAATACTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGt                                                                                                                                >  2:484409/1‑149 (MQ=255)
                    tAAATACTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCg                                                                                                                            >  2:275903/1‑150 (MQ=255)
                     aaaTACTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGAt                                                                                                                          >  1:605286/1‑151 (MQ=255)
                         aCTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCtt                                                                                                                      >  1:245873/1‑151 (MQ=255)
                              aaCATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTtttgt                                                                                                                 <  2:521714/151‑1 (MQ=255)
                                cATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTtttgttt                                                                                                               >  2:1843/1‑151 (MQ=255)
                                 aTTAGCTGGGACTCCGCTCAAAACGTGGAGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGATCAGGGTGATCACCAGGTCAATCGCCGCCGCCACAATGAACCAGGTACATGGTTACTGGGTACGAGTAGCGATGCttt                                                                                                                     >  2:722329/1‑144 (MQ=255)
                                    aGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTgc                                                                                                             >  2:646755/1‑149 (MQ=255)
                                       tGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCAt                                                                                                        >  1:629163/1‑151 (MQ=255)
                                       tGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCAt                                                                                                        >  1:81656/1‑151 (MQ=255)
                                         ggCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCAtt                                                                                                       >  1:822108/1‑150 (MQ=255)
                                         ggCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTa                                                                                                      <  2:1977/151‑1 (MQ=255)
                                          gCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTaa                                                                                                     <  2:498594/151‑1 (MQ=255)
                                          gCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTaa                                                                                                     <  2:217735/151‑1 (MQ=255)
                                           ccTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTaaa                                                                                                    <  1:292672/151‑1 (MQ=255)
                                            cTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAAt                                                                                                   >  2:82296/1‑151 (MQ=255)
                                             tCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATc                                                                                                  >  1:186683/1‑151 (MQ=255)
                                              ccGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATc                                                                                                  <  2:219240/150‑1 (MQ=255)
                                                gCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCAtt                                                                                               >  1:80/1‑151 (MQ=255)
                                                  tCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTc                                                                                              >  2:757921/1‑150 (MQ=255)
                                                           gtCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATc                                                                                     <  2:653255/150‑1 (MQ=255)
                                                              gcgcgcAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCaaa                                                                                  <  2:118204/150‑1 (MQ=255)
                                                              gcgcgcAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAAc                                                                                 <  1:873580/151‑1 (MQ=255)
                                                              gcgcgcAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAAc                                                                                 <  2:49840/151‑1 (MQ=255)
                                                              gcgcgcAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAAc                                                                                 >  2:375327/1‑151 (MQ=255)
                                                                     ggTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGAcc                                                                            >  1:202854/1‑149 (MQ=255)
                                                                      gTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAgc                                                                         >  2:277804/1‑151 (MQ=255)
                                                                        cAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCGGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAgcg                                                                        >  2:360744/1‑150 (MQ=255)
                                                                        cAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAgcgc                                                                       >  2:260830/1‑151 (MQ=255)
                                                                          aTGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCgg                                                                     >  2:896095/1‑151 (MQ=255)
                                                                           tGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGt                                                                    >  1:79362/1‑151 (MQ=255)
                                                                            gCGTTTGAGAATGTACTGACTGAGCAGGGTGATCCCCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGt                                                                    <  1:823368/150‑1 (MQ=255)
                                                                             cGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTAt                                                                  >  1:531570/1‑151 (MQ=255)
                                                                               tttGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCt                                                                >  2:770858/1‑151 (MQ=255)
                                                                               tttGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCt                                                                >  1:281849/1‑151 (MQ=255)
                                                                                ttGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCtt                                                               >  1:495191/1‑151 (MQ=255)
                                                                                ttGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCt                                                                <  1:468227/150‑1 (MQ=255)
                                                                                    gaATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTca                                                            <  2:682279/150‑1 (MQ=255)
                                                                                             gactgaGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCa                                                   <  2:511562/150‑1 (MQ=255)
                                                                                             gactgaGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAg                                                  <  2:898075/151‑1 (MQ=255)
                                                                                                   gCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAGCCACTg                                            >  1:12878/1‑151 (MQ=255)
                                                                                                      gggTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAGCCACTGGtt                                         >  1:754658/1‑151 (MQ=255)
                                                                                                      gggTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAGCCACTGGt                                          <  2:452200/150‑1 (MQ=255)
                                                                                                                     aaaTCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAGCCACTGGTTACCGAGGCCAGGCa                           <  1:310622/150‑1 (MQ=255)
                                                                                                                       aTCGCCGCCGCCACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAGCCACTGGTTACCGAGGCCAGGCAgcg                        >  2:626824/1‑151 (MQ=255)
                                                                                                                           ccgccgccACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAGCCACTGGTTACCGAGGCCAGGCAGCGCAt                     <  2:112372/150‑1 (MQ=255)
                                                                                                                              ccgccACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAGCCACTGGTTACCGAGGCCAGGCAGCGCATGACg                 <  1:98670/151‑1 (MQ=255)
                                                                                                                               cgccACAATGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAGCCACTGGTTACCGAGGCCAGGCAGCGCATGACGc                <  1:82296/151‑1 (MQ=255)
                                                                                                                                     aaTGTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAGCCACTGGTTACCGAGGCCAGGCAGCGCATGACGCCACATc          <  2:313447/151‑1 (MQ=255)
                                                                                                                                        gTACCAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAGCCACTGGTTACCGAGGCCAGGCAGCGCATGACGCCACATCTg        >  2:262259/1‑150 (MQ=255)
                                                                                                                                           ccAGGTACATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAGCCACTGGTTACCGAGGCCAGGCAGCGCATGACGCCACATCTGCGGc    >  1:321263/1‑151 (MQ=255)
                                                                                                                                             aGGTAAATGGTTCCTGGGTACGCGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAGCCACTGGTTACCGAGGCCAGGCAGCGCATGACGCCACATCTGCGGCAt  >  2:507830/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGCCCTTTCATCAGTTCGGGTAAATACTCAAACATTAGCTGGGCCTCCGCTCAAAACGTGTCGCGCGCAGGTCAATGCGTTTGAGAATGTACTGACTGAGCAGGGTGATCACCAGGTAAATCGCCGCCGCCACAATGTACCAGGTAAATGGTTCCTGGGTACGAGTAGCGATGCTTTTTGTTTGCAGCATTAAATCATTCACACTAATCAAACTGACCAGCGCGGTATCTTTCAGCAGCACCAGCCACTGGTTACCGAGGCCAGGCAGCGCATGACGCCACATCTGCGGCAT  >  NC_000913/901500‑901791

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