Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,119,486 C→A 7.6% P726T (CCG→ACG)  PP_3622 → isoquinoline 1‑oxidoreductase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,119,4860CA 7.6% 243.0 / 8.6 79P726T (CCG→ACG) PP_3622isoquinoline 1‑oxidoreductase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (36/37);  new base A (4/2);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

GGTGTACGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGCTGGCGCTGCCGTGCC  >  NC_002947/4119338‑4119633
                                                                                                                                                    |                                                                                                                                                   
ggTGTACGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGccc                                                                                                                                                     >  7:78164/1‑149 (MQ=255)
   gtACGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCcgg                                                                                                                                                  >  1:376567/1‑149 (MQ=255)
      cGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCcggcg                                                                                                                                                <  8:208974/148‑1 (MQ=255)
        tGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTg                                                                                                                                             <  3:267522/149‑1 (MQ=255)
        tGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTg                                                                                                                                             <  2:10391/149‑1 (MQ=255)
        tGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTg                                                                                                                                             >  7:194965/1‑149 (MQ=255)
             ggCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGTGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCaa                                                                                                                                        >  6:365166/1‑149 (MQ=255)
              gCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCaa                                                                                                                                        <  2:376567/148‑1 (MQ=255)
                aaGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAAcgc                                                                                                                                     >  7:62815/1‑149 (MQ=255)
                  gCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCgg                                                                                                                                   <  5:34954/149‑1 (MQ=255)
                       gcgcACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTcgcg                                                                                                                              <  8:349786/149‑1 (MQ=255)
                         gcacgcaATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTcgcg                                                                                                                              >  2:405636/1‑147 (MQ=255)
                           acgcaATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAg                                                                                                                           >  6:438341/1‑148 (MQ=255)
                           acgcaATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGc                                                                                                                          >  7:109767/1‑149 (MQ=255)
                           acgcaATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGc                                                                                                                          >  6:429127/1‑149 (MQ=255)
                            cgcaATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCt                                                                                                                         >  4:262376/1‑149 (MQ=255)
                             gcaATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTc                                                                                                                        >  8:354489/1‑149 (MQ=255)
                               aaTTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCa                                                                                                                       >  4:419236/1‑148 (MQ=255)
                                        cTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAgcg                                                                                                             >  6:126176/1‑149 (MQ=255)
                                             cccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTgc                                                                                                        <  6:362023/149‑1 (MQ=255)
                                             cccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTgc                                                                                                        <  1:405636/149‑1 (MQ=255)
                                             cccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTgc                                                                                                        >  8:112674/1‑149 (MQ=255)
                                              ccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTgcg                                                                                                       >  5:385543/1‑149 (MQ=255)
                                                 aTCCTACCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAg                                                                                                    <  3:247232/149‑1 (MQ=255)
                                                          ccTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCAtt                                                                                           <  3:297378/149‑1 (MQ=255)
                                                            ttCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGa                                                                                         <  8:78164/149‑1 (MQ=255)
                                                               gcgGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGAT‑GGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAgcc                                                                                      <  8:327215/148‑1 (MQ=255)
                                                                 ggATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTTGCCAGCGCGTGCGCAGTTTGCCATTGAgccgc                                                                                    >  6:422361/1‑149 (MQ=255)
                                                                 ggATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAgccgc                                                                                    <  2:267805/149‑1 (MQ=255)
                                                                      gcgcGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCcaca                                                                                <  5:219564/148‑1 (MQ=255)
                                                                         cgGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCtt                                                                            <  8:62815/149‑1 (MQ=255)
                                                                          gggTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCtt                                                                            >  6:21795/1‑148 (MQ=255)
                                                                           ggTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCa                                                                          <  6:385543/149‑1 (MQ=255)
                                                                            gTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAg                                                                         >  8:52714/1‑149 (MQ=255)
                                                                             tGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGc                                                                        <  5:424989/149‑1 (MQ=255)
                                                                                aTGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGa                                                                     >  2:167321/1‑149 (MQ=255)
                                                                                    gcgcATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGc                                                                 >  1:402155/1‑149 (MQ=255)
                                                                                        atcatcGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGcac                                                               <  6:227257/147‑1 (MQ=255)
                                                                                        atcatcGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGcac                                                               >  5:227257/1‑147 (MQ=255)
                                                                                           atcGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCTCTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAgg                                                           <  1:118511/148‑1 (MQ=255)
                                                                                                  gCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTc                                                   <  8:34375/149‑1 (MQ=255)
                                                                                                   cGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGa                                                                     >  2:159661/1‑130 (MQ=255)
                                                                                                   cGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGa                                                                     <  1:159661/130‑1 (MQ=255)
                                                                                                     gggCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCCCCACACAGGCCCGCGCCTGCGCAGTTTGCCATTGAGCCCCCACACCTTAAACGGATAGCACAAAGGACGCCTCAt                                                 >  2:426763/1‑148 (MQ=255)
                                                                                                       gCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGAc                                              <  5:365166/149‑1 (MQ=255)
                                                                                                       gCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGAc                                              <  6:170485/149‑1 (MQ=255)
                                                                                                         gAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGAccc                                            <  2:234515/149‑1 (MQ=255)
                                                                                                                ggcggcATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGttt                                     <  7:354489/149‑1 (MQ=255)
                                                                                                                 gcggcATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTc                                    <  5:126176/149‑1 (MQ=255)
                                                                                                                  cggcATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCg                                   >  2:264013/1‑149 (MQ=255)
                                                                                                                   ggcATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGc                                  <  5:422361/149‑1 (MQ=255)
                                                                                                                     cATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACTCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCaa                                <  5:438341/149‑1 (MQ=255)
                                                                                                                     cATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCAGTGGGCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGc                                  >  1:405259/1‑147 (MQ=255)
                                                                                                                     cATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCaa                                <  3:419236/149‑1 (MQ=255)
                                                                                                                      aTCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTgc                                                                                                        >  5:437176/1‑76 (MQ=255)
                                                                                                                      aTCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTgc                                                                                                        <  6:437176/76‑1 (MQ=255)
                                                                                                                      aTCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCCCCCCTTCAGCTGCTCGCACCAGGGCCGCCTCATGACCCATCGGTTCGc                                  >  7:159600/1‑146 (MQ=255)
                                                                                                                        cGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGaa                             >  1:294642/1‑149 (MQ=255)
                                                                                                                         ggCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAAc                            >  6:143715/1‑149 (MQ=255)
                                                                                                                         ggCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCCCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAAc                            >  2:24728/1‑149 (MQ=255)
                                                                                                                         ggCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCACCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTACGCTGATCGCACAAGGGACGCCTCATGACCCAACGTTTCGCAAGAAc                            >  6:42002/1‑149 (MQ=255)
                                                                                                                          gCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCCCACCTTCAGCTGATCGCACACGGGACGCCTCATGACCCATCGTTTCGCCAGAAcc                           >  5:18348/1‑149 (MQ=255)
                                                                                                                          gCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGttt                                     >  4:425824/1‑139 (MQ=255)
                                                                                                                          gCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGttt                                     <  3:425824/139‑1 (MQ=255)
                                                                                                                                 accGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgctggct                    <  1:426763/149‑1 (MQ=255)
                                                                                                                                 accGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgctggct                    >  4:412358/1‑149 (MQ=255)
                                                                                                                                 accGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgctggct                    <  1:264013/149‑1 (MQ=255)
                                                                                                                                    gCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgctggctggc                 >  3:238138/1‑149 (MQ=255)
                                                                                                                                     gccaccgccgGGGCCACGGCGGTGGCCAAAGAGGTCGCGCAGCTCACTGGCAAGCGAGTGCGCAGTTTGCCATTGAGACGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgctggctggct                >  2:310518/5‑149 (MQ=255)
                                                                                                                                      ccaccgccgGGGCCACGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgctgg                      <  2:291751/139‑1 (MQ=255)
                                                                                                                                      ccaccgccgGGGCCACGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgctgg                      >  1:291751/4‑142 (MQ=255)
                                                                                                                                        cccGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCtggctggctggc             <  3:262376/149‑1 (MQ=255)
                                                                                                                                        cccGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGACTCATGACCCATCGTTTAGCAAGAACCgctggctgtctg               >  8:197678/1‑147 (MQ=255)
                                                                                                                                         ccgccgTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGCTGgcg            >  7:205940/1‑149 (MQ=255)
                                                                                                                                           gccgTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGCTGGCGCt          <  8:205940/149‑1 (MQ=255)
                                                                                                                                            ccgTGGCCCCGGCGGTGGCCAACGCGGTCTCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGCTGGCGCTg         <  7:197678/149‑1 (MQ=255)
                                                                                                                                             gttGGCCACCGCCGGGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgctggctggc                 <  7:77554/138‑1 (MQ=255)
                                                                                                                                             gttGGCCACCGCCGGGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgctggctggc                 >  8:77554/3‑140 (MQ=255)
                                                                                                                                             gttGGCCACCGCCGGGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGCTGGCGCTGc        >  8:380782/3‑149 (MQ=255)
                                                                                                                                             cgTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGCTGGCGCTGc        <  5:429127/149‑1 (MQ=255)
                                                                                                                                                   ccccGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGACATTGAGACGCCCCACCTTCAGCTGATCGCACAAGGGACTCCTCATGACCCATAGTTTCGCAAGCACCGCTGGCTGGCTGGCGCTGCCGTGcc  >  1:324728/1‑149 (MQ=255)
                                                                                                                                                   cccGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCGAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGCTGGCGCTGCCGTGcc  <  2:101838/149‑1 (MQ=255)
                                                                                                                                                   cccGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGCTGGCGCTGCCGTGcc  >  8:118666/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGTGTACGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGCTGGCGCTGCCGTGCC  >  NC_002947/4119338‑4119633

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