Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 642,936 G→A 9.3% L283L (CTG→CTA ybdR → putative Zn‑dependent NAD(P)‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913642,9360GA 9.3% 261.5 / 14.4 86L283L (CTG→CTAybdRputative Zn‑dependent NAD(P)‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (44/34);  new base A (4/4);  total (48/38)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

CCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAGGGTTGTCGTTTAAAATGG  >  NC_000913/642790‑643081
                                                                                                                                                  |                                                                                                                                                 
ccGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAATGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGa                                                                                                                                                 <  3:522274/149‑1 (MQ=255)
ccGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGa                                                                                                                                                 >  4:236751/1‑149 (MQ=255)
ccGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGa                                                                                                                                                 >  2:655689/1‑149 (MQ=255)
  gATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAgg                                                                                                                                               >  1:598993/1‑149 (MQ=255)
   aTTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAggg                                                                                                                                              >  1:450112/1‑149 (MQ=255)
    ttAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGgc                                                                                                                                             <  1:153964/149‑1 (MQ=255)
     tAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGgca                                                                                                                                            <  1:344223/149‑1 (MQ=255)
        ttttGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGgcagca                                                                                                                                         >  1:163407/1‑149 (MQ=255)
            gatgaAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCgg                                                                                                                                     <  1:655689/149‑1 (MQ=255)
            gatgaAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCgg                                                                                                                                     >  2:260798/1‑149 (MQ=255)
            gatgaAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCgg                                                                                                                                     <  3:236751/149‑1 (MQ=255)
             atgaAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGt                                                                                                                                    <  3:33638/149‑1 (MQ=255)
                aaGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTaaa                                                                                                                                 >  4:612700/1‑149 (MQ=255)
                aaGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTAGAGGGCAGCAGCGGTaaa                                                                                                                                 >  2:620583/1‑149 (MQ=255)
                aaGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTAGAGGGCAGCAGCGGTaaa                                                                                                                                 >  2:563428/1‑149 (MQ=255)
                aaGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTAGAGGGCAGCAGCGGTaaa                                                                                                                                 >  2:40024/1‑149 (MQ=255)
                aaGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTAGAGGGCAGCAGCGGTaaa                                                                                                                                 >  4:416342/1‑149 (MQ=255)
                  gACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGc                                                                                                                               >  3:61494/1‑149 (MQ=255)
                     aGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAACGGTTCCTGCCTACCCGGAAACGGGAGGGCACCGCGGGAAAGccgtg                                                                                                                           >  4:602227/1‑145 (MQ=255)
                     aGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGtt                                                                                                                            <  3:45253/149‑1 (MQ=255)
                     aGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGtt                                                                                                                            <  1:625770/149‑1 (MQ=255)
                     aGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTAGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGAGCAGCAGCGGTAAAGCGtt                                                                                                                            <  3:637340/149‑1 (MQ=255)
                                  aGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTAtt                                                                                                               <  1:260798/149‑1 (MQ=255)
                                  aGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTAtt                                                                                                               <  1:241558/149‑1 (MQ=255)
                                    tCAATTATTGAACAAACGGGAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTAGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGa                                                                                                                                                 <  3:141530/113‑1 (MQ=255)
                                    tCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGa                                                                                                                                                 >  1:36928/1‑113 (MQ=255)
                                    tCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGa                                                                                                                                                 >  4:141530/1‑113 (MQ=255)
                                    tCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGa                                                                                                                                                 <  2:36928/113‑1 (MQ=255)
                                    tCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGt                                                                                                                             >  1:457360/1‑133 (MQ=255)
                                    tCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGt                                                                                                                             <  2:457360/133‑1 (MQ=255)
                                    tCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTgc                                                                                                             >  2:804734/1‑149 (MQ=255)
                                    tCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTgc                                                                                                             >  2:220400/1‑149 (MQ=255)
                                        ttattGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTg                                                                                                                                                   <  3:296948/107‑1 (MQ=255)
                                        ttattGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTg                                                                                                                                                   >  4:296948/1‑107 (MQ=255)
                                        ttattGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAgt                                                                                                                    >  3:683442/1‑138 (MQ=255)
                                        ttattGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAgt                                                                                                                    <  4:683442/138‑1 (MQ=255)
                                          attGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGt                                                                                                       <  3:20448/149‑1 (MQ=255)
                                          attGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGt                                                                                                       <  1:775920/149‑1 (MQ=255)
                                              aaCAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAgg                                                                                                   >  3:577219/1‑149 (MQ=255)
                                              aaCAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAgg                                                                                                   >  3:375258/1‑149 (MQ=255)
                                              aaCAAACGGCAGGTCACCGGGGCGAGGATGCAGTAATAGACGCTGTCGGTTTTGAAGCGAAAGGCAGCACCCCGGAAACGGTGCGGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAgg                                                                                                   >  3:145145/1‑149 (MQ=255)
                                                cAAACGGCAGGTCACGGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCg                                                                                                 >  1:591519/1‑149 (MQ=255)
                                                 aaaCGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGt                                                                                                >  4:341640/1‑149 (MQ=255)
                                                        aGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCa                                                                                         >  1:627466/1‑149 (MQ=255)
                                                         ggTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCAt                                                                                        <  2:450112/149‑1 (MQ=255)
                                                          gTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCAtt                                                                                       >  4:263971/1‑149 (MQ=255)
                                                             aCCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTAGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGtt                                                                                    <  1:40024/149‑1 (MQ=255)
                                                             aCCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTAGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGtt                                                                                    <  1:563428/149‑1 (MQ=255)
                                                             aCCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTAGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGtt                                                                                    <  1:620583/149‑1 (MQ=255)
                                                             aCCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTAGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGtt                                                                                    <  3:416342/149‑1 (MQ=255)
                                                                 gggCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGtt                                                                                    >  3:6686/1‑145 (MQ=255)
                                                                 gggCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGtt                                                                                    <  4:6686/145‑1 (MQ=255)
                                                                 gggCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGc                                                                                 <  2:89628/148‑1 (MQ=255)
                                                                 gggCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGc                                                                                 >  1:89628/1‑148 (MQ=255)
                                                                             aGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGTAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCt                                                                    <  4:577219/149‑1 (MQ=255)
                                                                             aGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCt                                                                    >  2:387071/1‑149 (MQ=255)
                                                                              gTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTa                                                                   >  2:86641/1‑149 (MQ=255)
                                                                                    gACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTgg                                                             >  2:806498/1‑149 (MQ=255)
                                                                                     aCGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGa                                                            >  1:395470/1‑149 (MQ=255)
                                                                                     aCGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGa                                                            >  3:646427/1‑149 (MQ=255)
                                                                                        ccGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAATATGGAGGGCAGCAGCAGTAAAGCGTTGCGTCAGTGTATGGCGGCGGTCAGGCGTGGCTGCATTGTTAGCGTACCGGGCGTCTACGCTGGAttt                                                         <  4:145145/149‑1 (MQ=255)
                                                                                        ccGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGAttt                                                         <  3:263971/149‑1 (MQ=255)
                                                                                        ccGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGAttt                                                         <  4:375258/149‑1 (MQ=255)
                                                                                          gTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTAt                                                       <  3:97092/149‑1 (MQ=255)
                                                                                                tttGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACgg                                                 <  1:806498/149‑1 (MQ=255)
                                                                                                     aGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTcc                                            <  1:804734/149‑1 (MQ=255)
                                                                                                      gCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCt                                           >  1:691274/1‑149 (MQ=255)
                                                                                                             gcagcaCCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTAtt                                                                                                               <  3:226819/74‑1 (MQ=255)
                                                                                                             gcagcaCCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTAtt                                                                                                               >  4:226819/1‑74 (MQ=255)
                                                                                                                 caccacGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACg                                >  4:704083/1‑149 (MQ=255)
                                                                                                                    cacGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCt                             >  4:76217/1‑149 (MQ=255)
                                                                                                                       ggAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGttt                                       >  4:634817/1‑136 (MQ=255)
                                                                                                                       ggAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGttt                                       <  3:634817/136‑1 (MQ=255)
                                                                                                                       ggAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTg                          >  1:715067/1‑149 (MQ=255)
                                                                                                                       ggAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTg                          >  1:280965/1‑149 (MQ=255)
                                                                                                                            cGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATaaa                     <  2:627466/149‑1 (MQ=255)
                                                                                                                            cGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATaaa                     <  1:220400/149‑1 (MQ=255)
                                                                                                                                  tGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAGGGTTg               >  4:343181/1‑149 (MQ=255)
                                                                                                                                   gACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAGGGTTGt              >  1:644848/1‑149 (MQ=255)
                                                                                                                                   gACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAGGGTTGt              >  1:187006/1‑149 (MQ=255)
                                                                                                                                        aCCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCttt                           <  2:29697/131‑1 (MQ=255)
                                                                                                                                        aCCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCttt                           >  1:29697/1‑131 (MQ=255)
                                                                                                                                         ccTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAGGGTTGTCGTTTa        >  2:425802/1‑149 (MQ=255)
                                                                                                                                          cTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAGGGTTGTCGTTTaa       >  1:626411/1‑149 (MQ=255)
                                                                                                                                           tGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAGGGTTGTCGTTTaaa      >  4:701501/1‑149 (MQ=255)
                                                                                                                                               aCTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTTAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATaaa                     <  3:520338/130‑1 (MQ=255)
                                                                                                                                               aCTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTTAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTATTTGCCGACGCCTTTGATaaa                     >  4:520338/1‑130 (MQ=255)
                                                                                                                                               aCTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAGGGTTGTCGTTTAAAATgg  <  1:607927/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCGATTAATTTTGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAGGGTTGTCGTTTAAAATGG  >  NC_000913/642790‑643081

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