Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 49,424:1 +C 5.8% coding (1226/1251 nt) czcB‑I ← cobalt‑zinc‑cadmium resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294749,4241.C 5.8% 351.3 / 11.5 86coding (1226/1251 nt)czcB‑Icobalt‑zinc‑cadmium resistance protein
Reads supporting (aligned to +/- strand):  ref base . (36/45);  new base C (2/3);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

ATAGCTGTAGATACCGATACCCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTACCGGGCGGGTGTCAAAGCCATCTGGGGTGCGTAC  >  NC_002947/49282‑49566
                                                                                                                                               |                                                                                                                                              
aTAGCTGTAGATTCCGATACCCTCCATGATCAGAACGGCGATCATTACTACGATGCGATGCTCGATGGAGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGcc                                                                                                                                            >  5:270093/1‑147 (MQ=255)
aTAGCTGTAGATACCGATACCCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCg                                                                                                                                           <  2:231608/148‑1 (MQ=255)
 tAGCTGTAGATACCGATACCCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGa                                                                                                                                          <  1:176154/148‑1 (MQ=255)
 tAGCTGTAGATACCGATACCCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAg                                                                                                                                         <  3:285145/149‑1 (MQ=255)
      gTAGATACCGATACCCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCg                                                                                                                                     >  8:75770/1‑148 (MQ=255)
       tAGATACCGATACCCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCg                                                                                                                                     <  2:273841/147‑1 (MQ=255)
       tAGATACCGATACCCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCg                                                                                                                                     <  5:267056/147‑1 (MQ=255)
       tAGATACCGATACCCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGa                                                                                                                                   <  7:75770/149‑1 (MQ=255)
        aGATACCGATACCAGCCATTATCAGAACGGCGATCATTACTACTATGCGCTGCTAGATGGAGAATCTGATGATACGTTCAAACATTAGAAATCATCCTGTCAGTTCGCATCTGTGACCGTGTTCGGCAGAAGCTTTTCCGAGCTCGGa                                                                                                                                   <  6:270093/148‑1 (MQ=255)
          ataccgataccCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACtt                                                                                                                                >  1:403129/1‑149 (MQ=255)
              cgataccCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACtt                                                                                                                                >  1:167978/1‑145 (MQ=255)
              cgataccCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACtt                                                                                                                                <  2:167978/145‑1 (MQ=255)
              cgataccCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAgt                                                                                                                            <  2:115349/149‑1 (MQ=255)
              cgataccCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAgt                                                                                                                            <  6:165126/149‑1 (MQ=255)
              cgataccCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAAAGCTTCTGCCGaa                                                                                                                                         >  2:187126/1‑135 (MQ=255)
              cgataccCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAAAGCTTCTGCCGaa                                                                                                                                         <  1:187126/136‑2 (MQ=255)
                    ccGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAg                                                                                                                      >  5:305264/1‑149 (MQ=255)
                       ccATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAgt                                                                                                                            >  7:374431/1‑140 (MQ=255)
                       ccATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAgt                                                                                                                            <  8:374431/140‑1 (MQ=255)
                       ccATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTg                                                                                                                   >  5:35674/1‑149 (MQ=255)
                            aTCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAAAGCTTCTGCCGAACACggt                                                                                                                                   <  6:47938/128‑2 (MQ=255)
                            aTCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAAAGCTTCTGCCGAACACggt                                                                                                                                   >  5:47938/1‑127 (MQ=255)
                             tCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTg                                                                                                             >  1:146932/1‑149 (MQ=255)
                             tCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATCATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTg                                                                                                             >  4:353685/1‑149 (MQ=255)
                               aGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGt                                                                                                           <  2:403129/149‑1 (MQ=255)
                                gAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGtt                                                                                                          >  1:270855/1‑149 (MQ=255)
                                gAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGtt                                                                                                          >  5:75075/1‑149 (MQ=255)
                                 aaCGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTg                                                                                                         >  8:380114/1‑149 (MQ=255)
                                   cGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCt                                                                                                       <  6:35674/149‑1 (MQ=255)
                                    ggCGATCATTACTTCGATGCGCTGCTCGATGGCGAATCTGTTGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGTAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTa                                                                                                      <  1:152503/149‑1 (MQ=255)
                                    ggCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTa                                                                                                      >  5:351717/1‑149 (MQ=255)
                                              actacGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTCTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTg                                                                                                             <  1:61128/133‑1 (MQ=255)
                                              actacGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTg                                                                                                             >  2:61128/1‑132 (MQ=255)
                                              actacGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGcc                                                                                            >  1:38615/1‑149 (MQ=255)
                                                tacGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGa                                                                                                                              >  7:22099/1‑113 (MQ=255)
                                                tacGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGa                                                                                                                              <  8:22099/113‑1 (MQ=255)
                                                   gATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCtt                                                                                       <  7:380114/149‑1 (MQ=255)
                                                    aTGCGCTGCTAGATTTCTAATCTGATGATACGTTCAAACATGTGAAATCATCCTGGCAGTTCGCATCAGTGACCGTGTTCGTCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGTTTGCTAGCTGGGCGCCGGCTTc                                                                                      <  2:181977/149‑1 (MQ=255)
                                                        gctgctCGATGGCGAATCTGATGATACGTTCAAACATGAGGAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATa                                                                                  >  6:26773/1‑149 (MQ=255)
                                                           gctCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCg                                                                               >  7:121919/1‑149 (MQ=255)
                                                           gctCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCg                                                                               >  6:80579/1‑149 (MQ=255)
                                                            ctCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGt                                                                              <  8:330741/149‑1 (MQ=255)
                                                                 tGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCgatg                                                                         <  2:146932/149‑1 (MQ=255)
                                                                   gCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCgatga                                                                        <  3:353685/148‑1 (MQ=255)
                                                                     gAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGAttt                                                                     >  5:86299/1‑149 (MQ=255)
                                                                           gatgatACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACg                                                                <  6:351717/148‑1 (MQ=255)
                                                                           gatgatACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACg                                                                >  8:2092/1‑148 (MQ=255)
                                                                            atgatACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCa                                                                  <  8:316576/145‑1 (MQ=255)
                                                                             tgatACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTAc                                                             >  8:338100/1‑149 (MQ=255)
                                                                             tgatACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTAc                                                             <  6:305264/149‑1 (MQ=255)
                                                                               atACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTATTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTAcc                                                            >  2:201970/1‑148 (MQ=255)
                                                                                tACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCa                                                          <  5:80579/149‑1 (MQ=255)
                                                                                 aCGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCAt                                                         >  2:249405/1‑149 (MQ=255)
                                                                                  cGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCAtt                                                        >  2:339297/1‑149 (MQ=255)
                                                                                   gTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACtt                                                                                                                                <  8:190090/76‑1 (MQ=255)
                                                                                   gTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACtt                                                                                                                                >  7:190090/1‑76 (MQ=255)
                                                                                   gTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTg                                                       >  7:259605/1‑149 (MQ=255)
                                                                                     tCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACtt                                                                                                                                >  7:232967/1‑74 (MQ=255)
                                                                                     tCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACtt                                                                                                                                <  8:232967/74‑1 (MQ=255)
                                                                                       aaaCATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATAc                                                                                 <  6:421842/119‑1 (MQ=255)
                                                                                       aaaCATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATAc                                                                                 >  5:421842/1‑119 (MQ=255)
                                                                                        aaCATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTcgcg                                                  <  2:270855/149‑1 (MQ=255)
                                                                                         aCATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCgcgg                                                 <  7:338100/149‑1 (MQ=255)
                                                                                           aTGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCgcggc                                                >  2:20069/1‑148 (MQ=255)
                                                                                             gagaAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGAc                                             >  3:264527/1‑149 (MQ=255)
                                                                                              agaAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGAc                                             <  6:201676/148‑1 (MQ=255)
                                                                                                aaaTCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGAcc                                            >  1:97409/1‑147 (MQ=255)
                                                                                                   tcatcCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGtt                                       <  5:26773/149‑1 (MQ=255)
                                                                                                     atcCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGtttt                                     <  1:20069/149‑1 (MQ=255)
                                                                                                        cTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTAc                                   <  6:86299/148‑1 (MQ=255)
                                                                                                           tCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTACcgg                                <  8:259605/148‑1 (MQ=255)
                                                                                                           tCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCCTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTGTTACcggg                               >  7:366372/1‑149 (MQ=255)
                                                                                                             aGTTCGCAAGCGTGACCGGGTTCCGCCGCAGCTTTGCCGAGCTCGGACTTGAGCGCGAAGCTGCCACTGGTAGCTACCTGGGCGCCGGCTTCTATACCGGCGATGAGTTGCAGGAACCCATTGTCGCGGCGACCCAGTTTTACcgggc                              <  6:303365/148‑1 (MQ=255)
                                                                                                                 cGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTACCGGGCGGgt                          <  2:97409/148‑1 (MQ=255)
                                                                                                                   cATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTACCGGGCGGgtgt                        >  1:365629/1‑148 (MQ=255)
                                                                                                                   cATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTACCGGGCGGGTGTc                       <  4:264527/149‑1 (MQ=255)
                                                                                                                     tCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTACCGGGCGGGTGTCa                      >  1:372596/1‑148 (MQ=255)
                                                                                                                       aGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTACCGGGCGGGTGTCaaa                    <  2:38615/148‑1 (MQ=255)
                                                                                                                       aGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTACCGGGCGGGTGTCaaa                    <  6:75075/148‑1 (MQ=255)
                                                                                                                       aGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTACCGGGCGGGTGTCAAAg                   <  2:365629/149‑1 (MQ=255)
                                                                                                                            ccGTGTTCGGCAGTAGATTTGCCGAGCTCGGACTTGAGTGTGCAGCTGCCACTGGTTGCTACCTGGGCGCCGGCCTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTACCGGGCGGGTGTCAAAGCCATc              <  8:366372/149‑1 (MQ=255)
                                                                                                                                gtTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCg                                                                            >  8:101576/1‑83 (MQ=255)
                                                                                                                                gtTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCg                                                                            <  7:101576/83‑1 (MQ=255)
                                                                                                                                 ttCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTACCGGGCGGGTGTCAAAGCCATCTggg          <  1:249405/148‑1 (MQ=255)
                                                                                                                                 ttCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTACCGGGCGGGTGTCAAAGCCATCTggg          <  1:339297/148‑1 (MQ=255)
                                                                                                                                        gAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTACCGGGCGGGTGTCAAAGCCATCTGGGGTGCGTAc  <  1:201970/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
ATAGCTGTAGATACCGATACCCGCCATGATCAGAACGGCGATCATTACTACGATGCGCTGCTCGATGGCGAATCTGATGATACGTTCAAACATGAGAAATCATCCTGTCAGTTCGCATCAGTGACCGTGTTCGGCAGAAGCTTTGCCGAGCTCGGACTTGAGTGTGAAGCTGCCACTGGTTGCTACCTGGGCGCCGGCTTCAATACCGTCGATGATTTCCACGTACCCATTGTCGCGGCGACCCAGTTTTACCGGGCGGGTGTCAAAGCCATCTGGGGTGCGTAC  >  NC_002947/49282‑49566

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: