Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,763,869 (T)7→8 coding (1285/1617 nt) yfjL ← CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,763,8621.T98.6% 237.8 / ‑3.7 69coding (1292/1617 nt)yfjLCP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (30/38);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

ATTAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑A‑TTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCACCTT  >  NC_000913/2763735‑2763991
                                                                                                                                 |                                                                                                                                 
aTTAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑A‑ttttttttgcg                                                                                                                        <  2:476600/139‑5 (MQ=255)
  tAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTgcgcg                                                                                                                       >  2:504060/1‑139 (MQ=255)
   agagTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTCATGTCGAGG‑ATTTTTTTTgcgcgc                                                                                                                      >  1:322870/1‑139 (MQ=255)
       tCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTgcgcg                                                                                                                       >  2:746471/1‑134 (MQ=255)
       tCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTgcgcg                                                                                                                       <  1:746471/134‑1 (MQ=255)
        cGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTAt                                                                                                                 <  1:127341/139‑1 (MQ=255)
         gTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATc                                                                                                                >  1:564384/1‑139 (MQ=255)
          tAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCg                                                                                                               >  1:1093433/1‑139 (MQ=255)
             cccACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATg                                                                                                            <  2:390244/139‑1 (MQ=255)
              ccACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGc                                                                                                           <  1:1061387/139‑1 (MQ=255)
               cACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCa                                                                                                          >  2:939514/1‑139 (MQ=255)
                       gaacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGaaa                                                                                                  <  2:443362/139‑1 (MQ=255)
                          caAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGCT‑TTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGt                                                                                               <  2:268005/139‑1 (MQ=255)
                           aaaCGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTc                                                                                              <  1:839513/139‑1 (MQ=255)
                            aaCGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                             <  1:504060/139‑1 (MQ=255)
                            aaCGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                             <  1:882735/139‑1 (MQ=255)
                            aaCGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                             <  2:51918/139‑1 (MQ=255)
                                  aaaCGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTAc                                                                                       >  1:922506/1‑139 (MQ=255)
                                                  cATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACt                                                                                   <  1:206538/127‑1 (MQ=255)
                                                  cATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACt                                                                                   >  2:206538/1‑127 (MQ=255)
                                                     cGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTgcgcg                                                                                                                       >  2:932629/1‑88 (MQ=255)
                                                     cGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTgcgcg                                                                                                                       <  1:932629/88‑1 (MQ=255)
                                                       aaaGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGc                                                                  <  2:84924/139‑1 (MQ=255)
                                                        aaGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCa                                                                 >  2:460968/1‑139 (MQ=255)
                                                              ggCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                             <  1:329071/105‑1 (MQ=255)
                                                              ggCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGa                                                           >  2:721266/1‑139 (MQ=255)
                                                              ggCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATATcc                                                                                             >  2:329071/1‑105 (MQ=255)
                                                               gCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGAc                                                          <  1:235963/139‑1 (MQ=255)
                                                               gCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGAc                                                          <  1:876293/139‑1 (MQ=255)
                                                                  gCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCTCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAACCTGACGGc                                                       <  2:158393/139‑1 (MQ=255)
                                                                   ccATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                      <  2:564384/139‑1 (MQ=255)
                                                                     aTCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAAt                                                    <  1:94625/139‑1 (MQ=255)
                                                                                cGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGaaa                                                                                                  >  1:758455/1‑82 (MQ=255)
                                                                                cGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGaaa                                                                                                  <  2:758455/82‑1 (MQ=255)
                                                                                cGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAatat                                            >  1:210282/1‑136 (MQ=255)
                                                                                cGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAatat                                            <  2:210282/136‑1 (MQ=255)
                                                                                cGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAg                                         >  2:542362/1‑139 (MQ=255)
                                                                                  tCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTc                                       >  1:1124267/1‑139 (MQ=255)
                                                                                     aTGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCAcc                                    >  2:343142/1‑139 (MQ=255)
                                                                                     aTGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCAcc                                    >  2:1017604/1‑139 (MQ=255)
                                                                                      tGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                      >  2:762343/1‑120 (MQ=255)
                                                                                      tGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                      <  1:762343/120‑1 (MQ=255)
                                                                                       gTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCaa                                  <  2:65476/139‑1 (MQ=255)
                                                                                        tCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAAc                                 <  1:939514/139‑1 (MQ=255)
                                                                                           tGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTGATATCAGTCACCAACTCg                              >  1:756275/1‑139 (MQ=255)
                                                                                                  gCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGc                       <  2:52096/139‑1 (MQ=255)
                                                                                                  gCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGAGCCGTTATTTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGc                       <  2:979512/139‑1 (MQ=255)
                                                                                                  gCAGGTAGTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATACGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGc                       >  1:801499/1‑139 (MQ=255)
                                                                                                         gTTATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCa                <  2:1124267/139‑1 (MQ=255)
                                                                                                          ttATCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCaa               <  2:922506/139‑1 (MQ=255)
                                                                                                            aTCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGt                                                                                               >  1:575543/1‑57 (MQ=255)
                                                                                                            aTCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGt                                                                                               <  2:575543/57‑1 (MQ=255)
                                                                                                            aTCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                      <  2:1007167/98‑1 (MQ=255)
                                                                                                            aTCGTCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                      >  1:1007167/1‑98 (MQ=255)
                                                                                                                tCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGa         >  1:986008/1‑139 (MQ=255)
                                                                                                                tCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGa         >  1:96287/1‑139 (MQ=255)
                                                                                                                tCTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGa         >  2:266476/1‑139 (MQ=255)
                                                                                                                tCTTTAATGTCCAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAAGCCGAAATGGCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGa         <  2:276498/139‑1 (MQ=255)
                                                                                                                 cTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGAt        <  1:330163/139‑1 (MQ=255)
                                                                                                                 cTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGAt        <  1:266476/139‑1 (MQ=255)
                                                                                                                 cTTTAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATG‑CCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATc       <  1:38307/139‑1 (MQ=255)
                                                                                                                  tttAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCa                <  1:495515/130‑1 (MQ=255)
                                                                                                                  tttAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCa                >  2:495515/1‑130 (MQ=255)
                                                                                                                  tttAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATc       <  2:756275/139‑1 (MQ=255)
                                                                                                                  tttAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATc       <  2:1093433/139‑1 (MQ=255)
                                                                                                                  tttAATGTCGAGG‑AGTGTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTGAATATCAGTGACCAACTCGCTGAAGCCAATCCAAAGCTGATc       <  2:801499/139‑1 (MQ=255)
                                                                                                                   ttAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCa      <  2:96287/139‑1 (MQ=255)
                                                                                                                   ttAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCa      <  2:986008/139‑1 (MQ=255)
                                                                                                                    tAATGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCAc     >  2:705530/1‑139 (MQ=255)
                                                                                                                       tGTCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCACCtt  >  2:886797/1‑139 (MQ=255)
                                                                                                                         tCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                      <  1:391691/85‑1 (MQ=255)
                                                                                                                         tCGAGG‑ATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                      >  2:391691/1‑85 (MQ=255)
                                                                                                                                 |                                                                                                                                 
ATTAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGG‑A‑TTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCACCTT  >  NC_000913/2763735‑2763991

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑