Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,352,926 C→T G235S (GGT→AGT)  sapD ← antimicrobial peptide ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,352,9260CT98.6% 267.0 / ‑4.1 73G235S (GGT→AGT) sapDantimicrobial peptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (22/50);  total (22/51)
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.38e-01

GTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACCGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGG  >  NC_000913/1352779‑1353072
                                                                                                                                                   |                                                                                                                                                  
gTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGc                                                                                                                                                  <  2:531323/150‑1 (MQ=255)
       aGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGc                                                                                                                                          <  1:1101399/151‑1 (MQ=255)
          gcCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACg                                                                                                                                       <  1:554178/151‑1 (MQ=255)
                    gTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTAt                                                                                                                             <  1:1061230/151‑1 (MQ=255)
                    gTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTAt                                                                                                                             <  1:1140204/151‑1 (MQ=255)
                     tATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATc                                                                                                                            >  2:1138164/1‑151 (MQ=255)
                     tATCGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGGATAAGGATGATGTGGCATCGTAACACACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTAt                                                                                                                             >  2:178529/1‑150 (MQ=255)
                      aTTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATcc                                                                                                                           <  2:949515/151‑1 (MQ=255)
                          aGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGccc                                                                                                                       <  1:1148481/151‑1 (MQ=255)
                          aGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGccc                                                                                                                       >  1:1050054/1‑151 (MQ=255)
                          aGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGccc                                                                                                                       <  1:901840/151‑1 (MQ=255)
                           ggCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACCGCAGTAAAGCACGTTAATTTTATCCGCCCa                                                                                                                      <  2:653022/151‑1 (MQ=255)
                            gCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCAt                                                                                                                     <  1:1101673/151‑1 (MQ=255)
                              gACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTg                                                                                                                   >  2:714318/1‑151 (MQ=255)
                                  tttATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGc                                                                                                                 <  2:286198/149‑1 (MQ=255)
                                  tttATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCt                                                                                                                <  2:1224926/150‑1 (MQ=255)
                                      tGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGc                                                                                                            >  1:674162/1‑150 (MQ=255)
                                        cGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCAtt                                                                                                         <  1:1141320/151‑1 (MQ=255)
                                         ggCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCAttt                                                                                                        <  1:714508/151‑1 (MQ=255)
                                              cgcgCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAg                                                                                                   <  2:700444/151‑1 (MQ=255)
                                               gcgcTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGt                                                                                                  <  2:1173426/151‑1 (MQ=255)
                                                    gCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATgg                                                                                             <  1:1154802/151‑1 (MQ=255)
                                                      cGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCt                                                                                           <  1:930496/151‑1 (MQ=255)
                                                      cGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCt                                                                                           >  2:729817/1‑151 (MQ=255)
                                                       gAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTg                                                                                          >  2:99142/1‑151 (MQ=255)
                                                              ggTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGc                                                                                    <  1:266505/150‑1 (MQ=255)
                                                              ggTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCa                                                                                   <  1:972451/151‑1 (MQ=255)
                                                                tATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGTGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCaa                                                                                  <  2:708276/150‑1 (MQ=255)
                                                                 aTCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCaaa                                                                                 <  2:954343/150‑1 (MQ=255)
                                                                   cgcgcgGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCaaaa                                                                                <  1:225750/149‑1 (MQ=255)
                                                                   cgcgcgGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAAt                                                                               <  1:878647/150‑1 (MQ=255)
                                                                   cgcgcgGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATa                                                                              <  1:123243/151‑1 (MQ=255)
                                                                    gcccccaTCAGCGCCTGGGCATAAGGATGATGTGGCATCGTCACCCACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAAt                                                                               >  2:393833/7‑149 (MQ=255)
                                                                      gcgGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTg                                                                           <  1:452955/151‑1 (MQ=255)
                                                                      gcgGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTg                                                                           >  2:310390/1‑151 (MQ=255)
                                                                            cAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTg                                                                     <  1:270341/151‑1 (MQ=255)
                                                                              gcgcCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGTCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCt                                                                   <  1:65427/151‑1 (MQ=255)
                                                                               cgcCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTg                                                                  <  1:799430/151‑1 (MQ=255)
                                                                               cgcCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTg                                                                  <  2:242636/151‑1 (MQ=255)
                                                                                gcCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGt                                                                 <  1:530363/151‑1 (MQ=255)
                                                                                   tGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGtttt                                                              >  1:565269/1‑151 (MQ=255)
                                                                                    gggTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGttt                                                               <  1:424715/149‑1 (MQ=255)
                                                                                    gggTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTg                                                             >  1:712804/1‑151 (MQ=255)
                                                                                    gggTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTg                                                             <  2:344820/151‑1 (MQ=255)
                                                                                     ggTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTgg                                                            <  2:1097902/151‑1 (MQ=255)
                                                                                      gTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGt                                                           <  2:622296/151‑1 (MQ=255)
                                                                                        ataAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTg                                                         >  1:159947/1‑151 (MQ=255)
                                                                                          aaGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTga                                                        <  1:1151953/150‑1 (MQ=255)
                                                                                          aaGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTg                                                         >  2:229929/1‑149 (MQ=255)
                                                                                            ggATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGAc                                                     <  1:129564/151‑1 (MQ=255)
                                                                                                gatgTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGAcgcg                                                  <  2:782964/150‑1 (MQ=255)
                                                                                                gatgTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGt                                                 >  2:462721/1‑151 (MQ=255)
                                                                                                 atgTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTc                                                <  1:279020/151‑1 (MQ=255)
                                                                                                   gtgGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTca                                               <  2:681340/150‑1 (MQ=255)
                                                                                                        aTCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGc                                         >  1:347837/1‑151 (MQ=255)
                                                                                                        aTCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGc                                         <  1:137004/151‑1 (MQ=255)
                                                                                                                     ccTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTgcg                             <  1:334693/150‑1 (MQ=255)
                                                                                                                        tACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGcc                           <  1:1206045/149‑1 (MQ=255)
                                                                                                                        tACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCt                          <  2:369099/150‑1 (MQ=255)
                                                                                                                         aCTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTgg                        >  1:288136/1‑151 (MQ=255)
                                                                                                                         aCTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTg                         <  2:159925/150‑1 (MQ=255)
                                                                                                                                 cgGTTTCCACTGTTTGACTGCAGAAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGCGGTATTGCTGATTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGttgttg                 >  1:1164783/1‑150 (MQ=255)
                                                                                                                                   gTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCt              >  2:804787/1‑151 (MQ=255)
                                                                                                                                    tttCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTc             >  2:703196/1‑151 (MQ=255)
                                                                                                                                    tttCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTc             >  2:167998/1‑151 (MQ=255)
                                                                                                                                          cTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTg        >  2:271794/1‑150 (MQ=255)
                                                                                                                                            gTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGa       <  2:697015/149‑1 (MQ=255)
                                                                                                                                            gTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAg      <  1:310444/150‑1 (MQ=255)
                                                                                                                                            gTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGt     <  1:723845/151‑1 (MQ=255)
                                                                                                                                              ttGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGtt    <  1:714449/150‑1 (MQ=255)
                                                                                                                                              ttGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGtt    <  1:707338/150‑1 (MQ=255)
                                                                                                                                               tGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTgg  >  2:666829/1‑151 (MQ=255)
                                                                                                                                               tGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTg   <  2:634241/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACCGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGG  >  NC_000913/1352779‑1353072

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: