Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,038,117 C→G 5.8% H126Q (CAC→CAG cbdA → cytochrome bd‑II oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,038,1170CG 5.8% 370.8 / 11.3 103H126Q (CAC→CAGcbdAcytochrome bd‑II oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base C (50/47);  new base G (3/3);  total (53/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.01e-01

Only 100 of 106 total aligned reads displayed.

TACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTC  >  NC_000913/1037973‑1038264
                                                                                                                                                |                                                                                                                                                   
tACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAcc                                                                                                                                                    >  3:213105/1‑146 (MQ=255)
tACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAcc                                                                                                                                                    <  4:213105/146‑1 (MQ=255)
 aCTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCCCCTTTGTCGGGCTGTTCTTATTCGGCTGGCCACGGCTGAATAAATACCAGCACctgc                                                                                                                                                 >  4:719409/1‑148 (MQ=255)
  cTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctgct                                                                                                                                                <  3:719409/148‑1 (MQ=255)
       tGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGa                                                                                                                                           >  1:681303/1‑148 (MQ=255)
          tCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGt                                                                                                                                        >  1:939222/1‑148 (MQ=255)
                          aTGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCgg                                                                                                                        <  2:939222/148‑1 (MQ=255)
                                  gATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTTGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtc                                                                                                                <  2:681303/148‑1 (MQ=255)
                                   atatTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGt                                                                                                                            >  2:26846/1‑135 (MQ=255)
                                   atatTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGt                                                                                                                            <  1:26846/135‑1 (MQ=255)
                                               cACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTCCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGta                                                                                                                                            <  2:784233/107‑2 (MQ=38)
                                               cACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTCCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGta                                                                                                                                            >  1:784233/1‑106 (MQ=38)
                                                        cGATGGAAGCATTAATGGCCTTCTTCCTCGAAGCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTtgt                                                                                                    >  3:970858/1‑138 (MQ=255)
                                                        cGACGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTtgt                                                                                                    <  4:970858/138‑1 (MQ=255)
                                                          aTGGAAGCATTAATGGCCTTCTTCCTCGAATCCACTTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGt                                                                                                                                             >  4:1079288/1‑95 (MQ=255)
                                                          aTGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGt                                                                                                                                             <  3:1079288/95‑1 (MQ=255)
                                                          aTGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGa                                                                                           <  1:403243/145‑1 (MQ=255)
                                                          aTGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGa                                                                                           >  2:403243/1‑145 (MQ=255)
                                                            ggAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGcc                                                                                      >  3:281388/1‑148 (MQ=255)
                                                             gAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGaa                                                                                          <  3:913430/143‑1 (MQ=255)
                                                             gAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGaa                                                                                          >  4:913430/1‑143 (MQ=255)
                                                              aaGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCaa                                                                                    >  4:544650/1‑148 (MQ=255)
                                                               aGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAAc                                                                                   >  2:899742/1‑148 (MQ=255)
                                                                  aTTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGt                                                                                <  4:281388/148‑1 (MQ=255)
                                                                   ttAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGc                                                                                                                                     >  3:612927/1‑94 (MQ=255)
                                                                   ttAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGc                                                                                                                                     <  4:612927/94‑1 (MQ=255)
                                                                       tGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGAt                                                                           <  3:544650/148‑1 (MQ=255)
                                                                        ggCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctctc                                                                                                            >  4:1138334/1‑114 (MQ=255)
                                                                         gCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGt                                                                                                                                        <  1:289650/85‑1 (MQ=255)
                                                                         gCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGt                                                                                                                                        >  2:289650/1‑85 (MQ=255)
                                                                            ttcttcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGta                                                                                                                                            >  2:280171/1‑77 (MQ=255)
                                                                            ttcttcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGta                                                                                                                                            <  1:280171/78‑2 (MQ=255)
                                                                            ttcttcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCtggt                                                                                                                                 >  3:970031/1‑89 (MQ=255)
                                                                            ttcttcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCtggt                                                                                                                                 <  4:970031/89‑1 (MQ=255)
                                                                            ttcttcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctct                                                                                                             <  2:896414/109‑1 (MQ=255)
                                                                            ttcttcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctct                                                                                                             >  1:896414/1‑109 (MQ=255)
                                                                                tcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctct                                                                                                             >  2:157865/1‑105 (MQ=255)
                                                                                tcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctct                                                                                                             <  1:157865/105‑1 (MQ=255)
                                                                                 ccTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCg                                                                 <  1:899742/148‑1 (MQ=255)
                                                                                   tCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctgc                                                                                                                                                 <  1:581175/66‑1 (MQ=255)
                                                                                   tCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctgc                                                                                                                                                 >  2:581175/1‑66 (MQ=255)
                                                                                   tCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGatat                                                                                              >  3:927569/1‑117 (MQ=255)
                                                                                   tCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGatat                                                                                              <  4:927569/117‑1 (MQ=255)
                                                                                    cGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGGACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtc                                                                                                                <  4:662102/98‑1 (MQ=255)
                                                                                    cGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtc                                                                                                                >  3:662102/1‑98 (MQ=255)
                                                                                    cGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGa                                                                                           <  4:848016/119‑1 (MQ=255)
                                                                                    cGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGa                                                                                           >  3:848016/1‑119 (MQ=255)
                                                                                     gAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGa                                                                                           <  3:498883/118‑1 (MQ=255)
                                                                                       aTCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGt                                                           >  1:142256/1‑148 (MQ=255)
                                                                                        tCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTg                                                          <  2:142256/148‑1 (MQ=255)
                                                                                         tcACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGtt                                                                                                                      >  4:791790/2‑87 (MQ=255)
                                                                                         ccACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctc                                                                                                              >  4:597657/1‑95 (MQ=255)
                                                                                         ccACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctc                                                                                                              <  3:597657/95‑1 (MQ=255)
                                                                                           aCCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCaa                                                                                                                   >  4:170497/1‑88 (MQ=255)
                                                                                           aCCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCaa                                                                                                                   <  3:170497/88‑1 (MQ=255)
                                                                                            ccTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctctc                                                                                                            >  4:428925/1‑94 (MQ=255)
                                                                                            ccTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctctc                                                                                                            <  3:428925/94‑1 (MQ=255)
                                                                                             cTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGt                                                                                >  1:1048524/1‑121 (MQ=255)
                                                                                             cTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGt                                                                                <  2:1048524/121‑1 (MQ=255)
                                                                                               ttGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctg                                                                                                                                                  >  2:587623/1‑53 (MQ=255)
                                                                                               ttGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctg                                                                                                                                                  <  1:587623/53‑1 (MQ=255)
                                                                                                   cGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACg                                                                                                                                         >  1:437368/1‑58 (MQ=255)
                                                                                                      gCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctctc                                                                                                            <  3:554565/84‑1 (MQ=255)
                                                                                                      gCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctctc                                                                                                            >  4:554565/1‑84 (MQ=255)
                                                                                                      gCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGa                                                                                           >  2:35700/1‑101 (MQ=255)
                                                                                                      gCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGa                                                                                           <  1:35700/101‑1 (MQ=255)
                                                                                                          ttcttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGaca                                        >  2:772771/1‑148 (MQ=255)
                                                                                                           tcttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctct                                                                                                             >  3:270959/1‑78 (MQ=255)
                                                                                                           tcttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctct                                                                                                             <  4:270959/78‑1 (MQ=255)
                                                                                                           tcttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTtgt                                                                                                    >  1:1068161/1‑87 (MQ=255)
                                                                                                           tcttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTtgt                                                                                                    <  2:1068161/87‑1 (MQ=255)
                                                                                                              tcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGa                                                >  3:987195/1‑136 (MQ=255)
                                                                                                                ttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCa                                                      >  2:757715/1‑128 (MQ=255)
                                                                                                                 tcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGatat                                                                                              >  1:889863/1‑87 (MQ=255)
                                                                                                                 tcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGatat                                                                                              <  2:889863/87‑1 (MQ=255)
                                                                                                                    gCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTa                               <  3:783104/147‑1 (MQ=255)
                                                                                                                    gCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTa                               >  4:783104/1‑147 (MQ=255)
                                                                                                                           aCGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTg                                   >  3:838357/1‑136 (MQ=255)
                                                                                                                           aCGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTg                                   <  4:838357/136‑1 (MQ=255)
                                                                                                                                    acgtcaCCAGCAGGTGCTGGTAACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTg                                                 >  1:70807/6‑113 (MQ=255)
                                                                                                                                    acgtcaCCAGCAGGTGCTGGTAACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTg                                                 <  2:70807/108‑1 (MQ=255)
                                                                                                                                     aaataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCa             >  2:364346/1‑148 (MQ=255)
                                                                                                                                      aataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATTCAATACCCGACCGGTGCGCATTTTg                                                 <  4:723417/111‑1 (MQ=255)
                                                                                                                                      aataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTg                                                 >  3:723417/1‑111 (MQ=255)
                                                                                                                                        taCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAGTCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATg                       <  2:634519/135‑1 (MQ=255)
                                                                                                                                        taCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGt                                                                                <  4:306663/78‑1 (MQ=255)
                                                                                                                                        taCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGt                                                                                >  3:306663/1‑78 (MQ=255)
                                                                                                                                        taCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTg                                   <  4:396727/123‑1 (MQ=255)
                                                                                                                                        taCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTg                                   >  3:396727/1‑123 (MQ=255)
                                                                                                                                        taCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATg                       >  1:634519/1‑135 (MQ=255)
                                                                                                                                         aCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGatat                                             >  3:127623/1‑112 (MQ=255)
                                                                                                                                         aCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGatat                                             <  4:127623/112‑1 (MQ=255)
                                                                                                                                         aCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGa         >  4:544708/1‑148 (MQ=255)
                                                                                                                                           cAGCACCTTCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTg                                   <  1:61177/120‑1 (MQ=255)
                                                                                                                                           aaGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTg                                   >  2:61177/2‑120 (MQ=255)
                                                                                                                                             gCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGa                      <  1:904036/131‑1 (MQ=255)
                                                                                                                                             gCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGa                      >  2:904036/1‑131 (MQ=255)
                                                                                                                                              cACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATc                                            >  1:922168/1‑108 (MQ=255)
                                                                                                                                              cACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATc                                            <  2:922168/108‑1 (MQ=255)
                                                                                                                                                ccTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTc  <  1:364346/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTC  >  NC_000913/1037973‑1038264

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