Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,169,100 A→G 5.6% N130S (AAC→AGC)  yegS → phosphatidylglycerol kinase, metal‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,169,1000AG 5.6% 219.0 / 3.0 71N130S (AAC→AGC) yegSphosphatidylglycerol kinase, metal‑dependent
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base G (2/2);  total (35/36)
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.45e-01

AGGGGGATGACATACCCGCGCTGGGAATTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACATCATTCATGGCTTAATGCGTATGGATACTCTGCAACCGGACCG  >  NC_000913/2168953‑2169247
                                                                                                                                                   |                                                                                                                                                   
aGGGGGATGACATACCCGCGCTGGGTATTTTGACATAAGTAAACGCCATTGATTTTGAAACCAGTGTAGGGATTGCTGAGGCTCTGGATAAGGCGATTAAAATGGAAAATGACTGTGACGGCATTGCGATAGATATGGCGGAGGTcaa                                                                                                                                                     >  4:47529/1‑148 (MQ=255)
  ggggATGACATACCCGCGCTGGGAATTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCCGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTcaaca                                                                                                                                                   >  8:208030/1‑148 (MQ=255)
   gggATGACATACCCGCGCTGGGAATTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCaacaaa                                                                                                                                                 <  1:239017/149‑1 (MQ=255)
    ggATGACATACCCTCGCTGGGAATTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGCCACTGGATAAGGCGCTGAAACTTGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCaacaaac                                                                                                                                                <  5:106335/149‑1 (MQ=255)
     gATGACATACCCGCGCTGGGAATTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCaacaaaca                                                                                                                                               <  7:267571/149‑1 (MQ=255)
          cATACCCGCGCTGGGAATTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCt                                                                                                                                          <  3:17279/149‑1 (MQ=255)
             aCCCGCGCTGGGAATTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGtt                                                                                                                                       >  4:186889/1‑149 (MQ=255)
                 gcgcTGGGAATTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTttat                                                                                                                                   <  7:216700/149‑1 (MQ=255)
                   gcTGGGAATTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTttatta                                                                                                                                 <  3:210783/149‑1 (MQ=255)
                    cTGGGAATTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTaa                                                                                                                                <  8:266446/149‑1 (MQ=255)
                          aTTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGc                                                                                                                          <  3:140150/149‑1 (MQ=255)
                               gCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAg                                                                                                                     >  6:192821/1‑149 (MQ=255)
                                 cATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGc                                                                                                                          <  8:42182/142‑1 (MQ=255)
                                 cATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGc                                                                                                                          >  7:42182/1‑142 (MQ=255)
                                 cATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGc                                                                                                                   >  7:100364/1‑149 (MQ=255)
                                  aTTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCg                                                                                                                  >  5:2104/1‑149 (MQ=255)
                                   ttAGGAACCGCCAATGATTTTGCCACCAGTGTAGGTATTCCTGAGGCACTGGATAAGGAGCAGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGAAATGGAGCAGGTCAACAAACAACCCTGTTTTAGTAATATGGCGACAGGCg                                                                                                                  >  1:265112/1‑148 (MQ=255)
                                   ttAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGCCAGGCg                                                                                                                  >  8:325986/1‑148 (MQ=255)
                                       gAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGAttt                                                                                                             <  2:79892/149‑1 (MQ=255)
                                             ccAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGAcg                                                                                                       >  4:204968/1‑149 (MQ=255)
                                              cAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGAcgc                                                                                                      >  2:15785/1‑149 (MQ=255)
                                               aaTGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGAc                                                                                                        >  7:196192/1‑146 (MQ=255)
                                               aaTGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGAc                                                                                                        <  8:196192/146‑1 (MQ=255)
                                                 tGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTa                                                                                                   <  8:244127/149‑1 (MQ=255)
                                                 tGATTATGCCACCATTGTAGTTATTCCTGAGTCACTGGATTAGGCGCTGAAACTGTCAATTGCCGGGGACGCCATTGCGATAGATATTGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTa                                                                                                   <  8:275214/149‑1 (MQ=255)
                                                    ttttGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCCGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTa                                                                                                <  7:208030/149‑1 (MQ=255)
                                                      ttGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTAc                                                                                               >  7:173218/1‑148 (MQ=255)
                                                      ttGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTAc                                                                                               <  3:204968/148‑1 (MQ=255)
                                                      ttGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTAc                                                                                               <  3:186889/148‑1 (MQ=255)
                                                        gCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACca                                                                                             <  5:155719/148‑1 (MQ=255)
                                                         ccaccaGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTcaac                                                                                                                                                    <  4:135607/92‑1 (MQ=255)
                                                         ccaccaGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACcaca                                                                                           <  1:15785/149‑1 (MQ=255)
                                                                tgtAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGc                                                                                    <  5:216266/149‑1 (MQ=255)
                                                                  tAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAGCAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCg                                                                                  >  8:82645/1‑149 (MQ=255)
                                                                  tAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGcc                                                                                   <  8:100364/148‑1 (MQ=255)
                                                                  tAGGAATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAGCAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCg                                                                                  >  4:320652/1‑149 (MQ=255)
                                                                    gggATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAAGATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCgg                                                                                 >  3:117123/1‑148 (MQ=255)
                                                                     ggATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGa                                                                                >  4:277807/1‑148 (MQ=255)
                                                                     ggATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGAGATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGAGGATTTGGGACGCGTATTGACACAGAAACGACGGa                                                                                >  8:284490/1‑148 (MQ=255)
                                                                      gATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAACCGCCGGaa                                                                               >  8:128956/1‑148 (MQ=255)
                                                                            tGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGAGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTa                                                                         >  2:42676/1‑148 (MQ=255)
                                                                              aGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTaaa                                                                       <  3:277807/148‑1 (MQ=255)
                                                                               ggCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAg                                                                     <  7:290017/149‑1 (MQ=255)
                                                                                         aaGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTggg                                                           >  7:15638/1‑149 (MQ=255)
                                                                                          aGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAGCAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGt                                                          <  7:82645/149‑1 (MQ=255)
                                                                                          aGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAGCAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGt                                                          <  3:320652/149‑1 (MQ=255)
                                                                                           ggCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTggg                                                           >  6:347218/1‑147 (MQ=255)
                                                                                           ggCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTa                                                         <  6:2104/149‑1 (MQ=255)
                                                                                               cTGCAACTGGCAAATGCCGGTGACGCCATTGCGATAGATATGGAGCTGGTAAACAAACAACCCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCTCAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGt                                                     <  3:47529/149‑1 (MQ=255)
                                                                                                  aaaCTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACcaca                                                                                           >  1:307495/1‑108 (MQ=255)
                                                                                                  aaaCTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACcaca                                                                                           <  2:307495/108‑1 (MQ=255)
                                                                                                   aaCTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTctct                                                 <  8:303797/149‑1 (MQ=255)
                                                                                                   aaCTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTctct                                                 >  2:59628/1‑149 (MQ=255)
                                                                                                      tGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTa                                               >  5:214218/1‑148 (MQ=255)
                                                                                                        gCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTAca                                             <  5:192821/148‑1 (MQ=255)
                                                                                                        gCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTAca                                             <  8:15638/148‑1 (MQ=255)
                                                                                                            ttGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTa                                               <  1:198481/142‑1 (MQ=255)
                                                                                                            ttGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTa                                               >  2:198481/1‑142 (MQ=255)
                                                                                                              gCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTAcc                                                                                              >  8:335140/1‑93 (MQ=255)
                                                                                                              gCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTAcc                                                                                              <  7:335140/93‑1 (MQ=255)
                                                                                                               ccGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACAtcattcat                                     >  1:19819/1‑149 (MQ=255)
                                                                                                               ccGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACAtcattca                                      <  7:325986/148‑1 (MQ=255)
                                                                                                                    gACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAACCCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAACCGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACATCCTTCATGGCt                                 >  2:224130/1‑148 (MQ=255)
                                                                                                                    gACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACATCATTCATGGCt                                 >  6:229267/1‑148 (MQ=255)
                                                                                                                     aCGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTa                                                                                                   <  6:48509/81‑1 (MQ=255)
                                                                                                                     aCGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTa                                                                                                   >  5:48509/1‑81 (MQ=255)
                                                                                                                           ttGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGTCGCGTATTACCACAGAATCGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACATCATTCATGGCTTAATGCGt                         <  4:123968/149‑1 (MQ=255)
                                                                                                                              ggatagattTGGCGCAGGTAAACAAAGTAACGTGTTTTATTAATATGGGGTCTGGCGGATTTGGTTGGGGTTTTACCACAGAATCGCGGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACATCATTCATGGCTTAATGCGTAt                       <  4:117123/147‑1 (MQ=255)
                                                                                                                              cGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCgcgc                                                               >  4:30890/1‑108 (MQ=255)
                                                                                                                              cGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCgcgc                                                               <  3:30890/108‑1 (MQ=255)
                                                                                                                               gatagataTGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACATCATTCATGGCTTAATGCGTATg                      >  5:130003/1‑148 (MQ=255)
                                                                                                                               gatagataTGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACATCATTCATGGCTTAATGCGTATg                      >  7:65340/1‑148 (MQ=255)
                                                                                                                                atagataTGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACATCATTCATGGCTTAATGCGTATGGa                    >  1:366822/1‑149 (MQ=255)
                                                                                                                                    atatGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACATCATTCATGGCTTAATGCGTATGGATAct                >  4:10558/1‑149 (MQ=255)
                                                                                                                                     tatGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACATCATTCATGGCTTAATGCGTATGGATActc               <  8:252052/149‑1 (MQ=255)
                                                                                                                                      atGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACAACAGAAACGACGGAAAAATTAAAAGCCGAGCTGGGTATCGTCTCTTACATCATTCATGGCTTAATGCGTATCGATActat              >  6:322907/1‑147 (MQ=255)
                                                                                                                                         gcgcAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACATCATTCATGGCTTAATGCGTATGGATACTCTGc            <  2:19819/148‑1 (MQ=255)
                                                                                                                                                 caacaaACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACATCATTCATGGCTTAATGCGTATGGATACTCTGCAACCGGAcc   >  8:35061/1‑149 (MQ=255)
                                                                                                                                                  aacaaacaAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCAGCGCTGGGTAGCGTCTCTTACATCATTCATGGCTTAATGCGTATGGATACTCTGCAACCGGACCg  >  7:142012/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGGGGGATGACATACCCGCGCTGGGAATTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGTAGGGATTCCTGAGGCACTGGATAAGGCGCTGAAACTGGCAATTGCCGGTGACGCCATTGCGATAGATATGGCGCAGGTCAACAAACAAACCTGTTTTATTAATATGGCGACAGGCGGATTTGGGACGCGTATTACCACAGAAACGCCGGAAAAATTAAAAGCCGCGCTGGGTAGCGTCTCTTACATCATTCATGGCTTAATGCGTATGGATACTCTGCAACCGGACCG  >  NC_000913/2168953‑2169247

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