Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,042,791 A→G 5.9% W120R (TGG→CGG)  yqfA ← hemolysin III family HylIII inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,042,7910AG 5.9% 196.6 / 3.3 68W120R (TGG→CGG) yqfAhemolysin III family HylIII inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base G (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAG  >  NC_000913/3042647‑3042939
                                                                                                                                                |                                                                                                                                                    
cGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCTGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCAAGGCTCCAGATaa                                                                                                                                                 >  2:545236/1‑150 (MQ=255)
cGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATaa                                                                                                                                                 <  2:887689/150‑1 (MQ=255)
  cTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCCCCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATCCCCAGCAATGCCAGGCTCCAGATAAcc                                                                                                                                               >  1:1190393/1‑149 (MQ=255)
   tGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCGCCAGATAACaa                                                                                                                                              <  2:989745/150‑1 (MQ=255)
    gCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGACTCGGTGCGCGAAGGTCAGTTTACACAGAATACCCAGCAATGCCAGGCTCCAGATAACAat                                                                                                                                             >  2:620897/1‑150 (MQ=255)
     cccGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAatc                                                                                                                                            <  1:1055832/150‑1 (MQ=255)
     cccGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAAAACACAGCAATTCCAGGCTCCAGATAACAata                                                                                                                                            >  2:626010/1‑149 (MQ=255)
       cgccgcGAGCTTAACTGCGATTTCATAATTTACCACCAGCGACAGCCAGCCCATCGGCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAatcat                                                                                                                                          <  2:1139120/150‑1 (MQ=255)
       cgccgcGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAACCAGAATACCCCGCAATGCCCGGCTCCAGATAACAatcat                                                                                                                                          >  2:941910/1‑150 (MQ=255)
         ccgcGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCCGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAAAGCCAGGCTCCAGATAACAAtcatca                                                                                                                                        >  2:106802/1‑150 (MQ=255)
         ccgcGAGCTTAACTGCCATATCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAATATTTTGAATCGGTGCGCGATGGTCAGTTTAACCAGAATACCCAGCTATGCCAGGCTCCAGATAACAAtcatca                                                                                                                                        <  1:630073/150‑1 (MQ=255)
           gcgAGCTTAACTGCTATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGTTCACCAGAGATAAAATTTTTAATCGGTGCGCGATTTTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAAc                                                                                                                                      <  1:33087/150‑1 (MQ=255)
            cgAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAAcc                                                                                                                                     >  2:1054769/1‑150 (MQ=255)
              aGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCACCCcg                                                                                                                                   >  1:1040298/1‑150 (MQ=255)
                 ttAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCcgcgc                                                                                                                                <  1:141088/150‑1 (MQ=255)
                   aaCTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACACAGCAATGCCAGGCTCCAGATAACAATCATCAACCCgcgcgc                                                                                                                              >  2:343407/1‑150 (MQ=255)
                    aCTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGGTTAAACAGCATACCCAGCAATGCCAGGCTCCAGATAACAATCATCCACCCGCGCGCt                                                                                                                             >  2:515626/1‑150 (MQ=255)
                            ttCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTGATTTTAATCAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGgaga                                                                                                                     >  1:230587/1‑150 (MQ=255)
                                tAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATcc                                                                                                                 <  1:490214/150‑1 (MQ=255)
                                   aTTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGc                                                                                                              <  1:1042810/150‑1 (MQ=255)
                                     tACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGccc                                                                                                            >  2:333811/1‑150 (MQ=255)
                                      accaccAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCGGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGcccc                                                                                                           >  1:685168/1‑150 (MQ=255)
                                      accaccAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCGGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGcccc                                                                                                           >  1:893982/1‑150 (MQ=255)
                                        caccaGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCccac                                                                                                         >  2:1116927/1‑150 (MQ=255)
                                         acccGCGCCAGCCAGCCCATCGCCAGATAGGTCACCAGCCATATAATTTTGAATCGGTGCGCGATGGTCAGGTTAAACAGAATACCCAGCACTGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAAAGGAGAATCCAGCCccacc                                                                                                        <  2:196348/150‑1 (MQ=255)
                                         accaGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCccacc                                                                                                        >  2:315373/1‑150 (MQ=255)
                                                  agccagccCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCaaaacc                                                                                               >  2:636660/1‑148 (MQ=255)
                                                           aTCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAg                                                                                      <  2:1071311/150‑1 (MQ=255)
                                                           aTCGACAGATAGGTCACCATATATAATATTTTGAATCGGTGCGCGATTGTCAGTTTAAACAGAATACCCTGAAATGACAGGCTCCATATAACTATCATCAACCCGCGGGCTAACGGAGAATCCAGCCCCACCAGCAAATACGGCGGGTAg                                                                                      <  1:634904/150‑1 (MQ=255)
                                                                   aTAGGTCAAAAGAGATAAAATTTTGAATCGGTGAGCGATGGTAAGTTTAAACAGAATACCCAGCCATTCCAGCAACCAGATAAGAATCAACATCCCGCGCGCTAACGGAGAATCAACCCACACCAGCAAAAACGGCGAGAAGGTTCCGGc                                                                              >  2:1273070/1‑150 (MQ=255)
                                                                    tAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCTATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTTACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCa                                                                             <  1:161265/150‑1 (MQ=255)
                                                                          aCCAGAGATAAAATTTTGAATCGGTGCGCGATGGGCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTACCGGAGACTCCAGCCCCACCAGCAAAAACGGCGGGTAGGTTCCGGCAATCAAc                                                                       <  2:819636/150‑1 (MQ=255)
                                                                            cAGAGATAAAATTATGAATCGGTGCGCGATGGTCAGTTTAGACAGAATACCCAGCCATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCACCAAAAACGGCGTGTAGGTTCCGGCAATCAACAg                                                                     <  2:587199/150‑1 (MQ=255)
                                                                              gagaTAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAACCAGAATACCCAGCAATGCCAGGCTCCAGCTAACAATCATCACCCCGCGCGCTAACGGAGAATCCCGCCCCACCAGCAAAAACGGCGTGGAGGTTCCGGCCATCAACAGGt                                                                   >  1:31865/1‑150 (MQ=255)
                                                                               agaTAATATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTa                                                                  <  2:936718/150‑1 (MQ=255)
                                                                                    aaaTTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAg                                                             <  1:302634/150‑1 (MQ=255)
                                                                                    aaaTTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACAGCGTGTAGGTTCCGGCAATCAACAGGTAAATAg                                                             <  1:924868/150‑1 (MQ=255)
                                                                                       ttttGAATCGGTGCGCGATGGTCAGTTAAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAgcgc                                                          <  2:88262/150‑1 (MQ=255)
                                                                                          tGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCAAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAAt                                                       <  2:6036/150‑1 (MQ=255)
                                                                                           gAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATg                                                      >  1:850402/1‑150 (MQ=255)
                                                                                              tCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTc                                                   <  2:410243/150‑1 (MQ=255)
                                                                                                 gTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTACATAGCGCAATGGTCaaa                                                >  1:1215775/1‑150 (MQ=255)
                                                                                                   gcgcGATGGTCAGTTTAAACAGAATAGCCAGCAATGCCAGGCTCCAGATAACTATCATCAACCCGCGCGCTAACGTAGAATCCAGCCCCACCAGCAAAAACGGCCTGTAGGTTCCGGCAATCAACAGGTAAATCGCGCAATGGGCAAAtt                                              <  2:697092/150‑1 (MQ=255)
                                                                                                   gcgcGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAAtt                                              <  1:541973/150‑1 (MQ=255)
                                                                                                    cgcgATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAAttt                                             <  2:1006469/150‑1 (MQ=255)
                                                                                                     gcgATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCCGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATttc                                            >  1:877872/1‑150 (MQ=255)
                                                                                                     gcgATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATttc                                            >  2:589464/1‑150 (MQ=255)
                                                                                                      cgATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATttct                                           >  2:941273/1‑150 (MQ=255)
                                                                                                      cgATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATttct                                           <  1:315373/150‑1 (MQ=255)
                                                                                                             cAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCTCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCca                                    <  1:358840/150‑1 (MQ=255)
                                                                                                              aGTTAAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACAGAGAATCCAGCCACACCAGCATAAACGTCGTGAAGGGTCCGGCTATCAACAGGTAAATCGAGCAATGGGCTACTTTCTTCAggccc                                   >  1:212131/1‑146 (MQ=255)
                                                                                                                  tAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCCCCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACAttt                               >  1:757034/1‑150 (MQ=255)
                                                                                                                      cAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGc                           <  2:739787/150‑1 (MQ=255)
                                                                                                                       aGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCTCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGcc                          <  2:230587/150‑1 (MQ=255)
                                                                                                                        gAATACCCAGCAATGCCAGGCTCCGGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGccc                         <  2:893982/150‑1 (MQ=255)
                                                                                                                        gAATACCCAGCAATGCCAGGCTCCGGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGccc                         <  2:685168/150‑1 (MQ=255)
                                                                                                                        gAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGccc                         >  2:300976/1‑150 (MQ=255)
                                                                                                                             cccAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTtga                    >  1:61011/1‑150 (MQ=255)
                                                                                                                             cccAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTtga                    >  1:450703/1‑150 (MQ=255)
                                                                                                                             cccAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCCCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCATCAGCCAAATATTTGCCAGttca                    >  1:540116/1‑148 (MQ=255)
                                                                                                                             cccAGCAATGCCAGGCTACAGATAACAATCATCAACCCGCGCGCTCACGGAGAATCAACCCCCACCAGCAAAAAAGGCCTAAAGCTTCCGCCAATCAACAGGAAAATAGGACAATGGTCAAATATCTTCAGCCACATTTTTGCCCTTtga                    >  1:803820/1‑150 (MQ=255)
                                                                                                                               cAGCAATGCCAGGCTCCAGATACCAATCGTCCACCCGCGCGCTAACGGAGCATCCAGCCCCACCAGCAAAACCGGCGTGTAGGTTCCGGCTATCAACAGGTAAATAGCGCCATGGTCAACTTTCTTCGGCGCCGTGTTTGCCCGGtgatg                  >  2:195970/1‑150 (MQ=255)
                                                                                                                               cAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTtgatg                  >  1:1132283/1‑150 (MQ=255)
                                                                                                                               cAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTtgatg                  >  1:1296580/1‑150 (MQ=255)
                                                                                                                                         aGGCTCCAGATAAAAATCATCAAACCGCGCGCTAACGGAGAATCCCGCCACACAAGTAAAAACGGCGTGTAGGTTCCGTCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCg        <  1:106802/150‑1 (MQ=255)
                                                                                                                                           gCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTg      <  2:486837/150‑1 (MQ=255)
                                                                                                                                               cAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGTTCAATTTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGAGATag  >  2:1049557/1‑150 (MQ=255)
                                                                                                                                               cAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATag  >  2:955070/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAG  >  NC_000913/3042647‑3042939

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