Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,460,624 A→G 6.4% intergenic (‑155/+26) fadI ← / ← yfcZ 3‑ketoacyl‑CoA thiolase FadI/DUF406 domain‑containing protein YfcZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,460,6240AG 6.4% 514.8 / 17.2 172intergenic (‑155/+26)fadI/yfcZ3‑ketoacyl‑CoA thiolase FadI/DUF406 domain‑containing protein YfcZ
Reads supporting (aligned to +/- strand):  ref base A (83/78);  new base G (8/3);  total (91/81)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

Only 100 of 172 total aligned reads displayed.

CCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTC  >  NC_000913/2460476‑2460772
                                                                                                                                                    |                                                                                                                                                    
ccccTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGAc                                                                                                                                                     <  2:1839775/150‑1 (MQ=255)
ccccTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGAc                                                                                                                                                     >  2:1917533/1‑150 (MQ=255)
     ttAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTATATTTATATTATTGGTGAAAATGGAATTTGTTATCTATTATACAGAGTTTAAATAAGTAGTGTAAAGATATCTCCCCTGCCAGCAACTGACAgggg                                                                                                                                                <  2:1692712/149‑1 (MQ=255)
       aaaaaaaTAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcg                                                                                                                                              <  2:423730/150‑1 (MQ=255)
        aaaaaaTAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgt                                                                                                                                             >  1:388367/1‑150 (MQ=255)
         aaaaaTAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGGGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgta                                                                                                                                            <  1:2094493/150‑1 (MQ=255)
          aaaaTAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgt                                                                                                                                             >  1:1399615/1‑148 (MQ=255)
            aataaataAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgtacgt                                                                                                                                         <  2:1861816/150‑1 (MQ=255)
             ataaataaAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgtacgta                                                                                                                                        >  1:333149/1‑150 (MQ=255)
                aataaaGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGAc                                                                                                                                                     >  1:498373/1‑134 (MQ=255)
                     aGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGagag                                                                                                                                <  1:436492/150‑1 (MQ=255)
                     aGTGGTCCGACCTGATCATAGTCTTAACCATTATATAACATTTACCCAAGGGGTGAAAAGGGAACGTGGTCGCTATGCCACAGAGAGAAAAGAAGCAGAGAAAAGTAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGagag                                                                                                                                <  2:2482324/150‑1 (MQ=255)
                            tGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTaa                                                                                                                         >  2:2243831/1‑150 (MQ=255)
                             tcccTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGACAAAAGGGAAAGTGGGACCCATGACACCGTGCCAAAACCACCAGCCAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAAc                                                                                                                        <  1:2398784/148‑1 (MQ=255)
                              aCCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACg                                                                                                                       <  1:1029055/150‑1 (MQ=255)
                               ccTGATCATAGTCTTAACCTTTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGc                                                                                                                      <  1:2058850/150‑1 (MQ=255)
                               ccTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGc                                                                                                                      >  2:590041/1‑150 (MQ=255)
                               ccTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGc                                                                                                                      >  2:2298884/1‑150 (MQ=255)
                                 tGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAg                                                                                                                    >  1:602267/1‑150 (MQ=255)
                                   aTCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGAc                                                                                                                  <  2:440392/150‑1 (MQ=255)
                                    tCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGAcc                                                                                                                 >  1:113391/1‑150 (MQ=255)
                                    tCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGAcc                                                                                                                 >  1:211698/1‑150 (MQ=255)
                                      aTAGTCTTAACCATTTTTTTACATTTAGCCAAGTGTAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCa                                                                                                               <  2:1532370/150‑1 (MQ=255)
                                              aaCCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGaaa                                                                                                       >  2:1662055/1‑150 (MQ=255)
                                              aaCCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGaaa                                                                                                       >  1:6948/1‑150 (MQ=255)
                                               aCCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGaaaa                                                                                                      <  2:645806/150‑1 (MQ=255)
                                                ccATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAACAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAAt                                                                                                     >  1:1944962/1‑150 (MQ=255)
                                                  aTTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCa                                                                                                   <  2:2390125/150‑1 (MQ=255)
                                                   tttttttACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAg                                                                                                  <  1:1917533/150‑1 (MQ=255)
                                                    ttttttACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGc                                                                                                 <  2:1080034/150‑1 (MQ=255)
                                                    ttttttACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGc                                                                                                 >  1:2355297/1‑150 (MQ=255)
                                                     tttttACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCa                                                                                                <  2:1246424/150‑1 (MQ=255)
                                                       tttACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGCCACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTCGAAAATCAGCAtt                                                                                              >  1:1730127/1‑150 (MQ=255)
                                                        ttACATTTAGCCAAGTGGAGAAAAGGGAAAGTGTGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTTCGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCAttt                                                                                             >  2:1868992/1‑150 (MQ=255)
                                                        ttACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCAttt                                                                                             <  1:1010661/150‑1 (MQ=255)
                                                             tttAGCCAAGTGGAGAAAAGGGAATGTGGGAGCTATGACACAGAGAGAAAAGAACTAGAGAATATAAAACGCCCCTGCCAGCAACTGACAGGTGCGTACGTACCCGATAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCt                                                                                        <  2:2173733/150‑1 (MQ=255)
                                                             tttAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCt                                                                                        >  1:2158461/1‑150 (MQ=255)
                                                                aGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGc                                                                                                                      >  1:873591/1‑117 (MQ=255)
                                                                aGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGc                                                                                                                      <  2:873591/117‑1 (MQ=255)
                                                                 gCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGc                                                                                    >  2:48424/1‑150 (MQ=255)
                                                                  ccAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCa                                                                                   >  1:1894889/1‑150 (MQ=255)
                                                                   cAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAg                                                                                  <  2:2073784/150‑1 (MQ=255)
                                                                   cAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGAACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAg                                                                                  >  2:1209216/1‑150 (MQ=255)
                                                                    aaGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAgg                                                                                 <  1:2449088/150‑1 (MQ=255)
                                                                         gagaAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCaaaa                                                                            >  1:1232376/1‑150 (MQ=255)
                                                                         gagaAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCaaaa                                                                            >  2:1684411/1‑150 (MQ=255)
                                                                           gaAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGt                                                                          >  2:111541/1‑150 (MQ=255)
                                                                           gaAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGATTCGCAGGCAAAAGt                                                                          >  2:1510139/1‑150 (MQ=255)
                                                                              aaGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGaa                                                                       >  1:938375/1‑150 (MQ=255)
                                                                                gggAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAAt                                                                     >  1:1658360/1‑150 (MQ=255)
                                                                                  gAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCg                                                                   >  2:2054668/1‑150 (MQ=255)
                                                                                     aGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGata                                                                <  1:212895/150‑1 (MQ=255)
                                                                                        gggAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACGCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCTCATGCAAAAGTGAAATCGATATcc                                                             <  2:1212283/150‑1 (MQ=255)
                                                                                           aGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTTCGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAgg                                                          <  1:1868992/150‑1 (MQ=255)
                                                                                           aGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAgg                                                          <  2:177083/150‑1 (MQ=255)
                                                                                           aGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAgg                                                          >  1:2301285/1‑150 (MQ=255)
                                                                                              tATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGgcgc                                                       >  1:1068267/1‑150 (MQ=255)
                                                                                              tATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGgcgc                                                       <  2:1880708/150‑1 (MQ=255)
                                                                                              tATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGgcgc                                                       >  2:1854799/1‑150 (MQ=255)
                                                                                              tATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGgcac                                                       >  1:1320507/1‑148 (MQ=255)
                                                                                                 gACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCTGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGcaca                                                    >  1:1086231/1‑150 (MQ=255)
                                                                                                  acacAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCacac                                                   >  1:2076786/1‑150 (MQ=255)
                                                                                                     caGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGt                                                >  1:565590/1‑150 (MQ=255)
                                                                                                      agagagAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTc                                               >  1:1860774/1‑150 (MQ=255)
                                                                                                        agagaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAg                                             >  1:2485165/1‑150 (MQ=255)
                                                                                                         gagaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGa                                            >  1:276192/1‑150 (MQ=255)
                                                                                                          agaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGAc                                           <  2:210411/150‑1 (MQ=255)
                                                                                                             aaaGAAGAATAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGActct                                        <  1:1333867/150‑1 (MQ=255)
                                                                                                             aaaGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGActct                                        <  2:1947670/150‑1 (MQ=255)
                                                                                                             aaaGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGActct                                        >  1:1016508/1‑150 (MQ=255)
                                                                                                             aaaGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGActct                                        <  1:468839/150‑1 (MQ=255)
                                                                                                                   gaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTg                                  >  1:1772910/1‑150 (MQ=255)
                                                                                                                      gagaAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGaaa                               >  1:508095/1‑150 (MQ=255)
                                                                                                                       agaaaagaaaACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAg                              <  2:2355297/150‑1 (MQ=255)
                                                                                                                       agaaaagaaaACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAg                              <  2:1384350/150‑1 (MQ=255)
                                                                                                                          aaagaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCg                           >  2:326700/1‑150 (MQ=255)
                                                                                                                           aagaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCgg                          <  1:2025356/150‑1 (MQ=255)
                                                                                                                             gaaaaCGCCCCTTCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCgggg                        <  1:944658/150‑1 (MQ=255)
                                                                                                                             gaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCgggg                        <  1:2259422/150‑1 (MQ=255)
                                                                                                                             gaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCgggg                        <  2:86295/150‑1 (MQ=255)
                                                                                                                              aaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGt                       >  1:279926/1‑150 (MQ=255)
                                                                                                                                taCGCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGa                     >  1:1965425/2‑150 (MQ=255)
                                                                                                                                aaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGa                     >  2:754169/1‑150 (MQ=255)
                                                                                                                                  cGCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtt                   >  2:2009465/1‑150 (MQ=255)
                                                                                                                                  cGCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtt                   >  2:200274/1‑150 (MQ=255)
                                                                                                                                  cGCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtt                   >  1:356261/1‑150 (MQ=255)
                                                                                                                                  cGCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtt                   >  2:46211/1‑150 (MQ=255)
                                                                                                                                  cGCCCCTGTCAGTTGCTGGCAGGGGCCTACGTACCCGAGAGGAATTAACGCAGAACCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtt                   >  2:1791180/1‑150 (MQ=255)
                                                                                                                                   gCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGGCAGACTCTTCAGTGAAAGGCGGGGTGAttt                  >  2:1578593/1‑150 (MQ=255)
                                                                                                                                   gCCCCTGTCAGTAGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAttt                  >  1:767270/1‑150 (MQ=255)
                                                                                                                                       cTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCa              <  2:1392580/150‑1 (MQ=255)
                                                                                                                                         gCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAg            <  1:2272993/150‑1 (MQ=255)
                                                                                                                                          ccAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAgg           >  1:2440620/1‑150 (MQ=255)
                                                                                                                                           cAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGt          >  2:2296556/1‑150 (MQ=255)
                                                                                                                                                aCTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGa     >  2:1425395/1‑150 (MQ=255)
                                                                                                                                                 cTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGAt    <  1:1791180/150‑1 (MQ=255)
                                                                                                                                                 cTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGAt    >  1:1652834/1‑150 (MQ=255)
                                                                                                                                                  tGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGAtt   >  1:1444005/1‑150 (MQ=255)
                                                                                                                                                   gACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTc  >  1:1570332/1‑150 (MQ=255)
                                                                                                                                                   gACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTc  >  1:1956731/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTC  >  NC_000913/2460476‑2460772

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: