Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,663,917 G→A V103V (GTG→GTA ynfH → oxidoreductase, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,663,9170GA98.3% 208.2 / ‑4.4 58V103V (GTG→GTAynfHoxidoreductase, membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (30/27);  total (31/27)
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.57e-01

TCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCT  >  NC_000913/1663775‑1664064
                                                                                                                                              |                                                                                                                                                   
tCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCAt                                                                                                                                             >  1:295387/1‑151 (MQ=255)
   ggTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCgg                                                                                                                                          <  1:646458/151‑1 (MQ=255)
      ttATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTa                                                                                                                                        >  2:598665/1‑150 (MQ=255)
       tATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTaa                                                                                                                                       >  2:725894/1‑150 (MQ=255)
          tGCCTCTGTCATGCTCCTGGGATCGCCACTGCGGGCGTTTAACTCACGGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTCCCATATTTTTTGCCGTTGCCGGATTATGGTGGATGGTAGGAGTCATCGGTAAAAt                                                                                                                                    <  1:253271/150‑1 (MQ=255)
           gCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGc                                                                                                                                  >  2:832385/1‑151 (MQ=255)
                   cATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGc                                                                                                                           <  1:415681/150‑1 (MQ=255)
                   cATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCg                                                                                                                          >  1:252619/1‑151 (MQ=255)
                     tGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGtt                                                                                                                        >  2:268060/1‑151 (MQ=255)
                      gCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGtt                                                                                                                        <  1:562013/150‑1 (MQ=255)
                      gCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTg                                                                                                                       <  1:523372/151‑1 (MQ=255)
                      gCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTg                                                                                                                       <  1:393562/151‑1 (MQ=255)
                       cACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTgg                                                                                                                      <  1:6278/151‑1 (MQ=255)
                        aCCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTggg                                                                                                                     >  1:722580/1‑151 (MQ=255)
                         ccTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTgggg                                                                                                                    <  2:76320/151‑1 (MQ=255)
                                 cGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTg                                                                                                            >  1:10942/1‑151 (MQ=255)
                                   ccACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGt                                                                                                          <  1:387549/151‑1 (MQ=255)
                                       tGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACt                                                                                                      >  1:607948/1‑151 (MQ=255)
                                        gCGGGCGTTTAACTCCCTGAATCGCATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGGGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACtt                                                                                                     >  1:502639/1‑151 (MQ=255)
                                        gCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACtt                                                                                                     >  1:655732/1‑151 (MQ=255)
                                         cGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACttt                                                                                                    <  2:11879/151‑1 (MQ=255)
                                             cGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTc                                                                                                  >  2:530212/1‑149 (MQ=255)
                                                 tAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCAt                                                                                             <  1:803510/150‑1 (MQ=255)
                                                  aaCTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATgg                                                                                           >  2:50474/1‑151 (MQ=255)
                                                       ccTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTc                                                                                      <  2:647233/151‑1 (MQ=255)
                                                          gaatcgaatcGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGc                                                                                   >  1:171391/1‑151 (MQ=255)
                                                           aatcgaatcgGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCg                                                                                  <  2:568650/151‑1 (MQ=255)
                                                              cgaatcgGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTc                                                                                <  2:171414/150‑1 (MQ=255)
                                                               gaatcgGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCAt                                                                              >  2:799767/1‑151 (MQ=255)
                                                                 atcgGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCAttt                                                                            <  2:234524/151‑1 (MQ=255)
                                                                   cgGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTc                                                                          <  2:735576/151‑1 (MQ=255)
                                                                   cgGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTc                                                                          <  1:741428/151‑1 (MQ=255)
                                                                     gAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGt                                                                        <  2:695415/151‑1 (MQ=255)
                                                                     gAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGt                                                                        <  2:216959/151‑1 (MQ=255)
                                                                          tCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTgg                                                                    >  1:32480/1‑150 (MQ=255)
                                                                           cgggTCTGAGTAACGAGCTTGCTGCCGGTTCCATATTCTTTGCCGTTGGCGGATGATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTgg                                                                    <  1:819830/147‑1 (MQ=255)
                                                                             ggTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGgat                                                                  <  2:99683/149‑1 (MQ=255)
                                                                              gTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGgatgat                                                               <  2:179898/151‑1 (MQ=255)
                                                                               tCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGatgatg                                                              >  1:774908/1‑151 (MQ=255)
                                                                                cTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGatgatg                                                              >  2:120502/1‑150 (MQ=255)
                                                                                  gAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGatgatg                                                              <  2:254608/148‑1 (MQ=255)
                                                                                   aGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCt                                                          >  2:791976/1‑151 (MQ=255)
                                                                                              tgctgcCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCa                                                >  1:149094/1‑150 (MQ=255)
                                                                                               gctgcCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCaaa                                              <  2:238579/151‑1 (MQ=255)
                                                                                                ctgcCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAAt                                             >  1:783583/1‑151 (MQ=255)
                                                                                                       ttCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACAcc                                      >  2:756620/1‑151 (MQ=255)
                                                                                                         ccATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGt                                    >  1:74395/1‑151 (MQ=255)
                                                                                                         ccATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGt                                    <  1:685304/151‑1 (MQ=255)
                                                                                                           atatTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGc                                  >  1:42125/1‑151 (MQ=255)
                                                                                                            tatTTTTTGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGcc                                 >  2:514869/1‑151 (MQ=255)
                                                                                                                  ttGCCGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCt                            <  1:181631/150‑1 (MQ=255)
                                                                                                                     ccGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGc                         >  2:289587/1‑150 (MQ=255)
                                                                                                                      cGTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCAt                       >  1:787743/1‑151 (MQ=255)
                                                                                                                       gTTGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCAt                       >  2:68949/1‑150 (MQ=255)
                                                                                                                         tGGCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAAc                    >  2:813340/1‑151 (MQ=255)
                                                                                                                           gCGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACgg                  >  1:477091/1‑151 (MQ=255)
                                                                                                                            cGGATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGt                 >  1:286403/1‑151 (MQ=255)
                                                                                                                             ggATTATGGTGGCTGGTAGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGtt                <  1:484895/151‑1 (MQ=255)
                                                                                                                                        gctggtagcAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTg     <  1:394151/151‑1 (MQ=255)
                                                                                                                                          tggtagcAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGc   >  1:307513/1‑151 (MQ=255)
                                                                                                                                           ggtagcAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCt  <  2:221343/151‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCT  >  NC_000913/1663775‑1664064

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