Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 650,244 C→A 6.2% E83* (GAA→TAA)  citE ← citrate lyase, citryl‑ACP lyase (beta) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,2440CA 6.2% 337.5 / 5.5 96E83* (GAA→TAA) citEcitrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base C (37/53);  new base A (3/3);  total (40/56)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCAT  >  NC_000913/650106‑650381
                                                                                                                                          |                                                                                                                                         
gTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCg                                                                                                                                          <  1:1611420/140‑1 (MQ=255)
 ttCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCgg                                                                                                                                         >  1:1624343/1‑140 (MQ=255)
 ttCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTAgg                                                                                                                                         >  2:334101/1‑140 (MQ=255)
 ttCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTAgg                                                                                                                                         >  2:334114/1‑140 (MQ=255)
 ttCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTAgg                                                                                                                                         >  2:334133/1‑140 (MQ=255)
 ttCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGATAACACCCCAGTCgg                                                                                                                                         >  1:125421/1‑140 (MQ=255)
           caGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAgcgc                                                                                                                               >  2:1138672/1‑140 (MQ=255)
           caGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAgcgc                                                                                                                               >  2:1195101/1‑140 (MQ=255)
           caGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAgcgc                                                                                                                               >  2:320605/1‑140 (MQ=255)
            aGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGcgcg                                                                                                                              <  1:401732/140‑1 (MQ=255)
            aGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGcgcg                                                                                                                              <  1:1627676/140‑1 (MQ=255)
                 ttttCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGac                                                                                                                         <  2:406677/140‑1 (MQ=255)
                 ttttCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGac                                                                                                                         <  1:1443599/140‑1 (MQ=255)
                 ttttCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGac                                                                                                                         <  2:406505/140‑1 (MQ=255)
                  tttCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGaca                                                                                                                        >  1:1074127/1‑140 (MQ=255)
                        tACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGcaca                                                                                                                  <  2:559068/140‑1 (MQ=255)
                        tACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGcaca                                                                                                                  <  2:1136086/140‑1 (MQ=255)
                        tACGCAGGATCTCTTATTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGcaca                                                                                                                  >  1:1019663/1‑140 (MQ=255)
                             aGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGt                                                                                                             <  1:1499538/140‑1 (MQ=255)
                             aGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGt                                                                                                             <  1:315428/140‑1 (MQ=255)
                             aGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGt                                                                                                             <  1:985518/140‑1 (MQ=255)
                                    ctTTTTCACTATCCCGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACTCCCCTTTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAata                                                                                                      <  2:1185487/140‑1 (MQ=255)
                                    ctTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCTGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAata                                                                                                      <  2:1209823/140‑1 (MQ=255)
                                       tttCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATcg                                                                                                   <  2:409095/140‑1 (MQ=255)
                                       tttCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATcg                                                                                                   >  1:675068/1‑140 (MQ=255)
                                       tttCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATcg                                                                                                   >  1:11911/1‑140 (MQ=255)
                                         tCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATcgcg                                                                                                 >  2:638184/1‑140 (MQ=255)
                                             tatCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACACCGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATAc                                                                                             >  2:1595827/1‑140 (MQ=255)
                                             tatCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCaa                                                                                                        >  2:357143/1‑129 (MQ=255)
                                             tatCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCaa                                                                                                        <  1:357143/129‑1 (MQ=255)
                                              atCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACa                                                                                            <  1:1025444/140‑1 (MQ=255)
                                              atCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACa                                                                                            <  2:1130798/140‑1 (MQ=255)
                                                ccAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGc                                                                                          <  2:630582/140‑1 (MQ=255)
                                                ccAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGc                                                                                          <  2:1334673/140‑1 (MQ=255)
                                                ccAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGc                                                                                          <  2:1336910/140‑1 (MQ=255)
                                                    taCATCCTTAGCGGTATCGGATTTCGGAAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTAAACACCCAATTCGGATTCCAGCGCGTTGACACGCAAAATGGGATCAATGTCGCGAGACAGCGGAt                                                                                      >  2:642473/2‑140 (MQ=255)
                                                         ccTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGc                                                                                 >  2:594586/1‑140 (MQ=255)
                                                         ccTGAGCGGTAACGGTTTTCGTCAGACGAACAAAGTCCGCACCACCGCGAACTACGGCTTCCAGGTAGTTAACACCCCATTCGGAATACAGCTCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGc                                                                                 >  2:594706/1‑140 (MQ=255)
                                                            gAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAgc                                                                              <  2:476295/140‑1 (MQ=255)
                                                             aGCGGTATCGGTTTTCGGCAGACGCACACCGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGACTCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAgcg                                                                             <  1:1645412/140‑1 (MQ=255)
                                                             aGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAgcg                                                                             <  1:1311870/140‑1 (MQ=255)
                                                              gCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCATCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAgcgc                                                                            >  2:937529/1‑140 (MQ=255)
                                                              gCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAgcgc                                                                            >  2:1215144/1‑140 (MQ=255)
                                                                   aTCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGt                                                                       >  2:660549/1‑140 (MQ=255)
                                                                   aTCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGt                                                                       <  1:149048/140‑1 (MQ=255)
                                                                       gTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTATCACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTGAAc                                                                   >  1:1384509/1‑140 (MQ=255)
                                                                       gTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGAAACAGCGGATGTTGCAGCGCGTGGTAAAc                                                                   >  1:916820/1‑140 (MQ=255)
                                                                          ttCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCAt                                                                >  1:1469754/1‑140 (MQ=255)
                                                                            cGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCATTATCGCGTTACAGCGGATGTTGCAGCGCGTGGTAAACCATgc                                                              >  1:1380922/1‑140 (MQ=255)
                                                                               cAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATgcggc                                                           <  1:1227667/140‑1 (MQ=255)
                                                                                 gACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGcggcgg                                                         >  2:1630819/1‑140 (MQ=255)
                                                                                  aCGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCggcggg                                                        >  1:1436216/1‑140 (MQ=255)
                                                                                                caccacCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCa                                          >  2:474733/1‑140 (MQ=255)
                                                                                                  ccaccTCGAACGACGGCTTCCATGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCCCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACg                                        <  2:1654951/140‑1 (MQ=255)
                                                                                                  ccaccGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACg                                        <  1:1088745/140‑1 (MQ=255)
                                                                                                  ccaccGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACg                                        <  1:1452788/140‑1 (MQ=255)
                                                                                                  ccaccGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACg                                        <  1:1452849/140‑1 (MQ=255)
                                                                                                  ccaccGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACg                                        <  1:834087/140‑1 (MQ=255)
                                                                                                    accGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCa                                      <  2:518287/140‑1 (MQ=255)
                                                                                                    accGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCa                                      <  2:518288/140‑1 (MQ=255)
                                                                                                        cgAACGACGGCTTCCAGGTCGTTAACACCCCGTTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGATGCAGCGCGAGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGc                                  >  2:466121/1‑140 (MQ=255)
                                                                                                        cgAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGc                                  >  2:1007638/1‑140 (MQ=255)
                                                                                                        cgAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGc                                  >  1:1429594/1‑140 (MQ=255)
                                                                                                        cgAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGACTTTTTCACGCacagc                                  >  1:971503/1‑136 (MQ=255)
                                                                                                          aaCGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTa                                <  1:1168956/140‑1 (MQ=255)
                                                                                                            cgacgGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACg                              <  2:148867/140‑1 (MQ=255)
                                                                                                             gacgGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACgg                             <  1:393923/140‑1 (MQ=255)
                                                                                                             gacgGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACgg                             <  1:1234668/140‑1 (MQ=255)
                                                                                                              acgGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGa                            <  2:1574974/140‑1 (MQ=255)
                                                                                                              acgGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGa                            <  1:736488/140‑1 (MQ=255)
                                                                                                                    tCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTc                      <  1:976473/140‑1 (MQ=255)
                                                                                                                    tCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTc                      <  1:894951/140‑1 (MQ=255)
                                                                                                                    tCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTc                      <  1:976503/140‑1 (MQ=255)
                                                                                                                    tCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTc                      <  1:894674/140‑1 (MQ=255)
                                                                                                                    tCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTc                      <  1:894939/140‑1 (MQ=255)
                                                                                                                           cGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCa               <  1:721588/140‑1 (MQ=255)
                                                                                                                           cGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCa               <  1:764118/140‑1 (MQ=255)
                                                                                                                           cGTTAACACCCCATTAGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCa               <  1:334101/140‑1 (MQ=255)
                                                                                                                           cGTTAACACCCCATTAGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCa               <  1:334114/140‑1 (MQ=255)
                                                                                                                           cGTTAACACCCCATTAGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCa               <  1:334133/140‑1 (MQ=255)
                                                                                                                             ttAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCaaa             <  2:842522/140‑1 (MQ=255)
                                                                                                                             ttAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCaaa             <  2:466937/140‑1 (MQ=255)
                                                                                                                                 cacCCCATGCGGAATCCAGCGCGTTGACACGCACAAGGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATc         <  2:74853/140‑1 (MQ=255)
                                                                                                                                   ccccATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAg       >  2:1506926/1‑140 (MQ=255)
                                                                                                                                   ccccATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGAAGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTGCGAGGTCAAACATCAg       >  2:1507007/1‑140 (MQ=255)
                                                                                                                                   ccccATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGAGGATGTTGCAGCGCGTGGTACACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGTGGTCAAACATctg       >  2:1507078/1‑138 (MQ=255)
                                                                                                                                      cATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGc    <  1:921743/140‑1 (MQ=255)
                                                                                                                                      cATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGc    >  1:1116736/1‑140 (MQ=255)
                                                                                                                                       aTTCTGAATCCAGAGCGTTGACACGCACAATGGTTTCAATATCTAGATACAGCGGATGATGCAGAGCGAGGTAAACCATGCTGCGGGCGGTGACTATTACACGAAATGCAACGGAGTCTTCTAGGTCAAACAACATGGCa   >  1:1394000/1‑140 (MQ=255)
                                                                                                                                       aTTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATCTCGCGATACAGCGGAGGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTACACGCACTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCa   >  2:1029006/1‑140 (MQ=255)
                                                                                                                                       aTTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCa   >  2:1029059/1‑140 (MQ=255)
                                                                                                                                       aTTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGACAAACATCAGGGCa   >  2:192279/1‑140 (MQ=255)
                                                                                                                                        ttCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCAt  <  1:941237/140‑1 (MQ=255)
                                                                                                                                        ttCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCAt  <  1:630456/140‑1 (MQ=255)
                                                                                                                                        ttCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCAt  <  1:664075/140‑1 (MQ=255)
                                                                                                                                        ctCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCAt  <  1:574924/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
GTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCAT  >  NC_000913/650106‑650381

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