Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,669,731 G→C 9.3% T21T (ACG→ACC nuoN → NADH‑quinone oxidoreductase subunit N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,669,7310GC 9.3% 188.1 / 7.9 64T21T (ACG→ACCnuoNNADH‑quinone oxidoreductase subunit N
Reads supporting (aligned to +/- strand):  ref base G (33/25);  new base C (2/4);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

ACCTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATC  >  NC_002947/4669587‑4669875
                                                                                                                                                 |                                                                                                                                                
acCTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtgg                                                                                                                                              >  5:157719/1‑149 (MQ=255)
  ctctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtg                                                                                                                                            <  2:235870/149‑1 (MQ=255)
  ctctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtg                                                                                                                                            <  6:45517/149‑1 (MQ=255)
   tctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCAGATGCTGATCACCACCATCACCACggtggtgg                                                                                                                                           >  2:271263/1‑149 (MQ=255)
     tGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACTCTTCATCGCATTGGCGCCGATGCTGATCACCAGCATCACCACggtggtgg                                                                                                                                           <  5:342362/147‑1 (MQ=255)
      gccgccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGa                                                                                                                                        >  8:104129/1‑149 (MQ=255)
       ccgccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtgg                                                                                                                                           <  4:263335/145‑1 (MQ=255)
       ccgccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtgg                                                                                                                                           >  3:263335/1‑145 (MQ=255)
         gccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGc                                                                                                                                     <  4:360329/149‑1 (MQ=255)
         gccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGc                                                                                                                                     >  5:176355/1‑149 (MQ=255)
            accaTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTgg                                                                                                                                  >  1:189432/1‑149 (MQ=255)
            accaTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTgg                                                                                                                                  <  5:117785/149‑1 (MQ=255)
            accaTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTgg                                                                                                                                  <  8:152789/149‑1 (MQ=255)
             ccaTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGc                                                                                                                                 >  2:361744/1‑149 (MQ=255)
              caTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCAGCATCACCACCACCGTGGTGAtggtggtg                                                                                                                                <  3:213057/149‑3 (MQ=255)
               aTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGa                                                                                                                               <  1:361744/149‑1 (MQ=255)
                tGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGAt                                                                                                                              >  1:195811/1‑149 (MQ=255)
                  aGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATc                                                                                                                             <  5:302971/148‑1 (MQ=255)
                      gtgtGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGTTGATGCTGGCGATCGCCTg                                                                                                                        >  2:52143/1‑149 (MQ=255)
                      gtgtGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCt                                                                                                                         <  5:72670/148‑1 (MQ=255)
                      gtgtGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCACCACCACCGTGGTGATGgtggtgatca                                                                                                                               >  1:287254/1‑138 (MQ=255)
                      gtgtGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCACCACCACCGTGGTGATGgtggtgatca                                                                                                                               <  2:287254/142‑5 (MQ=255)
                       tgtgCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtg                                                                                                                                               >  5:115987/1‑125 (MQ=255)
                       tgtgCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtg                                                                                                                                               <  6:115987/125‑1 (MQ=255)
                         tgCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGaa                                                                                                                     <  2:205426/149‑1 (MQ=255)
                         tgCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGaa                                                                                                                     <  2:195811/149‑1 (MQ=255)
                            agcagTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAgcg                                                                                                                  >  5:306999/1‑149 (MQ=255)
                             gcagTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAgcgc                                                                                                                 <  6:176355/149‑1 (MQ=255)
                                 tGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAAcc                                                                                                             <  2:189432/149‑1 (MQ=255)
                                   gCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCa                                                                                                            <  5:372126/148‑1 (MQ=255)
                                   gCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCAc                                                                                                           <  7:226711/149‑1 (MQ=255)
                                     tCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTg                                                                                                                                   <  2:315009/123‑1 (MQ=255)
                                     tCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTg                                                                                                                                   >  1:315009/1‑123 (MQ=255)
                                                 tCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACCACCGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTc                                                                                               <  1:214385/147‑1 (MQ=255)
                                                    cTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCctgct                                                                                          <  3:73702/149‑1 (MQ=255)
                                                    cTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCACCACCgtggtg                                                                                                                                               >  2:115760/1‑96 (MQ=255)
                                                     tCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCg                                                                                                        >  2:182281/1‑134 (MQ=255)
                                                     tCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCg                                                                                                        <  1:182281/134‑1 (MQ=255)
                                                     tCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCctgctg                                                                                         >  5:284778/1‑149 (MQ=255)
                                                     tCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCctgctg                                                                                         >  6:350006/1‑149 (MQ=255)
                                                     tCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCctgctg                                                                                         >  4:163457/1‑149 (MQ=255)
                                                              ttCGGCCCGGTAAGAGCGCTATGGAATTCATCCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGACGATGGCGATCGCCTGGAAgc                                                                                                                   >  1:226002/1‑115 (MQ=255)
                                                              ttCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGACGCTGGCGATCGCCTGGAAgc                                                                                                                   <  2:226002/114‑1 (MQ=255)
                                                                cGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGAGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTggg                                                                              >  7:80769/1‑149 (MQ=255)
                                                                  gCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGcc                                                                            >  1:110729/1‑149 (MQ=255)
                                                                        tAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAAcc                                                                      <  7:104129/149‑1 (MQ=255)
                                                                                     gAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGCCATGAACCTGCCCctcattac                                                         >  8:236290/1‑143 (MQ=255)
                                                                                         tCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTc                                                     >  4:105923/1‑149 (MQ=255)
                                                                                              cacTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGc                                                 >  2:296536/1‑149 (MQ=255)
                                                                                              cacTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGc                                                 >  5:348400/1‑149 (MQ=255)
                                                                                               acTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGcc                                                >  5:115472/1‑149 (MQ=255)
                                                                                                   aaCACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGgcg                                            >  1:326766/1‑149 (MQ=255)
                                                                                                     cacTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCCCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCt                                          >  6:7587/1‑149 (MQ=255)
                                                                                                        ttCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTg                                         >  5:290618/1‑147 (MQ=255)
                                                                                                        ttCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCACAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGaa                                       <  5:7587/149‑1 (MQ=255)
                                                                                                           aTCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAgg                                     <  5:350006/148‑1 (MQ=255)
                                                                                                           aTCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGt                                    >  7:78475/1‑149 (MQ=255)
                                                                                                            tCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctg                                                              <  1:131269/122‑1 (MQ=255)
                                                                                                            tCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctg                                                              >  2:131269/1‑122 (MQ=255)
                                                                                                                    gcgcCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCg                                                        <  7:312032/120‑1 (MQ=255)
                                                                                                                    gcgcCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCg                                                        >  8:312032/1‑120 (MQ=255)
                                                                                                                      gcCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTgg                         >  7:114196/1‑149 (MQ=255)
                                                                                                                           tGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTc                    <  2:326766/149‑1 (MQ=255)
                                                                                                                              tGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTcaca                 >  4:65260/1‑149 (MQ=255)
                                                                                                                              tGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTcaca                 >  8:110643/1‑149 (MQ=255)
                                                                                                                              tGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTcaca                 >  8:234576/1‑149 (MQ=255)
                                                                                                                               gATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACAt                >  7:263335/1‑149 (MQ=255)
                                                                                                                                   gccagcaTCACCACCACCGTGGTGATGGTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCgcgc                               <  2:185484/129‑1 (MQ=38)
                                                                                                                                   gccagcaTCACCACCACCGTGGTGATGGTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCgcgc                               >  1:185484/2‑130 (MQ=38)
                                                                                                                                   accaccATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCg              <  8:78475/147‑1 (MQ=255)
                                                                                                                                     caccaTCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGa          >  2:334613/1‑149 (MQ=255)
                                                                                                                                       ccaTCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGAt         >  4:204128/1‑148 (MQ=255)
                                                                                                                                            accacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCa    <  2:37784/148‑1 (MQ=255)
                                                                                                                                             ccacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATc  <  3:204128/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACCTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATC  >  NC_002947/4669587‑4669875

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