Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,817,957 T→A 8.8% D66E (GAT→GAA dinD → DNA damage‑inducible protein D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,817,9570TA 8.8% 149.6 / 5.7 62D66E (GAT→GAAdinDDNA damage‑inducible protein D
Reads supporting (aligned to +/- strand):  ref base T (25/29);  new base A (2/4);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

CAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCC  >  NC_000913/3817815‑3818099
                                                                                                                                              |                                                                                                                                              
cAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCATGTGAGAACAGCGGACACACAATTGATGACCATTTCGAGGAGATCCTCGat                                                                                                                                                >  1:884651/1‑143 (MQ=255)
cAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTGATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGCACCTCGat                                                                                                                                                >  2:199101/1‑143 (MQ=255)
  ggaaaaTGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTAGACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTTTGACCATATCGAGGATCTCCTCGAAATg                                                                                                                                             <  1:839711/139‑1 (MQ=38)
   aaCAATGGTCAGCAAGACAACTGTGGAAAGTGCTGGTCGACACAAAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATg                                                                                                                                             <  2:628445/143‑1 (MQ=255)
     cAATGGACAGCAAGACAACCTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATGGt                                                                                                                                           <  1:490565/143‑1 (MQ=255)
        tGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCaa                                                                                                                                        >  2:771642/1‑143 (MQ=255)
        cgggcAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATGGTCat                                                                                                                                        <  1:31021/139‑2 (MQ=255)
         ggTCAGCAATACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATGGTCatc                                                                                                                                       <  1:345585/143‑3 (MQ=39)
          gTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGCCaaaa                                                                                                                                      >  2:709384/1‑143 (MQ=255)
           tCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCCAAGAAGCCTGTGAGAACAGTGGTCAGACAATTGATGACCATTTTGAAGAGATCCGCGAAATGATAAAAAt                                                                                                                                     >  2:374233/1‑143 (MQ=255)
             aGCCAGACGCCTTGGGACACTACGGGGTTATTCAGAGTATCGTCACTTTATACCTGTAGTCACGCGCGCCAACGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTg                                                                                                                                   <  2:808129/143‑1 (MQ=255)
               cAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGc                                                                                                                                 >  2:790089/1‑143 (MQ=255)
               cAAGACAACTTGGGAAACGACTGGGTTATTCAGAGTATCGTCACTTTCTACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGc                                                                                                                                 <  1:5303/143‑1 (MQ=255)
                     aaCTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGCTTCTCGATATGGTCAAAATTGGCTCAAAt                                                                                                                           >  1:16878/1‑143 (MQ=255)
                        ttGGTAAACTACTTGGTTATTCAGAGTAGCGTCACTTTATACCTGTATTAACGCGCGCCAATGAAGCCTTTGAGAACAGTGGTCACACTATTGATGACCATTTCGATGAGATCCTCGATATGGTCAAAATTGGCTCAAATGcc                                                                                                                        <  2:40950/143‑1 (MQ=255)
                        ttGGGAAACTCCTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGACTGCGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGcc                                                                                                                        >  1:189564/1‑143 (MQ=255)
                            gAAACGACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCaaaa                                                                                                                    <  1:156999/143‑1 (MQ=255)
                                 tactGGATTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCa                                                                                                               <  1:10012/143‑1 (MQ=255)
                                      ggTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaa                                                                                                          >  2:559022/1‑143 (MQ=255)
                                       gTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaaa                                                                                                         <  1:269274/143‑1 (MQ=255)
                                            tCAGAGTATCGTCACTTTATACCTGTCTTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATAGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAACAGAGCATTAAAAGACAt                                                                                                    >  2:643628/1‑143 (MQ=255)
                                            tCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGGACAGTGGTGAGACAATTGATGACCATTTAGAGGAGATGCTCGATATGGTCTGAGTTGGCTGAAATGCCACAGGAGGAGTAGAAGAcgg                                                                                                    >  1:174784/1‑141 (MQ=255)
                                            tCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACAt                                                                                                    >  1:567646/1‑143 (MQ=255)
                                             gagagTATCGTCAATTTATACCTGTATTGACGGCCGCCAAAGCAGTCTGTGAGAACAGTGGTCACACACTTGATCTCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATc                                                                                                   <  2:268631/142‑1 (MQ=255)
                                              agagTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACCCAAGTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAACTGCCAAAAGGGCATTAAAGGGCATCg                                                                                                  >  1:10659/1‑143 (MQ=255)
                                                agTATCGTCACTTGATACCTGTATTAACGCGCGCCAAAGAAGCCTGGGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCGCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTa                                                                                                <  2:643339/143‑1 (MQ=255)
                                                    tCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTActct                                                                                            >  1:488817/1‑143 (MQ=255)
                                                     cGTCACTGGCTACCTGTAGTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACtctc                                                                                           <  2:480994/143‑1 (MQ=255)
                                                       tCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTccc                                                                                         <  1:771642/143‑1 (MQ=255)
                                                       tCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGAAATGGACAAAACTGAATCAAATGCCAAAAAAGCAGTAAAAGACATCGTACTCTccc                                                                                         >  2:382843/1‑143 (MQ=255)
                                                           tttATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTa                                                                                     >  2:409704/1‑143 (MQ=255)
                                                               tatcTGTATTATCGGGCGCCAAAGTGCCCTATGAGAACAGTGGTCGCAGACTTGATGTCGATTTATAGGCGACCCCCTTTATTGTCGAACTTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGcc                                                                                 <  2:826720/140‑1 (MQ=255)
                                                                aCCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCt                                                                                <  1:304765/143‑1 (MQ=255)
                                                                 ccTGTATTGACGCGCGCCACAGCACCCTGTGACGACGTTGGTCACACAATTGATGACCTTTGTGAGGAGATCCGCGATATGGTCAAAGCTTGCTCACATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTg                                                                               <  2:656963/143‑1 (MQ=255)
                                                                  cTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGt                                                                              >  1:300901/1‑143 (MQ=255)
                                                                   tGTATTAACGCGCGCCAAAGGGGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGtt                                                                             <  1:120104/143‑1 (MQ=255)
                                                                    gTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTa                                                                            >  2:296148/1‑143 (MQ=255)
                                                                        tAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTg                                                                        <  1:496298/143‑1 (MQ=255)
                                                                           cgcgcgCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTACGCCTGTTACCTGgta                                                                     >  2:635980/1‑143 (MQ=255)
                                                                            gcgcgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCACGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCAGTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgcag                                                                    >  2:21814/1‑143 (MQ=255)
                                                                               agcCAAAGAAGGCTGTGACCACAGTGGTCACACAATTGATGCCCGTTGCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTAc                                                                 <  1:374233/142‑1 (MQ=255)
                                                                                 ccAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACaa                                                               >  1:716560/1‑143 (MQ=255)
                                                                                  cAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACaaa                                                              >  1:106518/1‑143 (MQ=255)
                                                                                   aaaGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACaaaa                                                             >  1:147142/1‑143 (MQ=255)
                                                                                      gaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACgg                                                          >  2:563626/1‑143 (MQ=255)
                                                                                      gaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACgg                                                          >  2:129888/1‑143 (MQ=255)
                                                                                      gaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACgg                                                          <  1:648753/143‑1 (MQ=255)
                                                                                           cTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGAcc                                                     <  2:361740/143‑1 (MQ=255)
                                                                                               gagaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGCAAATACAAAACGGCGACCCTGc                                                 >  2:570161/1‑143 (MQ=255)
                                                                                                 gaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGa                                               <  1:296148/143‑1 (MQ=255)
                                                                                                              acaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGCCATTGc                                  >  2:503937/1‑143 (MQ=255)
                                                                                                              acaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGGTACCTGGTAGTACCAAACGGCGACCCTGCGAAACCGGTCATTGc                                  >  2:643280/1‑143 (MQ=255)
                                                                                                                    gctggccATTTCGAGGAGAGCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATCAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCcgcagg                            <  2:742822/138‑1 (MQ=255)
                                                                                                                       ggccATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagggca                         <  2:884651/141‑1 (MQ=255)
                                                                                                                         ccATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCGCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGa                       <  2:189564/143‑1 (MQ=255)
                                                                                                                                gaggagATCCTCGAGATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAttt                <  2:573224/143‑1 (MQ=255)
                                                                                                                                  gtagaTCCTCGAGATGTGCAAAATTGCCTCAAATGCCAAAAGAGCATTAAAAGACTTCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTg              <  1:753610/141‑1 (MQ=255)
                                                                                                                                  ggagATCCTCGATATGGTCAACCTCGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTCCTCCCCCGCTCTGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTg              <  1:12120/143‑1 (MQ=255)
                                                                                                                                  ggagATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTg              <  2:488817/143‑1 (MQ=255)
                                                                                                                                  ggagATCCTCGATATGGTCAAAATTGGCTCAAAGGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTg              <  1:89186/143‑1 (MQ=255)
                                                                                                                                    ggaTCCTCGATATGGTCAAAATTAGCTCAAATGCCAAACGAGCATTAAAAGACATCGTACTCTCCCACTACGCCTGGTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCt            <  2:778083/142‑1 (MQ=255)
                                                                                                                                         cTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCa       >  1:468283/1‑143 (MQ=255)
                                                                                                                                            gATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACGCTCCCGCTATGCCTGCCCCCGGGTGGGACGAGACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGAc    <  2:356628/143‑1 (MQ=255)
                                                                                                                                            gATATGGTCAAAATTGGCGCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGAc    <  1:340528/143‑1 (MQ=255)
                                                                                                                                             atatGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACATACCTTGCTACACAGAcc   >  2:655959/1‑143 (MQ=255)
                                                                                                                                              tatGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAACAGACATCGTACTCTCCCGCTATGCCTGTTAGATGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCGGGGCAGACGTATTTGGCGATCAAGActc  >  2:813013/1‑141 (MQ=255)
                                                                                                                                              |                                                                                                                                              
CAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCC  >  NC_000913/3817815‑3818099

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: