Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,148,485 T→C 5.4% Y289H (TAT→CAT)  plsX → putative phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,148,4850TC 5.4% 231.9 / 3.3 74Y289H (TAT→CAT) plsXputative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (31/39);  new base C (2/2);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

ACGATGGAAGGTGTTGTCAGGATGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGC  >  NC_000913/1148341‑1148618
                                                                                                                                                |                                                                                                                                     
aCGATGGAAGGTGTTGTCAGGATGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGtat                                                                                                                                     >  8:51159/1‑147 (MQ=255)
     ggAAGGTGTTGTCAGGATGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGtata                                                                                                                                    <  7:40774/143‑1 (MQ=255)
     ggAAGGTGTTGTCAGGATGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGtata                                                                                                                                    >  8:40774/1‑143 (MQ=255)
      gAAGGTGTTGTCAGGATGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGgcg                                                                                                                              <  3:276395/148‑1 (MQ=255)
        aGGTGTTGTCAGGATGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGcc                                                                                                                            <  4:392984/148‑1 (MQ=255)
           tgttgtCAGGATGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCctgt                                                                                                                         <  5:280581/148‑1 (MQ=255)
              tgtCAGGATGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGtata                                                                                                                                    <  3:91688/134‑1 (MQ=255)
              tgtCAGGATGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGtata                                                                                                                                    >  4:91688/1‑134 (MQ=255)
                tCAGGATGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCctgtctgt                                                                                                                     >  8:126222/1‑147 (MQ=255)
                 cAGGATGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCACCCCCGACCAGTATAGCGGCGCctgtctgg                                                                                                                     >  2:112579/1‑145 (MQ=255)
                   ggATGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAACCGGGCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAgg                                                                                                                 <  4:232812/148‑1 (MQ=255)
                      tGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACAGCGCCTGTCTGTTAGGAtt                                                                                                              <  1:374483/148‑1 (MQ=255)
                      tGTTCCTTTCTCTGCTGAAATCTCAAGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGAtt                                                                                                              <  7:24012/148‑1 (MQ=255)
                          ccTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTgcgc                                                                                                          >  5:5645/1‑148 (MQ=255)
                             ttCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGtt                                                                                                                    >  7:141621/1‑135 (MQ=255)
                             ttCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGtt                                                                                                                    <  8:141621/135‑1 (MQ=255)
                             ttCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCgg                                                                                                        >  3:299582/1‑147 (MQ=255)
                              tctcTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGCTTCAGTCACCTCAACCCCGACCAGTATAAAGGCGCCTGTCTGTTAGGATTGCGCGGCa                                                                                                      >  7:139273/1‑148 (MQ=255)
                                 ctgctgAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATaa                                                                                                                                   >  2:300748/1‑116 (MQ=255)
                                 ctgctgAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATaa                                                                                                                                   <  1:300748/116‑1 (MQ=255)
                                  tgctgAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACg                                                                                                    <  1:45023/146‑1 (MQ=255)
                                    ctgAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGtata                                                                                                                                    <  4:354334/112‑1 (MQ=255)
                                    ctgAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGtata                                                                                                                                    >  3:354334/1‑112 (MQ=255)
                                    ctgAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGa                                                                                                >  5:100485/1‑148 (MQ=255)
                                      gAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCCGTCTCCTCAACCCCGCCCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGAt                                                                                               >  7:315230/1‑147 (MQ=255)
                                        aaTCTCAGGGTGAAGGGAAAAACCGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATaaa                                                                                            <  2:323521/148‑1 (MQ=255)
                                         aTCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATaaaa                                                                                           <  6:5645/148‑1 (MQ=255)
                                         aTCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATaaaa                                                                                           <  1:91110/148‑1 (MQ=255)
                                         aTCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATaaaa                                                                                           <  1:19221/148‑1 (MQ=255)
                                           ctcAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGAt                                                                                                               >  8:49132/1‑126 (MQ=255)
                                           ctcAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGAt                                                                                                               <  7:49132/126‑1 (MQ=255)
                                            tcAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAACGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGt                                                                                        >  5:218427/1‑148 (MQ=255)
                                             cAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCAc                                                                                                     <  5:336340/134‑1 (MQ=255)
                                             cAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCAc                                                                                                     >  6:336340/1‑134 (MQ=255)
                                              aGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATaaaaa                                                                                          <  4:299582/144‑1 (MQ=255)
                                               gggTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCAt                                                                                     <  7:51159/148‑1 (MQ=255)
                                                ggTGAAGGGAAAAACCGGTCTTGGTGGCTACTGTTATTAAAGCGGTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCAt                                                                                     <  8:315230/147‑1 (MQ=255)
                                                ggTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATg                                                                                    <  8:40124/148‑1 (MQ=255)
                                                       ggAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGc                                                                             >  4:144653/1‑148 (MQ=255)
                                                        gAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAACGAGCCTGACGAGGCGATTCAGTCCCCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGc                                                                             >  1:311523/1‑147 (MQ=255)
                                                        gAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTa                                                                                                                   <  6:341232/109‑1 (MQ=255)
                                                        gAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTa                                                                                                                   >  5:341232/1‑109 (MQ=255)
                                                         aaaaaaCGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCa                                                                           <  1:112579/148‑1 (MQ=255)
                                                           aaaaCGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAAt                                                                         >  3:171574/1‑148 (MQ=255)
                                                             atcGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCATTCACCTCAACCCCGACCAGTATAACGCCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCa                                                                       <  2:151983/146‑1 (MQ=255)
                                                             aaCGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGCCCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAAt                                                                         >  6:428933/1‑146 (MQ=255)
                                                                     gggtggCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCtt                                                               <  1:394434/147‑1 (MQ=255)
                                                                                gTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTcgc                                                     >  4:302274/1‑147 (MQ=255)
                                                                                gTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTcgc                                                     >  1:425342/1‑147 (MQ=255)
                                                                                 ttattaAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTcgcg                                                    <  4:171574/147‑1 (MQ=255)
                                                                                 ttattaAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTcgcg                                                    <  3:144653/147‑1 (MQ=255)
                                                                                 ttattaAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGa                                                   >  8:120533/1‑148 (MQ=255)
                                                                                  tattaAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTg                                                            >  5:340170/1‑138 (MQ=255)
                                                                                  tattaAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTg                                                            <  6:340170/138‑1 (MQ=255)
                                                                                  tattaAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGAt                                                  >  8:54218/1‑148 (MQ=255)
                                                                                     taAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATaaaa                                                                                           >  3:410793/1‑104 (MQ=255)
                                                                                      aaaGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGa                                               <  3:302274/147‑1 (MQ=255)
                                                                                        aGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGCATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCAt                                                                                     >  3:138143/1‑107 (MQ=255)
                                                                                        aGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGCATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCAt                                                                                     >  3:144159/1‑107 (MQ=255)
                                                                                        aGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGCATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCAt                                                                                     <  4:138143/107‑1 (MQ=255)
                                                                                        aGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGCATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCAt                                                                                     <  4:144159/107‑1 (MQ=255)
                                                                                           gTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCGTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGc                                         <  7:54218/148‑1 (MQ=255)
                                                                                               gCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGt                                      <  7:120533/147‑1 (MQ=255)
                                                                                                cTACAAAAGAGCCTGTCGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGc                                    <  6:218427/148‑1 (MQ=255)
                                                                                                 tACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCACCCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTg                                     >  1:179916/1‑146 (MQ=255)
                                                                                                     aaaGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTc                                                                                       >  5:130677/1‑92 (MQ=255)
                                                                                                     aaaGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTc                                                                                       <  6:130677/92‑1 (MQ=255)
                                                                                                                cGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTg                                                                                                 >  4:103433/1‑71 (MQ=255)
                                                                                                                cGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTg                                                                                                 <  3:103433/71‑1 (MQ=255)
                                                                                                                 gAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACa                    <  2:311523/147‑1 (MQ=255)
                                                                                                                  aGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAg                  >  2:233974/1‑148 (MQ=255)
                                                                                                                   ggCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGt                 <  1:233974/148‑1 (MQ=255)
                                                                                                                            gTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGc                                                                             >  5:116836/1‑79 (MQ=255)
                                                                                                                            gTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGc                                                                             <  6:116836/79‑1 (MQ=255)
                                                                                                                            gTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCg        <  6:100485/148‑1 (MQ=255)
                                                                                                                             tCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGa       <  2:425342/148‑1 (MQ=255)
                                                                                                                                  tCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGc  >  6:418610/1‑148 (MQ=255)
                                                                                                                                         cGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGa                                               >  2:407986/1‑96 (MQ=255)
                                                                                                                                         cGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGa                                               <  1:407986/96‑1 (MQ=255)
                                                                                                                                           aCCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCa                                        <  5:376578/101‑1 (MQ=255)
                                                                                                                                           aCCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCa                                        >  6:376578/1‑101 (MQ=255)
                                                                                                                                                |                                                                                                                                     
ACGATGGAAGGTGTTGTCAGGATGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGTGGTGGCTACTGTTATTAAAGCGTTGGCTACAAAAGAGCCTGACGAGGCGATTCAGTCACCTCAACCCCGACCAGTATAACGGCGCCTGTCTGTTAGGATTGCGCGGCACGGTGATAAAAAGTCATGGTGCAGCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGC  >  NC_000913/1148341‑1148618

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