Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,130,849 C→A 6.5% Q339K (CAA→AAA)  metL → Bifunctional aspartokinase/homoserine dehydrogenase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,130,8490CA 6.5% 269.2 / 4.0 80Q339K (CAA→AAA) metLBifunctional aspartokinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base A (4/2);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCG  >  NC_000913/4130703‑4130991
                                                                                                                                                  |                                                                                                                                              
cGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCAAAAAGAGATCGAcc                                                                                                                                                >  7:257647/1‑147 (MQ=255)
            gTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcg                                                                                                                                   <  5:285920/148‑1 (MQ=255)
            gTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcg                                                                                                                                   >  4:192361/1‑148 (MQ=255)
              tCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgcg                                                                                                                                 >  5:292352/1‑148 (MQ=255)
                 aCGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACgcgcgc                                                                                                                              <  8:168177/148‑1 (MQ=255)
                      cATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTa                                                                                                                         >  4:124320/1‑148 (MQ=255)
                        ttGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACg                                                                                                                       <  4:1926/148‑1 (MQ=255)
                         tGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTTTTCTCACCAGCCACGATCATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACAAAATCCAGAAACGCGCGCAGGTACGc                                                                                                                      >  1:149319/1‑148 (MQ=255)
                         tGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGc                                                                                                                      >  3:99624/1‑148 (MQ=255)
                         tGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACAAAATCCTGAAACGCGCGCAGGTACGc                                                                                                                      >  1:154836/1‑148 (MQ=255)
                                  gCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCg                                                                                                             >  2:20324/1‑148 (MQ=255)
                                          ccGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGCTCGACCAAATc                                                                                                                                           <  8:171789/110‑1 (MQ=255)
                                          ccGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATc                                                                                                                                           >  7:171789/1‑110 (MQ=255)
                                          ccGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGaa                                                                                                                                      >  1:170999/1‑115 (MQ=255)
                                          ccGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGaa                                                                                                                                      <  2:170999/115‑1 (MQ=255)
                                                   gTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACg                                                                                             <  1:20324/147‑1 (MQ=255)
                                                     gcgcGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATc                                                                                          >  2:355373/1‑148 (MQ=255)
                                                      cgcgTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCg                                                                                         >  3:144732/1‑148 (MQ=255)
                                                       gcgTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTTGCGTACATAACGATCGc                                                                                        >  5:128540/1‑148 (MQ=255)
                                                           aTTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTg                                                                                                         >  1:272051/1‑127 (MQ=255)
                                                           aTTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTg                                                                                                         <  2:272051/127‑1 (MQ=255)
                                                           aTTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACAGAACGATCGCCAGt                                                                                    >  8:248630/1‑148 (MQ=255)
                                                              gTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTtgc                                                                                 >  4:111100/1‑148 (MQ=255)
                                                               tCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTtgct                                                                                >  5:368118/1‑148 (MQ=255)
                                                                 aCCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTtgctgc                                                                              <  3:111100/148‑1 (MQ=255)
                                                                  ccAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTgg                                                                                                        >  1:159513/1‑121 (MQ=255)
                                                                  ccAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTgg                                                                                                        <  2:159513/121‑1 (MQ=255)
                                                                    aGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAAt                                                                           >  3:330027/1‑148 (MQ=255)
                                                                     gCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCAAAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAAtt                                                                          <  8:257647/148‑1 (MQ=255)
                                                                        aCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCAGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTg                                                                       >  1:204910/1‑148 (MQ=255)
                                                                         agatgatGTCTGTTTGATTGTTTTTCAGGTGCCCGCAAGTCATTATTTAAAACTGGCTCATAAAGAGATCGAACAAATCCTTAAACGCGCTCATGTTCGCCCTCTGGCGGTTGGCGTACATAAAGATCGCCATTTGCTGCAATTTTGc                                                                      <  4:330027/147‑1 (MQ=255)
                                                                          gatgatGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCt                                                                     <  7:248630/148‑1 (MQ=255)
                                                                            tgatgTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTa                                                                    >  7:136599/1‑147 (MQ=255)
                                                                                gTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCa                                                                                                                             <  5:22514/86‑1 (MQ=255)
                                                                                gTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCa                                                                                                                             >  6:22514/1‑86 (MQ=255)
                                                                                     ttagattgaGTTTCAGGTGCCCGCCAGTCAGGATATCCAACTGGCGCATAAAGAGATCCACAAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTgaga                                                           <  2:149319/144‑1 (MQ=255)
                                                                                     tgtgattgaGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACAAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGa                                                           <  2:154836/145‑1 (MQ=255)
                                                                                        gattgaGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGt                                                        <  8:136599/147‑1 (MQ=255)
                                                                                              gTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGAc                                                 <  4:144732/148‑1 (MQ=255)
                                                                                                ttCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGa                                                  >  1:296201/1‑145 (MQ=255)
                                                                                                 tCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGAc                                                 <  2:296201/145‑1 (MQ=255)
                                                                                                    ggTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTgcg                                           >  5:219635/1‑148 (MQ=255)
                                                                                                      tGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGAc                                                 >  1:356786/1‑140 (MQ=255)
                                                                                                      tGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGAc                                                 <  2:356786/140‑1 (MQ=255)
                                                                                                             caggcaggATTTCAAACTGGCGCATAAAGAGATCGACCACATCCTGACACGCGCGCCGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATc                                  <  6:292352/148‑1 (MQ=255)
                                                                                                               gTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCt                                <  7:220912/148‑1 (MQ=255)
                                                                                                                tCATGATTTAAAACTTGTGCATAACGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGc                                                     <  4:254661/126‑1 (MQ=255)
                                                                                                                tCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTc                               >  2:17009/1‑148 (MQ=255)
                                                                                                                 cAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCATGTACGCCCGCTGGCGGTTGGCGTACATAACGAGCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGcc                                                    >  3:254661/1‑126 (MQ=255)
                                                                                                                 cAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTc                               >  5:198798/1‑147 (MQ=255)
                                                                                                                 cAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTc                               >  7:26487/1‑147 (MQ=255)
                                                                                                                  aGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTcga                             >  2:301947/1‑148 (MQ=255)
                                                                                                                     aTTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTcgacg                           >  4:164715/1‑147 (MQ=255)
                                                                                                                       ttCAAACTGGCGCATAAAGAGATCGACAAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGATGGCGTACATAACGATCGCCAGATGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAg                        >  3:5084/1‑148 (MQ=255)
                                                                                                                        tctaaCTGGCGCTTAAAGTGATCGTCCAAATCCTGTAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATGGCCAGTTGCTTCAATTTTGCTACACCTCAGAAGTGGCCGACa                                                <  4:107709/120‑1 (MQ=255)
                                                                                                                         cAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTg                                                                       >  3:97200/1‑99 (MQ=255)
                                                                                                                         cAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTg                                                                       <  4:97200/99‑1 (MQ=255)
                                                                                                                         cAAACTGGCGCATAAAGAGATCGACAAAATCCTGAAACGCGCGCAGGTACGCACGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACAACTCAGAAGTGGCCGACAg                                               >  3:107709/1‑123 (MQ=255)
                                                                                                                          aaaCTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCg                      <  3:192361/147‑1 (MQ=255)
                                                                                                                          aaaCTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCg                      <  3:164715/147‑1 (MQ=255)
                                                                                                                          aaaCTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCg                      <  3:124320/147‑1 (MQ=255)
                                                                                                                          aaaCTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCg                      <  1:301947/147‑1 (MQ=255)
                                                                                                                           aaCTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCggg                    <  6:198798/148‑1 (MQ=255)
                                                                                                                                   cATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCt                                                                     >  1:383518/1‑91 (MQ=255)
                                                                                                                                   cATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCt                                                                     >  5:245063/1‑91 (MQ=255)
                                                                                                                                   cATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCt                                                                     <  2:383518/91‑1 (MQ=255)
                                                                                                                                   cATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCt                                                                     <  6:245063/91‑1 (MQ=255)
                                                                                                                                     ttaaGAGATCGACCAATTCCGGAAACGCGCGCAGGTACGCCCGCTGGCGGTTTGCGTACATAACGATCGCCAGTTGGTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGc          <  4:5084/146‑1 (MQ=255)
                                                                                                                                      aaaGAGATCGACCAAATCTTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTa                <  8:321880/141‑1 (MQ=255)
                                                                                                                                      aaaGAGATCGACCAAATCTTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTa                >  7:321880/1‑141 (MQ=255)
                                                                                                                                      aaaGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGc          >  3:340857/1‑147 (MQ=255)
                                                                                                                                      aaaGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGc          >  8:159328/1‑147 (MQ=255)
                                                                                                                                       aaGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTTGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGa        <  6:128540/148‑1 (MQ=255)
                                                                                                                                        agagATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTac                                                                   <  7:48363/88‑1 (MQ=255)
                                                                                                                                        agagATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTac                                                                   >  8:48363/1‑88 (MQ=255)
                                                                                                                                        agagATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGaa                                                          <  3:320758/97‑1 (MQ=255)
                                                                                                                                        agagATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGaa                                                          >  4:320758/1‑97 (MQ=255)
                                                                                                                                        agagATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGaa       <  6:219635/148‑1 (MQ=255)
                                                                                                                                             tCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCCCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTgcg  >  4:309001/1‑148 (MQ=255)
                                                                                                                                              cGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGaaa                                     >  2:365804/1‑112 (MQ=255)
                                                                                                                                              cGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGaaa                                     <  1:365804/112‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCG  >  NC_000913/4130703‑4130991

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