Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,263,969 A→T 7.8% L93Q (CTG→CAG)  zinU ← P‑loop guanosine triphosphatase‑dependent zinc chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,263,9690AT 7.8% 186.7 / 5.1 64L93Q (CTG→CAG) zinUP‑loop guanosine triphosphatase‑dependent zinc chaperone
Reads supporting (aligned to +/- strand):  ref base A (27/32);  new base T (3/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01

GTCTGGGTCAGGTGCAGTTCGGCATGCGCGGCATCGACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGACCTGCACCGGCTGTTCCCCCAGCAGTTGGCTGTCGAT  >  NC_002947/5263825‑5264117
                                                                                                                                                |                                                                                                                                                    
gTCTGGGTCAGGTGCAGTTCGGCATGCGCGGCATCGACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGc                                                                                                                                                   >  3:339490/1‑148 (MQ=255)
       tCAGGTGCAGTTCGGCATGCGCGGCATCGACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCGCCTGGTGATCGaat                                                                                                                                           <  6:171775/149‑2 (MQ=255)
       tCAGGTGCAGTTCGGCATGCGCGGCATCGACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAg                                                                                                                                           <  5:347505/149‑1 (MQ=255)
        cAGGTGCAGTTCGGCATGCGCGGCATCGACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCGCCTGGTGATCGaatg                                                                                                                                          <  1:163019/149‑3 (MQ=255)
           gTGCAGTTCGGCATGCGCGGCATCGACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGa                                                                                                                                       <  8:1666/149‑1 (MQ=255)
           gTGCAGTTCGGCATGCGCGGCATCGACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGa                                                                                                                                       >  2:25717/1‑149 (MQ=255)
           gTGCAGTTCCGCATGAGCGTCATCAACAAGGCTGATGATGACGACGACGAACTACAGCTCACGCAACTCTTCTTCCATGAAGAAACTTTGCGCGACCGGCGCCGGGTCGGAAAGCCCGCTACATACGATCACCAGGCGATCGAAGGcta                                                                                                                                       >  8:246260/1‑147 (MQ=255)
             gCAGTTCGGCATGCGCGGCATCGACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGAtt                                                                                                                                     >  4:118890/1‑149 (MQ=255)
               aGTTCGGCATGCGCGGCATCGACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTc                                                                                                                                   <  1:242636/149‑1 (MQ=255)
                        tGCGCGGCCTCGACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATc                                                                                                                          >  5:239901/1‑149 (MQ=255)
                               cATCGACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGt                                                                                                                   <  1:25717/149‑1 (MQ=255)
                                   gACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTc                                                                                                                 >  5:185291/1‑147 (MQ=255)
                                   gACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCa                                                                                                                <  4:339490/148‑1 (MQ=255)
                                   gACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATAACCAGGCGATCGAAGGCGATTTCCCCGGCATCCATGCGTTc                                                                                                                 >  2:89272/1‑147 (MQ=255)
                                    aCAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAag                                                                                                              >  6:374644/1‑149 (MQ=255)
                                           tgatgaTGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAgg                                                                                                       <  6:185291/149‑1 (MQ=255)
                                           tgatgaTGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAgg                                                                                                       <  7:246260/149‑1 (MQ=255)
                                            gatgatGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGt                                                                                                      <  2:227928/149‑1 (MQ=255)
                                              tgatgCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTa                                                                                                     <  3:245370/148‑1 (MQ=255)
                                               gatgCCGTCGAGGATGTGGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTa                                                                                                     >  5:42456/1‑147 (MQ=255)
                                                atgCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGGCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACa                                                                                                   <  2:250104/148‑1 (MQ=255)
                                                    cGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCGCCAGGCGCTCGAGAGCGCTTTGCCCGGCAGCCCGGCGATCAAGCAGAAGGTGCAgctcg                                                                                              >  8:151880/1‑149 (MQ=255)
                                                     gTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAgcgcgc                                                                                             <  3:127594/149‑1 (MQ=255)
                                                         aGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGt                                                                                         <  3:274143/149‑1 (MQ=255)
                                                              gTAGCGCTCACGCAGCTCTTCGTCGATGAAGCACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGt                                                                                    >  5:365008/1‑149 (MQ=255)
                                                               tAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTc                                                                                   >  1:231843/1‑149 (MQ=255)
                                                                  cgcTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGAt                                                                                                                                                <  6:228110/85‑1 (MQ=255)
                                                                  cgcTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGAt                                                                                                                                                >  5:228110/1‑85 (MQ=255)
                                                                   gcTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCa                                                                               >  4:183775/1‑149 (MQ=255)
                                                                   gcTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCa                                                                               <  5:374644/149‑1 (MQ=255)
                                                                     tCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATg                                                                             >  1:73832/1‑149 (MQ=255)
                                                                       aCGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCAAGTCAGGTCGCCATGGa                                                                           <  1:15594/149‑1 (MQ=255)
                                                                        cGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCATCAGGTACAGCGCGCGAGTCAGGTCGCCATGGAt                                                                          <  1:43741/149‑1 (MQ=255)
                                                                         gCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGt                                                                                                                   >  4:113756/1‑107 (MQ=255)
                                                                         gCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGt                                                                                                                   <  3:113756/107‑1 (MQ=255)
                                                                         gCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGAtg                                                                         <  7:55962/149‑1 (MQ=255)
                                                                               ctTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTgcag                                                                   >  6:338296/1‑149 (MQ=255)
                                                                                   gtcgATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATa                                                               >  5:280315/1‑149 (MQ=255)
                                                                                   gtcgATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATa                                                               <  6:42456/149‑1 (MQ=255)
                                                                                      gATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATaca                                                             >  5:71728/1‑148 (MQ=255)
                                                                                           ataaCGTTTGCGCGACAGGAGCCTGGTCGGACAGGCAGTGGAATTCGCTAACCAGTCGATCGAAGGCGATTTCCCAGGCATCAAGGCGTGCAAGCAGCAGGTACAGCGCGCGAGTCAGGGCGCCATGGATGCTGCAGCATACACAAcc                                                        <  7:156066/146‑1 (MQ=255)
                                                                                             aaCGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGtt                                                     <  5:27472/149‑1 (MQ=255)
                                                                                                    gcgcGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCaa                                               <  6:71728/148‑1 (MQ=255)
                                                                                                    gcgcGACCGGCGCCGGGTCGGACAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAg                                              >  2:67178/1‑149 (MQ=255)
                                                                                                     cgcgACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAgg                                             <  6:365008/149‑1 (MQ=255)
                                                                                                        gACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCa                                          >  4:133566/1‑149 (MQ=255)
                                                                                                          ccGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATg                                        <  3:118890/149‑1 (MQ=255)
                                                                                                            ggCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAATGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGAc                                      <  8:211242/149‑1 (MQ=255)
                                                                                                            ggCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGAc                                      <  5:338296/149‑1 (MQ=255)
                                                                                                             gcgcCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGAcc                                     >  7:39175/1‑149 (MQ=255)
                                                                                                              cgcCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGACCt                                    <  6:239901/149‑1 (MQ=255)
                                                                                                                     tCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGACCTGCACCgg                             <  3:133566/149‑1 (MQ=255)
                                                                                                                         ccAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGACCTGCACCGGCTGt                         >  4:65861/1‑149 (MQ=255)
                                                                                                                            ggCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGACCTGCACCGGCTGTTcc                      <  2:231843/149‑1 (MQ=255)
                                                                                                                               cGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGAAGGTACAGAGCGCGAGTCAGGTCGACATGGATGCTGCAGCATACACAACGGTTGGCCAAGGTCATGACCTGCAGCGGCTGTTcccac                   >  4:224519/1‑147 (MQ=255)
                                                                                                                                ggTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGACCTGCACCGGCTGTTCCCCca                  <  6:280315/149‑1 (MQ=255)
                                                                                                                                    cATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGACCTGCACCGGCTGTTCCCCcagca               <  2:73832/148‑1 (MQ=255)
                                                                                                                                     cttCGATCGCCTGGTGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGTGCGAGTCAGGTCGCCATGGATGCTGCATCATACAAAACCGTTGGCCAAGGTCAGGACCTGCACCGGCTGTTCCCCcagccgg             >  8:351631/2‑146 (MQ=255)
                                                                                                                                     cttCGATCGCCTGGTGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTTGCCAAGGTCATGACCTGCACCGGCTGTTCCCCCAGCAGt             >  7:240535/2‑149 (MQ=255)
                                                                                                                                     cttCGATCGCCTGGTGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGACCTGCACCGGCTGTTCCCCCAGCAGt             >  6:248087/2‑149 (MQ=255)
                                                                                                                                      ttCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGACCTGCACCGGCTGTTCCCCCAGCAGtt            >  2:18343/1‑149 (MQ=255)
                                                                                                                                         gATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGACCTGCACCGGCTGTTCCCCCAGCAGTTgg          <  1:18343/148‑1 (MQ=255)
                                                                                                                                             aCCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGACCTGCACCGGCTGTTCCCCCAGCAGTTGGCTGTc     >  3:306861/1‑149 (MQ=255)
                                                                                                                                               cAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGACCTGCACCGGCTGTTCCCCCAGCAGTTGGCTGTCGa   >  6:91877/1‑149 (MQ=255)
                                                                                                                                                aGGCGAGCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGGAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGACCTGCACCGGCTGTTCCCCCAGCAGTTGGCTGTCGAt  <  6:379952/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GTCTGGGTCAGGTGCAGTTCGGCATGCGCGGCATCGACAAGGGTGATGATGCCGTCGAGGATGTAGCGCTCACGCAGCTCTTCGTCGATGAAGAACGTTTGCGCGACCGGCGCCGGGTCGGCCAGGCCGGTGCATTCGATCACCAGGCGATCGAAGGCGATTTCCCCGGCATCCAGGCGTTCAAGCAGCAGGTACAGCGCGCGAGTCAGGTCGCCATGGATGCTGCAGCATACACAACCGTTGGCCAAGGTCATGACCTGCACCGGCTGTTCCCCCAGCAGTTGGCTGTCGAT  >  NC_002947/5263825‑5264117

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: