Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,985,515 G→A A205V (GCA→GTA)  araF ← arabinose ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,985,5150GA97.1% 208.6 / ‑4.5 68A205V (GCA→GTA) araFarabinose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (34/32);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGATATCGGCCGCTTTAAAGCCCTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATCCAGTTCGTTGGCGGTAATCGCCA  >  NC_000913/1985380‑1985647
                                                                                                                                       |                                                                                                                                    
atgatATCGGCCGCTTTAAAGCCCTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATAcccc                                                                                                                                  >  1:1403855/1‑140 (MQ=255)
atgatATCGGCCGCTNTAAAGCCCTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATAcccc                                                                                                                                  >  1:1077992/1‑140 (MQ=255)
    tatCGGCCGCTTTAAAGCCCTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCggg                                                                                                                              <  2:1313489/140‑1 (MQ=255)
       cGGCCGCTTTAAAGCCCTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATg                                                                                                                           >  1:618928/1‑140 (MQ=255)
           cGCTTTAAAGCCCTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGt                                                                                                                       <  2:1403855/140‑1 (MQ=255)
           cGCTTTAAAGCCCTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGt                                                                                                                       >  2:1234766/1‑140 (MQ=255)
                 aaaGCCCTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCGTTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGAtt                                                                                                                 <  2:537940/140‑1 (MQ=255)
                 aaaGCCCTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGAtt                                                                                                                 <  2:851293/140‑1 (MQ=255)
                     cccTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTggt                                                                                                             <  2:204634/140‑1 (MQ=255)
                       cTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTggtag                                                                                                           >  2:1043131/1‑140 (MQ=255)
                        tGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTggtagg                                                                                                          >  1:67758/1‑140 (MQ=255)
                          aCCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGt                                                                                                                       <  2:98800/125‑1 (MQ=255)
                          aCCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGt                                                                                                                       >  1:98800/1‑125 (MQ=255)
                          aCCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTggtaggta                                                                                                        >  2:1407546/1‑140 (MQ=255)
                          aCCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTggtaggta                                                                                                        >  1:1348068/1‑140 (MQ=255)
                           ccTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTAc                                                                                                       <  1:1043131/140‑1 (MQ=255)
                                    gcgcgTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATATTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAt                                                                                              >  2:1342023/1‑140 (MQ=255)
                                        gTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGt                                                                                          <  2:1348068/140‑1 (MQ=255)
                                         tACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGtt                                                                                         >  2:766455/1‑140 (MQ=255)
                                          aCGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGttt                                                                                        >  2:631332/1‑140 (MQ=255)
                                                          tGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTcc                                                                                    >  1:28217/1‑128 (MQ=255)
                                                          tGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTcc                                                                                    <  2:28217/128‑1 (MQ=255)
                                                          tGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGc                                                                         <  2:966556/139‑1 (MQ=255)
                                                          tGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGc                                                                         >  1:966556/1‑139 (MQ=255)
                                                                  ttCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                                                                <  2:952392/140‑1 (MQ=255)
                                                                  ttCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                                                                >  1:482634/1‑140 (MQ=255)
                                                                  ttCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                                                                <  2:891925/140‑1 (MQ=255)
                                                                  ttCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                                                                <  2:482634/140‑1 (MQ=255)
                                                                    cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGTAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                                                              <  1:717823/140‑1 (MQ=255)
                                                                    cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGTAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                                                              >  2:717823/1‑140 (MQ=255)
                                                                    cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                                                              <  2:483351/140‑1 (MQ=255)
                                                                    cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                                                              <  2:1434560/140‑1 (MQ=255)
                                                                    cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                                                              >  1:1335672/1‑140 (MQ=255)
                                                                    cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                                                              >  1:483351/1‑140 (MQ=255)
                                                                    cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                                                                >  1:1139324/1‑138 (MQ=255)
                                                                    cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                                                                <  1:1407888/138‑1 (MQ=255)
                                                                    cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                                                                <  2:1139324/138‑1 (MQ=255)
                                                                    cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                                                                >  2:1407888/1‑138 (MQ=255)
                                                                         cgacgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATcc                                                         >  1:637531/1‑140 (MQ=255)
                                                                           acgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                                                                >  2:614090/1‑131 (MQ=255)
                                                                           acgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                                                                <  1:614090/131‑1 (MQ=255)
                                                                            cgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATa                                                      >  1:1422934/1‑140 (MQ=255)
                                                                             gaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                                                              <  1:1471867/131‑1 (MQ=255)
                                                                             gaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                                                              >  2:1471867/1‑131 (MQ=255)
                                                                              aTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGa                                                    <  1:1407546/140‑1 (MQ=255)
                                                                              aTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGa                                                    <  2:1133693/140‑1 (MQ=255)
                                                                                cAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATc                                                  <  2:637531/140‑1 (MQ=255)
                                                                                 aGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgc                                                               >  1:904478/1‑126 (MQ=255)
                                                                                 aGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATa                                                      <  1:1228228/135‑1 (MQ=255)
                                                                                 aGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATa                                                      >  2:1228228/1‑135 (MQ=255)
                                                                                 aGCCAATGTTTAACTTACGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgc                                                               <  2:904478/126‑1 (MQ=255)
                                                                                   ccAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                                                              >  2:418500/1‑125 (MQ=255)
                                                                                   ccAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                                                              <  1:418500/125‑1 (MQ=255)
                                                                                      aTGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCgtag                                            >  1:974179/1‑140 (MQ=255)
                                                                                               ttCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTAcggcggc                                   <  2:974179/140‑1 (MQ=255)
                                                                                                 ccGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTAcggcggcgg                                 >  1:560885/1‑140 (MQ=255)
                                                                                                 ccGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTAcggcggcgg                                 >  1:119778/1‑140 (MQ=255)
                                                                                                 ccGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTAcggcggc                                   >  1:1186230/1‑138 (MQ=255)
                                                                                                 ccGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTAcggcggc                                   <  2:1186230/138‑1 (MQ=255)
                                                                                                  cGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCggcggg                                <  2:277125/140‑1 (MQ=255)
                                                                                                         tGAACCAGCATTGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATc                         <  1:1342023/140‑1 (MQ=255)
                                                                                                                    tGAGTTGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATCCAGTTCGTTgg              <  2:1077992/140‑1 (MQ=255)
                                                                                                                        ttGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATCCAGTTCGTTGGCGGt          >  1:444753/1‑140 (MQ=255)
                                                                                                                         tGGCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATCCAGTTCGTTGGCGGTa         <  1:1082554/140‑1 (MQ=255)
                                                                                                                          ggCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATCCAGTTCGTTGGCGGTaa        >  1:697197/1‑140 (MQ=255)
                                                                                                                           gCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATCCAGTTCGTTGGCGGTa         >  2:170115/1‑138 (MQ=255)
                                                                                                                           gCAGCGTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATCCAGTTCGTTGGCGGTa         <  1:170115/138‑1 (MQ=255)
                                                                                                                                gTCAAATACCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATCCAGTTCGTTGGCGGTAATCGCCa  >  1:499992/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                    
ATGATATCGGCCGCTTTAAAGCCCTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATCCAGTTCGTTGGCGGTAATCGCCA  >  NC_000913/1985380‑1985647

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