Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 324,575 T→C 8.3% intergenic (+122/+121) ykgG → / ← ykgH LutC family protein; putative electron transport chain YkgEFG component/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913324,5750TC 8.3% 275.3 / 12.0 84intergenic (+122/+121)ykgG/ykgHLutC family protein; putative electron transport chain YkgEFG component/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (38/39);  new base C (4/3);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.58e-01

TATCTGATTATTGAGGATTGTTGAGAGAGATATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGCCCAGAGTTTTG  >  NC_000913/324430‑324713
                                                                                                                                                 |                                                                                                                                          
tATCTGATTATTGAGGATTGTTGAGAGAGATATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTAt                                                                                                                                          <  4:818355/148‑1 (MQ=255)
   cTGATTATTGAGGATTGTTGAGAGAGATATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTAc                                                                                                                                       >  4:267549/1‑148 (MQ=255)
      attattGAGGATTGTTGAGAGAGATATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGCAttggtaa                                                                                                                                    >  2:833062/1‑143 (MQ=255)
      attattGAGGATTGTTGAGAGAGATATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTATAAATATCCTGGTATTACCaa                                                                                                                                    >  4:826835/1‑148 (MQ=255)
         attGAGGATTGTTGAGAGAGATATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGc                                                                                                                                 <  3:267549/148‑1 (MQ=255)
          ttGAGGATTGTTGAGAGAGATATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGcc                                                                                                                                <  3:826835/148‑1 (MQ=255)
           tGAGGATTGTTGAGAGAGATATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCa                                                                                                                               >  3:962672/1‑148 (MQ=255)
                     tgAGAGAGATATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTAc                                                                                                                     <  4:962672/148‑1 (MQ=255)
                        gagagaTATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCaa                                                                                                                                    <  4:349481/130‑1 (MQ=255)
                        gagagaTATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCaa                                                                                                                                    >  3:349481/1‑130 (MQ=255)
                                   aGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTAc                                                                                                                                       <  3:1055180/116‑1 (MQ=255)
                                   aGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTAc                                                                                                                                       >  4:1055180/1‑116 (MQ=255)
                                          aCATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCaa                                                                                                                                    >  1:701122/1‑112 (MQ=255)
                                          aCATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCaa                                                                                                                                    <  2:701122/112‑1 (MQ=255)
                                          aCATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGcc                                                                                                                                >  4:24258/1‑116 (MQ=255)
                                          aCATAAACCTGAAGAAAGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGcc                                                                                                                                <  3:24258/116‑1 (MQ=255)
                                            aTAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCa                                                                                                                                     >  2:684013/1‑109 (MQ=255)
                                            aTAAAACTGAAGAACGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCa                                                                                                                                     <  1:684013/109‑1 (MQ=255)
                                                                      tGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCa                                                                                                    >  4:1120523/1‑116 (MQ=255)
                                                                      tGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCa                                                                                                    <  3:1120523/116‑1 (MQ=255)
                                                                        aaaaTAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGaaa                                                                                       >  2:1068090/1‑127 (MQ=255)
                                                                        aaaaTAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGaaa                                                                                       <  1:1068090/127‑1 (MQ=255)
                                                                         aaaTAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCaa                                                                                                                                    >  3:530363/1‑81 (MQ=255)
                                                                         aaaTAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCaa                                                                                                                                    <  4:530363/81‑1 (MQ=255)
                                                                           aTAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATAttt                                                                                                                       >  3:607372/1‑92 (MQ=255)
                                                                           aTAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATAttt                                                                                                                       <  4:607372/92‑1 (MQ=255)
                                                                               cGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACt                                                           >  3:1145307/1‑148 (MQ=255)
                                                                                  ttAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGAc                                                                                                         >  4:108984/1‑99 (MQ=255)
                                                                                  ttAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGAc                                                                                                         <  3:108984/99‑1 (MQ=255)
                                                                                  ttAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACa                                                                                                        <  4:1085918/100‑1 (MQ=255)
                                                                                  ttAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACa                                                                                                        >  3:1085918/1‑100 (MQ=255)
                                                                                       aCAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTAc                                                                                                                     <  3:1122985/82‑1 (MQ=255)
                                                                                       aCAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTAc                                                                                                                     >  4:1122985/1‑82 (MQ=255)
                                                                                            cTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACa                                                                                                                    >  1:164091/1‑78 (MQ=255)
                                                                                            cTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACa                                                                                                                    <  2:164091/78‑1 (MQ=255)
                                                                                             tGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCtt                                                    <  2:436606/141‑1 (MQ=255)
                                                                                             tGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCtt                                                    >  1:436606/1‑141 (MQ=255)
                                                                                             tGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGCATTGGTAATACCAGGata                                                                                                                          <  2:628077/71‑1 (MQ=21)
                                                                                             tGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGCATTGGTAATACCAGGata                                                                                                                          >  1:628077/1‑71 (MQ=21)
                                                                                               aCTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGc                                                      >  4:56939/1‑137 (MQ=255)
                                                                                               aCTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGc                                                      <  3:56939/137‑1 (MQ=255)
                                                                                                cTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACg                                                                                    >  1:304544/1‑106 (MQ=255)
                                                                                                cTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACg                                                                                    <  2:304544/106‑1 (MQ=255)
                                                                                                    aCCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCAtt                                                                                              >  1:283096/1‑92 (MQ=255)
                                                                                                    aCCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCAtt                                                                                              <  2:283096/92‑1 (MQ=255)
                                                                                                    aCCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCActct                                      >  3:167504/1‑148 (MQ=255)
                                                                                                      cccGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACg                                                                                    <  2:489909/100‑1 (MQ=255)
                                                                                                      cccGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACg                                                                                    >  1:489909/1‑100 (MQ=255)
                                                                                                       ccGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCAtt                                                                 >  1:351412/1‑118 (MQ=255)
                                                                                                       ccGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCAtt                                                                 <  2:351412/118‑1 (MQ=255)
                                                                                                        cGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACa                                              <  4:693349/136‑1 (MQ=255)
                                                                                                        cGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACa                                              >  3:693349/1‑136 (MQ=255)
                                                                                                                  tGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAg                                                                                                           <  4:692652/65‑1 (MQ=255)
                                                                                                                  tGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAg                                                                                                           >  3:692652/1‑65 (MQ=255)
                                                                                                                  tGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAAtgttg                        <  4:1145307/148‑1 (MQ=255)
                                                                                                                     ttATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGAc                                                                                                         >  2:488101/1‑64 (MQ=255)
                                                                                                                     ttATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGAc                                                                                                         <  1:488101/64‑1 (MQ=255)
                                                                                                                     ttATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCActctct                                    >  3:590462/1‑133 (MQ=255)
                                                                                                                     ttATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATACATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCActctct                                    <  4:590462/133‑1 (MQ=255)
                                                                                                                            gctCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGttcccttc                                                                          >  4:970160/1‑88 (MQ=255)
                                                                                                                            gctCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGttcccttc                                                                          <  3:970160/88‑1 (MQ=255)
                                                                                                                                   tGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGCCCag         >  1:429244/1‑146 (MQ=255)
                                                                                                                                   tGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGCCCag         <  2:429244/146‑1 (MQ=255)
                                                                                                                                    gTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATAc                                               >  1:21724/1‑107 (MQ=255)
                                                                                                                                    gTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATAc                                               <  2:21724/107‑1 (MQ=255)
                                                                                                                                    gTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCActct                                      >  1:249187/1‑116 (MQ=255)
                                                                                                                                    gTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCActct                                      <  2:249187/116‑1 (MQ=255)
                                                                                                                                    gTAAATATCCTGGCATTGGTAATACCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTc                                                              >  3:263398/1‑92 (MQ=25)
                                                                                                                                    gTAAATATCCTGGCATTGGTAATACCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTc                                                              >  3:443758/1‑92 (MQ=25)
                                                                                                                                    gTAAATATCCTGGCATTGGTAATACCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTc                                                              <  4:443758/92‑1 (MQ=25)
                                                                                                                                    gTAAATATCCTGGCATTGGTAATACCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTc                                                              <  4:263398/92‑1 (MQ=25)
                                                                                                                                     tAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGCCCAGAGtt     >  1:424476/1‑148 (MQ=255)
                                                                                                                                      aaaTATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTc                                                              <  2:28074/90‑1 (MQ=255)
                                                                                                                                      aaaTATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTc                                                              >  1:28074/1‑90 (MQ=255)
                                                                                                                                      aaaTATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGCCCAGAGttt    <  2:424476/148‑1 (MQ=255)
                                                                                                                                       aaTATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTg                                                                <  1:409314/87‑1 (MQ=255)
                                                                                                                                       aaTATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTg                                                                >  2:409314/1‑87 (MQ=255)
                                                                                                                                        atatCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTGCCTTCCCGGTCATTGTCACTTCAGCTTTATACAg                                             <  2:293312/105‑1 (MQ=255)
                                                                                                                                        atatCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAg                                             >  1:293312/1‑105 (MQ=255)
                                                                                                                                        atatCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGCCCAGAGTTttg  <  4:167504/148‑1 (MQ=255)
                                                                                                                                               ggTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGaa                             >  2:897948/1‑114 (MQ=255)
                                                                                                                                               ggTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGaa                             <  1:897948/114‑1 (MQ=255)
                                                                                                                                               ggTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTttat                   <  3:13073/124‑1 (MQ=255)
                                                                                                                                               ggTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTttat                   >  4:13073/1‑124 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
TATCTGATTATTGAGGATTGTTGAGAGAGATATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGCCCAGAGTTTTG  >  NC_000913/324430‑324713

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