Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,039,479 G→A 17.1% V62I (GTC→ATC) ‡ cbdB → cytochrome bd‑II oxidase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,039,4790GA17.1% 194.1 / 29.3 87V62I (GTC→ATC) cbdBcytochrome bd‑II oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base G (42/30);  new base A (6/9);  total (48/39)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

GACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTG  >  NC_000913/1039380‑1039577
                                                                                                   |                                                                                                  
gACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAgg                                                                                                    <  2:1654942/100‑1 (MQ=255)
gACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAgg                                                                                                    <  2:1654427/100‑1 (MQ=255)
  cATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTc                                                                                                  <  1:1907253/100‑1 (MQ=255)
    tGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTg                                                                                                >  2:1925797/1‑100 (MQ=255)
     ggggATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTgg                                                                                               <  1:1259541/100‑1 (MQ=255)
       ggATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGtt                                                                                             <  2:367441/100‑1 (MQ=255)
       ggATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGtt                                                                                             >  1:1510104/1‑100 (MQ=255)
       ggATCGGCCGTGTGCTGCCGCTGGGGGCGCGTAATGATGATGAGCGCCGTACAGCGCTAAACAGCGTGGGTGCATACTGGGAAGGCAACCAGGTCTGGt                                                                                              >  2:989312/1‑99 (MQ=255)
          tCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGTCAACCAGACCTGGTTgcc                                                                                          <  1:424427/100‑3 (MQ=255)
          tCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGAt                                                                                          <  2:1868723/100‑1 (MQ=255)
           cGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATc                                                                                         >  2:732108/1‑100 (MQ=255)
           cGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTgcct                                                                                         <  1:1584325/100‑4 (MQ=255)
            ggCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTgcctt                                                                                        <  1:2214810/100‑5 (MQ=255)
             gCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCt                                                                                       <  2:1012349/100‑1 (MQ=255)
             gCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCt                                                                                       >  1:92055/1‑100 (MQ=255)
             gCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTgccttc                                                                                       <  1:294435/100‑6 (MQ=255)
             gCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTgccttc                                                                                       <  1:947225/100‑6 (MQ=255)
             gCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGGGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATca                                                                                        >  1:1319248/1‑98 (MQ=255)
              cTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTgccttcc                                                                                      <  1:825109/100‑7 (MQ=255)
               tGTCTGCTGCCGGTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTgccttccc                                                                                     <  2:503767/100‑8 (MQ=255)
               tGTCTGCTGCCGCTGGTTGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGCTGCCCACTGGGAAGGCAACCAGACCTGGTTgccttccc                                                                                     <  2:477166/100‑8 (MQ=255)
               tGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTgccttccc                                                                                     <  1:1491173/100‑8 (MQ=255)
                 tCTGCTGCCACTGGTGGCGCGTAATTATGATGATCGCCGGATAGTGATAAACCGCTTTGGTGCACACTGGGATGGCAACCAGGTCTGGTTGATCGTCGCt                                                                                   <  2:993163/100‑1 (MQ=255)
                   tgctgcCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTATGGTTGATCCTCGCtgg                                                                                 >  2:1939911/1‑100 (MQ=255)
                    gctgcCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCtggt                                                                                >  1:2063766/1‑100 (MQ=255)
                     ctgcCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCtggtg                                                                               >  2:118644/1‑100 (MQ=255)
                       gcCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTggg                                                                             >  1:1241562/1‑100 (MQ=255)
                        ccGCGGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTgggg                                                                            <  2:530653/100‑1 (MQ=255)
                          gCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCa                                                                          >  2:955771/1‑100 (MQ=255)
                          gCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAAGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCa                                                                          >  2:449440/1‑100 (MQ=255)
                             ggtggCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCatta                                                                       <  2:1820978/100‑1 (MQ=255)
                             ggtggCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCatta                                                                       >  1:600491/1‑100 (MQ=255)
                                 gcgcGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTg                                                                   <  2:2378932/100‑1 (MQ=255)
                                 gcgcGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTg                                                                   <  2:2280032/100‑1 (MQ=255)
                                 gcgcGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTg                                                                   <  1:2409577/100‑1 (MQ=255)
                                  cgcgTAGTGATGATGAGCGCCGGATAGTGATAGACTGCGTTGGTGCACACTGGGGAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTgc                                                                  <  2:2514754/100‑1 (MQ=255)
                                   gcgTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTgcc                                                                 >  2:672829/1‑100 (MQ=255)
                                     gTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTgccgc                                                               >  1:1576412/1‑100 (MQ=255)
                                      tAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTgccgcc                                                              >  1:2343621/1‑100 (MQ=255)
                                      tAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTgccgcc                                                              <  2:1241554/100‑1 (MQ=255)
                                        atgatgatgAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTg                                                            >  2:1217908/1‑100 (MQ=255)
                                         tgatgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTgg                                                           >  1:389972/1‑100 (MQ=255)
                                         tgatgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTgg                                                           >  1:750769/1‑100 (MQ=255)
                                         tgatgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTgg                                                           >  1:1986072/1‑100 (MQ=255)
                                           atgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGcc                                                         >  1:599628/1‑100 (MQ=255)
                                           atgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGcc                                                         <  2:900758/100‑1 (MQ=255)
                                           atgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGc                                                          >  2:6136/1‑99 (MQ=255)
                                            tgatgaACGCCGGATAGTGATAAATAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGccc                                                        >  2:2298517/1‑100 (MQ=255)
                                            tgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGccc                                                        >  1:1088256/1‑100 (MQ=255)
                                            tgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGccc                                                        >  1:1136557/1‑100 (MQ=255)
                                            tgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGccc                                                        >  1:181665/1‑100 (MQ=255)
                                             gatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCa                                                       >  1:1136217/1‑100 (MQ=255)
                                               tgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCaga                                                     >  2:1931997/1‑100 (MQ=255)
                                               tgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCaga                                                     >  1:2198314/1‑100 (MQ=255)
                                                gaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCagag                                                    >  1:1231335/1‑100 (MQ=255)
                                                  aCGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAgtg                                                  <  2:1106473/100‑1 (MQ=255)
                                                  aCGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAgtg                                                  <  1:169981/100‑1 (MQ=255)
                                                   cGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAgtg                                                  <  1:1176537/99‑1 (MQ=255)
                                                     ccGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCTCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTAt                                               <  2:1576401/100‑1 (MQ=255)
                                                      cGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATg                                              <  1:1843928/100‑1 (MQ=255)
                                                      cGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATg                                              <  1:401138/100‑1 (MQ=255)
                                                       ggATAGTGGTAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGc                                             <  1:732115/100‑1 (MQ=255)
                                                        tatatagATTAACTGCGTGGGTTCACACTGTGAAGGCAACTAGGTATGGATGATCCTCGCTGGTGGGGCATTATTAGCCGCCAGGCTCAGAGTGTAtcc                                             <  2:472416/93‑3 (MQ=255)
                                                            gTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAgcgg                                        >  2:76105/1‑100 (MQ=255)
                                                            gTGACAAACGGCGTTGGTGGACACTGGGGAGGCATCCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTTTATGCAgcgg                                        <  2:537845/100‑1 (MQ=255)
                                                             tGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAgcggc                                       >  2:51085/1‑100 (MQ=255)
                                                               aTAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGt                                     >  2:102847/1‑100 (MQ=255)
                                                                 aaaCAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGttt                                   >  2:2496621/1‑100 (MQ=255)
                                                                   aCAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGCGTGTATGCAGCGGCGTTTTc                                 >  2:664179/1‑100 (MQ=255)
                                                                    cAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTc                                 >  2:834381/1‑99 (MQ=255)
                                                                     aGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCg                               <  1:1276621/100‑1 (MQ=255)
                                                                        gTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCt                            <  1:1932011/100‑1 (MQ=255)
                                                                            gTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTTCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCttn                          <  2:2324172/98‑2 (MQ=255)
                                                                               cacaCTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTAtgt                     >  2:1940607/1‑100 (MQ=255)
                                                                                 cacTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTgg                   <  2:2384981/100‑1 (MQ=255)
                                                                                  acTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGc                  <  2:207593/100‑1 (MQ=255)
                                                                                  acTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGc                  <  1:57808/100‑1 (MQ=255)
                                                                                   cTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATCATGTGCTGCTTGTGCCAGGGTGTATACAGCGGTGTTTTCCGGCTTTTAtgt                     >  2:251527/1‑96 (MQ=255)
                                                                                   agcgaggatcAACCAGACCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCg                 >  1:1009888/10‑100 (MQ=255)
                                                                                     gggAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCgat               >  2:671209/1‑100 (MQ=255)
                                                                                     gggAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGAGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCgat               >  1:1806216/1‑100 (MQ=255)
                                                                                      gaggatcAACCAGACCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCgatg              >  1:465147/7‑100 (MQ=255)
                                                                                       aggatcAACCAGACCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCga                >  2:2142941/6‑97 (MQ=255)
                                                                                        ggatcAACCAGACCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCgatgat            >  2:1300023/5‑100 (MQ=255)
                                                                                        ggatcAACCAGACCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCgatgat            >  2:2040631/5‑100 (MQ=255)
                                                                                        ggatcAACCAGACCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCgatg              >  2:1631108/5‑98 (MQ=255)
                                                                                          ggCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATcc          >  1:2201318/1‑100 (MQ=255)
                                                                                              aCCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgt      <  1:1642624/100‑1 (MQ=255)
                                                                                               ccAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgtg     >  1:30487/1‑100 (MQ=255)
                                                                                                  ggTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgtgctg  >  1:583882/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTG  >  NC_000913/1039380‑1039577

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