Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 786,462 (A)5→4 coding (830/1053 nt) aroG → 3‑deoxy‑7‑phosphoheptulonate synthase, Phe‑sensitive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913786,4580A.98.3% 255.0 / ‑3.3 60coding (826/1053 nt)aroG3‑deoxy‑7‑phosphoheptulonate synthase, Phe‑sensitive
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (31/28);  minor base G (1/0);  total (32/28)
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.09e-01

CCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACG  >  NC_000913/786319‑786593
                                                                                                                                           |                                                                                                                                        
ccATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGGAGGGCTGAACAAAGCAGGCCTGCCAGCAAAGGTGATGATCGCATACAGCCATGCTAACTCGTCGAAACAATTCagaa                                                                                                                                       >  2:639042/1‑143 (MQ=255)
ccATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCaaaa                                                                                                                                       >  1:598180/1‑143 (MQ=255)
ccATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATACTAACTCGTCCAAAGAATTCagaa                                                                                                                                       >  1:604595/1‑143 (MQ=255)
cactatCATTCTGCGCGGCGGTAACGACCCTAACTACGGCGCGAACCACGTTGCTCAAGTGCAAGAAGGCCTGAACAACGCAGGCCTGCCAGCACAGGCGCTGATCGATCTCACCCATGCTAACTCGTCCAAACAATTCaaac                                                                                                                                       <  2:131135/140‑2 (MQ=255)
    atCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATCCTAACCAGTCCAAACAAGTCaaaagcaa                                                                                                                                   >  2:1010976/1‑139 (MQ=255)
     tCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCaaaagcagc                                                                                                                                  >  2:97373/1‑138 (MQ=255)
        ttCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCGGCCAGCACAGGTGATGATCGATTTCAGCCATGCTCACTCGTCCAAACACTTCAAAAGCAgatgg                                                                                                                              >  2:868239/1‑143 (MQ=255)
          cggcgcGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACCGGTCTCAATCGGAATCAGCCTTGCTAACTCGTCCAAACAATTCAAAAGCAgatggat                                                                                                                            <  1:737114/141‑1 (MQ=255)
            gcgcGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAgatggatga                                                                                                                          >  2:1007372/1‑142 (MQ=255)
                   ggTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGACGGGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATGCAAAAGCAGATGGATGTTTGTGCt                                                                                                                   >  2:1050688/1‑143 (MQ=255)
                     tgaggggCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGa                                                                                                                 <  1:137807/137‑1 (MQ=255)
                       aaGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACg                                                                                                               >  2:668375/1‑143 (MQ=255)
                            ccTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGGTAACTCGTCCAAACAATTCAAAGGCAGATGGATGTTTGTGCTGGCGTGTGc                                                                                                          >  1:935285/1‑143 (MQ=255)
                               aaCTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGCTCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCcag                                                                                                       <  1:653969/143‑1 (MQ=255)
                                  tACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTATGTGCTGACGTTTGCcagcag                                                                                                    >  2:449671/1‑143 (MQ=255)
                                    cAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCAAGCAGAt                                                                                                  >  2:1007008/1‑143 (MQ=255)
                                     aGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCGGACCAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCGAACAATTCAAAAGCAGATGGAGGCTCGTGCCGCCGTTTGCCGGCAGAtt                                                                                                 >  2:192482/1‑143 (MQ=255)
                                       ggcaaaGCACGTTTCCGAAGTGAAAGAAGGGCTGAACAAAGCAGGCTTGCCAGCACAGGGGATGATCGATGTCAGGCATGCTAACTCGGCCAATCAATTCAAAAGCAGATGGATGTTTGTGTTGACGTTTGCCAGCAGATTGc                                                                                               <  2:938111/139‑1 (MQ=255)
                                       cgcgAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGc                                                                                               >  1:299947/1‑143 (MQ=255)
                                            aGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGACATCAGATTGCCGGTg                                                                                          >  1:951805/1‑143 (MQ=255)
                                             gCACGGTGCTGAAGGGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTgg                                                                                         <  2:1057828/143‑1 (MQ=255)
                                                   tGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGaaaa                                                                                   <  1:919700/143‑1 (MQ=255)
                                                      tGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGc                                                                                >  1:683138/1‑143 (MQ=255)
                                                          gTGAAAGAAGGGCTGAACAAAGCAGACCTGCCAGCACAAGTGGTGGTCGATTTCAGCCATGCTAACTCGTGCAACCAATACAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGACGATTGCAGTTGGCGCACAGGcgact                                                                            >  1:599658/1‑139 (MQ=255)
                                                           tGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGGGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGCTGGATGTTTGTGCTGACGGTTGCCAGCAGAGTGCCGGTGGCGAAGAGGccgata                                                                           >  2:451395/1‑139 (MQ=255)
                                                             aaaGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCattatt                                                                         >  1:609458/1‑143 (MQ=255)
                                                              aagaagGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTg                                                                        <  1:658586/143‑1 (MQ=255)
                                                                 aagGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCg                                                                    >  1:609217/1‑143 (MQ=255)
                                                                     gcgtaccAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCTGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCTAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAt                                                                <  2:877998/137‑1 (MQ=255)
                                                                     gCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAt                                                                >  1:646888/1‑143 (MQ=255)
                                                                     gCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGACGGCAGTTTGTGCCGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAt                                                                <  1:588137/143‑1 (MQ=255)
                                                                         caCAAAGCAGGCCTGCCAGCACAGCTGATGATCGATTTCGGCCAGGCTAACTCGTCCAAACAAGTCGAAAGCAGATGGATGCGTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAtggtg                                                            <  2:521550/142‑1 (MQ=255)
                                                                         aaCAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAACGGCCATTATTGGCGTGAtggtg                                                            >  2:961348/1‑143 (MQ=255)
                                                                         aaCAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAtggtg                                                            >  2:353227/1‑143 (MQ=255)
                                                                         aaCAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCACTGCTGGAAAGACCAAATGATTAGTTAAACTTTCTTCGGTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAtggtg                                                            >  1:788167/1‑143 (MQ=255)
                                                                             gaGCAGGTCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATGCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGaaa                                                        <  1:190636/142‑1 (MQ=255)
                                                                             aggCAGCCTTGTCAGCACAGGTGATCATCGATTTCAGGCATGCTACCTCGTCCTAACAATTCAAAAGCAGATGTGTGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGaaa                                                        <  2:96258/141‑1 (MQ=255)
                                                                             aaGTAGGCCTGCCAGCACAGGTCATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAACAGCAGCTGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGaaa                                                        <  2:569193/143‑1 (MQ=255)
                                                                             aaGCAGGCCCGCCAGCACAGGGGATGAGCGAGGGCAGCCACGCTAACTCGGCCAAACAATGCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGaaa                                                        <  1:14899/143‑1 (MQ=255)
                                                                                 aGGCCTGCCAGCACAGGTGATGATCGATTTCAGTCATGCTAACTCGTCCAAATAATTCAAAAGCAGATGGATGTTTGTGATGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCa                                                    <  1:82109/143‑1 (MQ=255)
                                                                                  ggCCTGCCAGCACAGGGGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCAt                                                   <  1:579552/143‑1 (MQ=255)
                                                                                     cTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCCTTATTCGCGTGATGGTGGCGAGCCATCtc                                                >  2:1071079/1‑142 (MQ=255)
                                                                                      tGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCtgg                                               <  2:788167/143‑1 (MQ=255)
                                                                                       gaccgCAAAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCtggt                                              <  1:663810/139‑1 (MQ=255)
                                                                                       gCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGGCCAAACAATGCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCtggt                                              <  2:862817/143‑1 (MQ=255)
                                                                                           gCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAAATCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGaa                                          <  1:381768/143‑1 (MQ=255)
                                                                                                     tgatCGATTTCAGCCATGCTAACTCGTCCAGACAATTCAAAAGCAGATGGATGTTTGTGCCGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCaga                                <  2:535945/143‑1 (MQ=255)
                                                                                                      gatCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCAATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCagag                               >  2:725148/1‑143 (MQ=255)
                                                                                                       atCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGc                              >  2:575130/1‑143 (MQ=255)
                                                                                                       atCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGGTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCACTCAGAGc                              >  1:357165/1‑143 (MQ=255)
                                                                                                                aGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATCGACGGTGGCGAAAACGCATTTTAGGCAGCGATAGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGc                     >  2:761089/1‑143 (MQ=255)
                                                                                                                aGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGGTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGc                     >  1:407279/1‑143 (MQ=255)
                                                                                                                 gCCATGCTAACTCGTCCAAACAATTCAAAAGCAGAGGGATGGTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCg                    <  2:852535/143‑1 (MQ=255)
                                                                                                                  cgaTGTTAACTCGGCCAAACAATTCAAAAGGAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCgg                   <  2:877049/141‑1 (MQ=255)
                                                                                                                  ccATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGCGAGCgg                   >  2:611621/1‑143 (MQ=255)
                                                                                                                   cATGCTATGTCGTCCAGACAACTCATAAGCAGATGGATGTTTGTGCTCACGTTTGCCAGGAGGTGGCCGGTGGCGAAAAGGCCATAATTGGCGTGATGGGTGAGGGTCATATGGTGGAAGGCAATCAGAGCTCCGAGAGCggg                  <  2:281756/143‑1 (MQ=255)
                                                                                                                    aTGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCgggg                 <  1:562257/143‑1 (MQ=255)
                                                                                                                      gCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTACTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAg               >  1:644753/1‑143 (MQ=255)
                                                                                                                         agcTCGGCCAAACAATTCGAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCg            <  2:249317/141‑1 (MQ=255)
                                                                                                                              gTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGc       <  1:1050688/143‑1 (MQ=255)
                                                                                                                              gTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGACCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCGCGAGAGCGGGGAGCCGCTGGc       >  2:437312/1‑143 (MQ=255)
                                                                                                                               tCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGcc      >  2:1019164/1‑143 (MQ=255)
                                                                                                                                ccAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAGGGTGGAAAGCCATCAGGTGGCCGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCGGGCCt     >  2:617603/1‑143 (MQ=255)
                                                                                                                                   aaCAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCAGTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGTCCTACg  >  2:859073/1‑143 (MQ=255)
                                                                                                                                      aaTTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGaa                                          >  2:304405/1‑100 (MQ=255)
                                                                                                                                      aaTTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGaa                                          <  1:304405/100‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                        
CCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACG  >  NC_000913/786319‑786593

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