Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,076,349 T→A 5.9% D845V (GAT→GTT)  putA ← fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,076,3490TA 5.9% 223.4 / 4.2 68D845V (GAT→GTT) putAfused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (29/35);  new base A (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

CGACAATCATCGCGTTCATGCCGCCGGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGTTTGTTCTGCCGGTTTTGCCAGCACGCTGTTAC  >  NC_000913/1076206‑1076496
                                                                                                                                               |                                                                                                                                                   
cGACAATCATCGCGTTCATGCCGCCGGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCAACAcc                                                                                                                                                 >  5:140317/1‑148 (MQ=255)
cGACAATCATCGCGTTCATGCCGCCGGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCAACAcc                                                                                                                                                 >  5:145058/1‑148 (MQ=255)
cGACAATCATCGCGTTCATGCCGCCGGTTTCAGCGATGAGCGGAATAGGGCGACCCAGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTACCGACTTCGGTTGAACCGGTACACATAACCCCGCGCACGCGATCATCAcc                                                                                                                                                 >  8:263367/1‑148 (MQ=255)
cGACAATCATCGCGTTCATGCCGCCGGTGTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGGTGAACCGGTAAACATCACCCCGCGCACGCGATCATCAcc                                                                                                                                                 >  6:121162/1‑148 (MQ=255)
 gACAATCATCGCGTTCATGCCGCCGGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCAcc                                                                                                                                                 <  6:147308/147‑1 (MQ=255)
   cAATCATCGCGTTCATGCCGCCGGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGt                                                                                                                                              >  2:96022/1‑148 (MQ=255)
      tcatcGCGTTCATGCCGCCGGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAg                                                                                                                                           >  7:125694/1‑148 (MQ=255)
         tcGCGTTCATGCCGCCGGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTg                                                                                                                                        <  7:317471/148‑1 (MQ=255)
          cgcgTTCATGCCGCCGGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTgc                                                                                                                                       >  4:232043/1‑148 (MQ=255)
              ttCATGCCGCCGGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTgcgcgc                                                                                                                                   >  5:241831/1‑148 (MQ=255)
                 aTGCCGCCGGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCa                                                                                                                                <  8:125694/148‑1 (MQ=255)
                 aTGCCGCCGGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCa                                                                                                                                <  3:86036/148‑1 (MQ=255)
                    ccgccgGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACgg                                                                                                                             >  7:4206/1‑148 (MQ=255)
                     cgccgGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGt                                                                                                                            <  4:48596/148‑1 (MQ=255)
                      gccgGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGtt                                                                                                                           >  1:59017/1‑148 (MQ=255)
                       ccgGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGttt                                                                                                                          <  8:4206/148‑1 (MQ=255)
                        cgGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTc                                                                                                                         <  4:280262/148‑1 (MQ=255)
                         ggTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCa                                                                                                                        >  2:167958/1‑148 (MQ=255)
                           tttCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCAcc                                                                                                                      >  7:126371/1‑148 (MQ=255)
                              cAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCg                                                                                                                   >  8:77388/1‑148 (MQ=255)
                               aGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCCCCCGTCAGTTGCGCGCCCACGGTTTCACCCCGa                                                                                                                  >  4:105497/1‑148 (MQ=255)
                               aGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGa                                                                                                                  <  1:286745/148‑1 (MQ=255)
                                  gATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCt                                                                                                               >  7:314450/1‑148 (MQ=255)
                                  gATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCt                                                                                                               >  2:26204/1‑148 (MQ=255)
                                   aTGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTg                                                                                                              <  1:57718/148‑1 (MQ=255)
                                    tGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTgg                                                                                                             <  1:26204/148‑1 (MQ=255)
                                         ggAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGgcagca                                                                                                        <  8:126371/148‑1 (MQ=255)
                                          gAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCaa                                                                                                       <  6:268674/148‑1 (MQ=255)
                                            aTAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAAtt                                                                                                     <  2:221515/148‑1 (MQ=255)
                                             tAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTg                                                                                                    <  1:96022/148‑1 (MQ=255)
                                                ggCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGtt                                                                                                                           >  1:222768/1‑122 (MQ=255)
                                                ggCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGtt                                                                                                                           <  2:222768/122‑1 (MQ=255)
                                                    aCCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACg                                                                                             >  8:112590/1‑148 (MQ=255)
                                                         gAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTgg                                                                                        >  8:96884/1‑148 (MQ=255)
                                                         gAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTgg                                                                                        >  4:37479/1‑148 (MQ=255)
                                                          aGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTg                                                                                         <  6:268402/146‑1 (MQ=255)
                                                           gCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCg                                                                                      <  7:77388/148‑1 (MQ=255)
                                                            cGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCgg                                                                                     <  2:59017/148‑1 (MQ=255)
                                                               ccAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTGCGGTTGAACCGGTAACCATCACCCCGCGCACGCGAGCATCACCCGTCAGTTGCGCGCCCACGGGGTCACCCCGACCTGGAAGCAATTGCACCACGCCTGGCGGTac                                                                                  >  5:96573/1‑148 (MQ=255)
                                                                 aGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTacac                                                                                <  6:241831/148‑1 (MQ=255)
                                                                  ggcggcTGGCGATATTGCGCTGCATTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTgcg                                                                                                                                      <  2:321642/93‑1 (MQ=255)
                                                                  ggcggcTGGCGATATTGCGCTGCACTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTgcg                                                                                                                                      >  1:321642/1‑93 (MQ=255)
                                                                    cggcTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTg                                                                                                    >  8:271095/1‑125 (MQ=255)
                                                                    cggcTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTg                                                                                                    <  7:271095/125‑1 (MQ=255)
                                                                       cTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCAACACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCtt                                                                          <  6:140317/148‑1 (MQ=255)
                                                                       cTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCAACACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCtt                                                                          <  6:145058/148‑1 (MQ=255)
                                                                           cGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAg                                                                      >  5:28945/1‑148 (MQ=255)
                                                                            gATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAg                                                                      >  5:124239/1‑147 (MQ=255)
                                                                                  gcgcTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATg                                                               >  7:270207/1‑148 (MQ=255)
                                                                                   cgcTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATgg                                                              <  7:96884/148‑1 (MQ=255)
                                                                                         aGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATcc                                                        <  1:167958/148‑1 (MQ=255)
                                                                                              cGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTgc                                                   <  6:28945/148‑1 (MQ=255)
                                                                                               gTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCAc                                                                                                                                                  <  3:176756/52‑1 (MQ=255)
                                                                                               gTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCAc                                                                                                                                                  >  4:176756/1‑52 (MQ=255)
                                                                                               gTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTgcg                                                  <  8:314450/148‑1 (MQ=255)
                                                                                                    gACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCa                                             <  7:112590/148‑1 (MQ=255)
                                                                                                          ggatGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAgc                                       <  5:121162/145‑1 (MQ=255)
                                                                                                                  cGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGt                                  <  2:27246/145‑1 (MQ=255)
                                                                                                                  cGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGt                                  >  1:27246/1‑145 (MQ=255)
                                                                                                                    gTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGTTTGtt                             <  3:232043/148‑1 (MQ=255)
                                                                                                                     tAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACCCCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGTTTGTTc                            >  7:315693/1‑148 (MQ=255)
                                                                                                                           tCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGTTTGTTCTGCCgg                      >  6:243797/1‑148 (MQ=255)
                                                                                                                             aCCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGTTTGTTCTGCCGGtt                    <  3:37479/148‑1 (MQ=255)
                                                                                                                             aCCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGTTTGTTCTGCCGGt                     <  6:280123/147‑1 (MQ=255)
                                                                                                                              ccccGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGTTTGTTCTGCCGGttt                   <  5:243797/148‑1 (MQ=255)
                                                                                                                              ccccGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGTTTGTTCTGCCGGttt                   <  6:124239/148‑1 (MQ=255)
                                                                                                                                     cACGCGATAATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCCACCTGGCACCAATTGCACCACGCCTGGCGGTACACCCGCGACCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGTTTGTTCTGCCGCTTTTGCCAg             >  3:197227/1‑147 (MQ=255)
                                                                                                                                       cgcgATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGTTTGTTCTGCCGGTTTTGCCAGCAc          <  8:270207/148‑1 (MQ=255)
                                                                                                                                              atcaCCCGTCAGTTGCGCGCCGACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGTTTGTTCTGCCGGTTTTGCCAGCACGCTGTTa   <  7:263367/148‑1 (MQ=255)
                                                                                                                                               tcaCCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGTTTGTTCTGCCGGTTTTGCCAGCACGCTGTTAc  >  4:38522/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CGACAATCATCGCGTTCATGCCGCCGGTTTCAGCGATGAGCGGAATAGGGCGACCCTGAGCGTCCAGGCGGCTGGCGATATTGCGCTGCAGTAACGTAGCGACTTCGGTTGAACCGGTAAACATCACCCCGCGCACGCGATCATCACCCGTCAGTTGCGCGCCCACGGTTTCACCCCGACCTGGCAGCAATTGCACCACGCCTGGCGGTACACCCGCTTCCAGCAAAATGGCGATCCCTTGCGCGGCAATCAGCGGCGTTTGTTCTGCCGGTTTTGCCAGCACGCTGTTAC  >  NC_000913/1076206‑1076496

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: ‑