Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,614,970 C→T 5.9% A304A (GCG→GCA yfgO ← putative UPF0118 family inner membrane permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,614,9700CT 5.9% 217.5 / 2.9 68A304A (GCG→GCAyfgOputative UPF0118 family inner membrane permease
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base T (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

GCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCG  >  NC_000913/2614830‑2615114
                                                                                                                                            |                                                                                                                                                
gcgATTTGCCCATCGGGCCAGGCATGAATCACTGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATAACCACCGATTAAATAATCACCAGCGGATGCAGGTTAAACGCTTCg                                                                                                                                            >  6:222067/1‑147 (MQ=255)
gcgATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCACNACAAACCACCGAAGCTCACCACCGATAANATAATCACCAGCGGATGCAGGTTAACCGCGTCg                                                                                                                                            >  6:216525/1‑147 (MQ=255)
       gCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGaaca                                                                                                                                     <  5:279440/147‑1 (MQ=255)
         ccATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAacaaca                                                                                                                                  >  5:90948/1‑148 (MQ=255)
                   aGGCATGAATCACGGCTTTGATCAGCGGAGCCAACGGAATGGCGAAGAATTTTTTTTTTAATGATAAGTATACAACGTAGATCTACACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGgaga                                                                                                                                        <  2:18009/132‑1 (MQ=255)
                    ggCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAAc                                                                                                                       <  5:222067/148‑1 (MQ=255)
                    ggCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGAGCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAAc                                                                                                                       <  5:216525/148‑1 (MQ=255)
                      cATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAg                                                                                                                     >  1:107041/1‑148 (MQ=255)
                        tGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGc                                                                                                                             <  6:313382/138‑1 (MQ=255)
                        tGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGc                                                                                                                             >  5:313382/1‑138 (MQ=255)
                        tGAATCACGGCTTTGATCAGCGTAGCCAACGGAATCGCTAAGAATACGCCCCAGAATCCCCACAAACAACCGAAGATCACCACCGATAAAATAATCACCAGGGGATGCAGGTTAAcc                                                                                                                                                  <  5:387900/117‑1 (MQ=255)
                         gAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGCATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGc                                                                                                                             >  5:303176/1‑137 (MQ=255)
                         gAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGgag                                                                                                                                         <  6:89681/125‑1 (MQ=255)
                         gAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGgag                                                                                                                                         >  5:89681/1‑125 (MQ=255)
                         gAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGc                                                                                                                             <  6:303176/137‑1 (MQ=255)
                          aaTCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTg                                                                                                                 >  6:71198/1‑148 (MQ=255)
                          aaTCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTg                                                                                                                 >  7:214463/1‑148 (MQ=255)
                              aCGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGt                                                                                                             <  2:107041/148‑1 (MQ=255)
                                   tttGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAgc                                                                                                        >  3:24542/1‑148 (MQ=255)
                                   tttGATCAGCGTCGCCAACGGAATGGCGAAGAATAAGCTCAAGAATGCCCACAATACAACGTAGATCAACACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTcc                                                                                                           <  2:31319/145‑1 (MQ=255)
                                    ttGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCa                                                                                                          >  1:31319/1‑145 (MQ=255)
                                     tGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTg                                                                                                                 >  2:254456/1‑137 (MQ=255)
                                     tGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTg                                                                                                                 <  1:254456/137‑1 (MQ=255)
                                       aTCAGCGTCGCCAACGGAATGGCGAAGAATAGGCCCGAGAATCCCCACAAACCCCCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTAGGAGAACAACACCGGCACTaa                                                                                                                        <  4:38621/128‑1 (MQ=255)
                                       aTCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGGTGCCGTCCAGCGcc                                                                                                     <  6:90948/147‑1 (MQ=255)
                                        tCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTaa                                                                                                                        <  4:33744/127‑1 (MQ=255)
                                        tCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTaa                                                                                                                        >  3:33744/1‑127 (MQ=255)
                                         cAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGga                                                                                                                                          >  1:233470/1‑108 (MQ=255)
                                           gCGTCGCCATCGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGaat                                                                                                <  8:214463/148‑1 (MQ=255)
                                               cGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGaataat                                                                                             >  2:62402/1‑147 (MQ=255)
                                                gCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCa                                                                                           <  7:205879/148‑1 (MQ=255)
                                                 ccAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAg                                                                                          <  1:62402/148‑1 (MQ=255)
                                                 ccAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAg                                                                                          <  5:71198/148‑1 (MQ=255)
                                                     cGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCa                                                                                           <  8:173566/143‑1 (MQ=255)
                                                     cGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCa                                                                                           >  7:173566/1‑143 (MQ=255)
                                                        aaTGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAg                                                                                          <  2:306520/141‑1 (MQ=255)
                                                        aaTGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAg                                                                                          >  1:306520/1‑141 (MQ=255)
                                                                  aaTACGCCCCAGAAGCCCCACAAACCACCGAAGATCACCCCCGATAAAATAATCACCAGCGGATGCAGGTTACCCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCATCGCCTGAATAATCAGATACCCCGCGAagcagc                                                                         >  3:273909/1‑148 (MQ=255)
                                                                   aTACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCt                                                                        >  6:169703/1‑148 (MQ=255)
                                                                            aGAATCCCCAGAACCCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTc                                                               <  6:9119/148‑1 (MQ=255)
                                                                            aGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACTGCTTCGGAGAACAACACCGGCACTa                                                                                                                         >  6:180990/1‑90 (MQ=255)
                                                                            aGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACTGCTTCGGAGAACAACACCGGCACTa                                                                                                                         >  6:186571/1‑90 (MQ=255)
                                                                            aGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACTGCTTCGGAGAACAACACCGGCACTa                                                                                                                         <  5:186571/90‑1 (MQ=255)
                                                                            aGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAAtt                                                                >  3:295219/1‑147 (MQ=255)
                                                                            aGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTc                                                               <  8:168125/148‑1 (MQ=255)
                                                                            aGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGATAACTGCTTCGGAGAACAACACCGGCACTa                                                                                                                         <  5:180990/90‑1 (MQ=255)
                                                                              aaTCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGa                                                                                                  <  3:78794/111‑1 (MQ=255)
                                                                              aaTCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGa                                                                                                  >  4:78794/1‑111 (MQ=255)
                                                                              aaTCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAg                                                                                          <  8:170116/119‑1 (MQ=255)
                                                                              aaTCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAg                                                                                          >  7:170116/1‑119 (MQ=255)
                                                                              aaTCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAAtt                                                                >  7:395125/1‑145 (MQ=255)
                                                                               aTCCCCACAACCCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTc                                                               <  8:395125/145‑1 (MQ=255)
                                                                                tCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTa                                                           <  4:295219/148‑1 (MQ=255)
                                                                                     acaAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTgc                                                      >  7:374747/1‑148 (MQ=255)
                                                                                           caccGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCaaa                                                <  5:129835/148‑1 (MQ=255)
                                                                                           caccGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCaaa                                                >  6:129835/1‑148 (MQ=255)
                                                                                             ccGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTACCCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCAAACt                                              >  5:132829/1‑148 (MQ=255)
                                                                                                 aGATCAACACCGATAAAATAATCAACAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCATGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCAATCTgg                                            <  6:132829/146‑1 (MQ=255)
                                                                                                          ccGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCAAACTGGAATAGcgccac                                 <  8:374747/148‑1 (MQ=255)
                                                                                                             aTAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAAtt                                                                <  5:221799/114‑1 (MQ=255)
                                                                                                             aTAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAAtt                                                                >  6:221799/1‑114 (MQ=255)
                                                                                                                   taatCACGAGCGGATGGAGGTTAAGCGCTTCGGAGAACAACACGGGCAGTAAGAGGTTTCCGTCCTGGGCCTGAATAATCAGATACACGGCGATGCAGCTCCAGTATTCAGTACCTGCGCGAAACTGGat                                          <  4:297424/130‑2 (MQ=255)
                                                                                                                    aatCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCAAACTGGAAt                                         >  3:297424/1‑130 (MQ=255)
                                                                                                                      tCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACCCCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCgg                     <  4:24542/148‑1 (MQ=255)
                                                                                                                             cGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGTTTGCCGTCCAGCGCCTGAATATTCAGATACACCGCGAAGCAGCTCCATAATTCAGTACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTc              <  1:45912/148‑1 (MQ=255)
                                                                                                                                       ttAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATg    >  8:231789/1‑148 (MQ=255)
                                                                                                                                       ttAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATg    >  6:250602/1‑148 (MQ=255)
                                                                                                                                          aCCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATgcg  <  7:231789/147‑1 (MQ=255)
                                                                                                                                          aCCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATgcg  <  5:250602/147‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
GCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCGGATGCAGGTTAACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCG  >  NC_000913/2614830‑2615114

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