Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,460,620 A→G 6.8% intergenic (‑151/+30) fadI ← / ← yfcZ beta‑ketoacyl‑CoA thiolase, anaerobic, subunit/UPF0381 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,460,6200AG 6.8% 333.3 / 11.9 118intergenic (‑151/+30)fadI/yfcZbeta‑ketoacyl‑CoA thiolase, anaerobic, subunit/UPF0381 family protein
Reads supporting (aligned to +/- strand):  ref base A (55/55);  new base G (5/3);  total (60/58)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

Only 100 of 123 total aligned reads displayed.

ACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAG  >  NC_000913/2460471‑2460768
                                                                                                                                                     |                                                                                                                                                    
         tGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAggg                                                                                                                                             >  1:214226/1‑150 (MQ=255)
          taaaaaaaaaaaaaaaAGTGTTCTCACCTGAAAATAGACTAAACCATTTTTTAAAATTTAGCCAAGTGGACAAAAGGTAAAGTGGTAGATATGACACAGAGATAAAATAACAAGAGAAAAGAAAACGCCCCTGTCAGTTGCTGGCAgggg                                                                                                                                            <  2:350927/148‑1 (MQ=255)
          gTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAgggg                                                                                                                                            <  1:950341/150‑1 (MQ=255)
            aaaaaaaTAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcg                                                                                                                                          <  2:295066/150‑1 (MQ=255)
            aaaaaaaTAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTAACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcg                                                                                                                                          >  1:1349291/1‑150 (MQ=255)
               aaaaTAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgtac                                                                                                                                       >  2:763297/1‑150 (MQ=255)
                aaataaatAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGTCAGTTGCTGGCAGGGGcgtttt                                                                                                                                      <  2:592526/150‑4 (MQ=255)
                  ataaataaCTTGGTCTGACCTGCTCATAGTCTTAACCATTTTTTTACATTTCGCCAAGTGGAGAAACGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgtacgta                                                                                                                                    <  2:1472143/150‑1 (MQ=255)
                  ataaataaAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGTCAGTTGCTGGCAGGGGcgttttct                                                                                                                                    <  2:1009592/150‑6 (MQ=255)
                      ataaaGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTAACATTTACCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAACGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCg                                                                                                                                <  1:1395267/150‑1 (MQ=255)
                         aaGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCgaga                                                                                                                             >  1:629664/1‑150 (MQ=255)
                          agaggTCTGACCTGATCATAGTCTTATCCATTTTTTTACATTTAGCCAAGTGGAGAAAATGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTAGCCTagag                                                                                                                            <  2:94294/147‑1 (MQ=255)
                          aGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGagag                                                                                                                            <  1:255605/150‑1 (MQ=255)
                                cTGACCTGATCTTATTCTTAACCATTTTTATACATTTAGGCAAGTGGAGAAAAGGGAAAGTGGGTGCAATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTa                                                                                                                      <  2:1257424/150‑1 (MQ=255)
                                cTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTa                                                                                                                      >  2:1211432/1‑150 (MQ=255)
                                  taCCTGATCAATGTCTAAACCATTTTTTTACATTTAGCCACGGGTAGAAAAGGTAAAGGGGTCGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAAc                                                                                                                    <  1:879588/149‑1 (MQ=255)
                                   aCCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACg                                                                                                                   >  2:1513881/1‑150 (MQ=255)
                                     cTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAATGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCa                                                                                                                 <  2:1068637/150‑1 (MQ=255)
                                     cTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCa                                                                                                                 <  2:503985/150‑1 (MQ=255)
                                     cTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGCGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCa                                                                                                                 >  1:546502/1‑150 (MQ=255)
                                     cTGATCATAGTCTTAACAATTTTTTTACATTTAGCCACGTGTAGAAAAGGGAATTTGGGAGCTATGACACAGAGAGAAAAGAAGAAGATAAAAGAAAACGCCACTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCa                                                                                                                 <  2:810533/150‑1 (MQ=255)
                                      tGATCTAAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAATAGGGAAAGTGGGATCTATGAAACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAg                                                                                                                <  1:837376/150‑1 (MQ=255)
                                      tGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGccg                                                                                                                >  2:1285298/1‑148 (MQ=255)
                                       gATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGa                                                                                                               >  2:1368478/1‑150 (MQ=255)
                                       gATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGa                                                                                                               >  1:1009845/1‑150 (MQ=255)
                                        aTCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGAc                                                                                                              <  2:214226/150‑1 (MQ=255)
                                         tCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGAcc                                                                                                             >  2:1114664/1‑150 (MQ=255)
                                          aattttCTTAACCTTTTTTTTAAATTTAGCAAATTGTATAAAAGGGAAAGTGGTAGCTATGACACAGATATAAAATAATAAGATAAAATAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGAccc                                                                                                            <  2:1661975/145‑1 (MQ=255)
                                           aTAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCa                                                                                                           >  2:1471384/1‑150 (MQ=255)
                                              gTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCt                                                                                                        <  1:763297/150‑1 (MQ=255)
                                                  tAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGCCAGGGCCGTACGTACACGAGAGGAATTAACGCAGACCCAGCTGGaa                                                                                                    >  1:241202/1‑150 (MQ=255)
                                                  tAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGaa                                                                                                    >  1:183424/1‑150 (MQ=255)
                                                    aCCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGaaaa                                                                                                  >  1:1688332/1‑150 (MQ=255)
                                                    aCCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGaaaa                                                                                                  >  2:1449969/1‑150 (MQ=255)
                                                     ccATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAAt                                                                                                 >  1:108231/1‑150 (MQ=255)
                                                      cATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATc                                                                                                >  1:311028/1‑150 (MQ=255)
                                                         ttttttACATTTAGCCATGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAATAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGc                                                                                             <  2:1718774/150‑1 (MQ=255)
                                                         tttttaAAATTTAGCCAATTGGATAAAAGGTAAAGTGGGAGCTATGACACATATATAAAATAAGAACACAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGc                                                                                             <  1:1437861/150‑1 (MQ=255)
                                                           ttttACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACTGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCAt                                                                                           <  2:776930/150‑1 (MQ=255)
                                                                cATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGCCCCAGCTGGAAAATCAGCATTTCAg                                                                                      >  1:1041783/1‑150 (MQ=255)
                                                                 aTTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGc                                                                                     <  1:588270/150‑1 (MQ=255)
                                                                     aGCCTAGTGGAGAAACGGGACGGTGGGAGCTATCACACAGAGAGAACACAAGACGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCg                                                                                 <  1:585149/150‑1 (MQ=255)
                                                                      gCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGc                                                                                <  1:941090/150‑1 (MQ=255)
                                                                         aaGTGGAGAAAAGGTAAAGTGGGAGCTATGACACAGAGATAAAATAATAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAgg                                                                             <  2:1442148/150‑1 (MQ=255)
                                                                          aGTGTATAAAAGGTAAAGTGGTAGCAATGACACTCAGAGAAAAGAAGAAGATAAAAGAACACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGc                                                                            <  2:1551690/150‑1 (MQ=255)
                                                                              gagaAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCATCAATTCAGCTTCCCAGGCaaaa                                                                        >  1:839359/1‑150 (MQ=255)
                                                                              gagaAAAGGGAAAGTGCGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATATCAGCTTCGCAGGCaaaa                                                                        >  1:1136402/1‑150 (MQ=255)
                                                                                gaAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGt                                                                      >  1:819777/1‑150 (MQ=255)
                                                                                 aaaaGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTg                                                                     <  1:1021293/150‑1 (MQ=255)
                                                                                  aaaGGGAATGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGGAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGCCCCAGCGGGAAAATCAGCATTTCACCGTCGCCGGCAAAAGTgg                                                                    >  1:301192/1‑149 (MQ=255)
                                                                                    aGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAAAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCTTTTCAGCTTCGCAGGCAAAAGTGaaa                                                                  >  2:1419769/1‑150 (MQ=255)
                                                                                     gggAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAAt                                                                 >  2:1476393/1‑150 (MQ=255)
                                                                                      ggAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGCCCCAGCAGGAAAATCAGCATTTCAGCGTCGCAGGCAAAAGTGAAATc                                                                >  1:1063753/1‑150 (MQ=255)
                                                                                       gAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCg                                                               <  2:1054679/150‑1 (MQ=255)
                                                                                          aGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTAAGCTTCGCAGGCAAAAGTGAAATCGata                                                            >  2:192735/1‑150 (MQ=255)
                                                                                             gggAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGCCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATcc                                                         >  2:844939/1‑150 (MQ=255)
                                                                                             gggAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATcc                                                         <  1:1114664/150‑1 (MQ=255)
                                                                                              ggAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCa                                                        <  1:497664/150‑1 (MQ=255)
                                                                                                aGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAgg                                                      <  2:981288/150‑1 (MQ=255)
                                                                                                  cTATGACACAGAGAGAAAAGAAGAAGAGAAACGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGgcg                                                    <  2:826980/150‑1 (MQ=255)
                                                                                                  cTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGgcg                                                    <  1:414171/150‑1 (MQ=255)
                                                                                                   tATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGgcgc                                                   >  2:1008711/1‑150 (MQ=255)
                                                                                                   tATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACAGACAGGGGCGTACGTACCCGAGAGCAATTAACGCAGACCCAGCTGGCAAATCAGCATATCACCATCGCAGGAAAAAGTGAAATCGATATCCAGgcgc                                                   >  2:1002264/1‑150 (MQ=255)
                                                                                                   aaTGACAAAGAGAGAATAGAAGAAGAGAAAAGAAAACGGCCCTGCCAGCAACTGACAGGGGCGTACGAACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGgcgc                                                   <  2:831031/149‑1 (MQ=255)
                                                                                                          caGAGAGCAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGCCCCAGCTGGAAAATCAGCATTTCCGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGt                                            >  1:881357/1‑150 (MQ=255)
                                                                                                          caGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCCGCTTCGCAGGCAAAAGTGAAATCGATTTCCAGGCGCACACCGt                                            >  1:875278/1‑150 (MQ=255)
                                                                                                            gagagaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCa                                          <  2:684775/150‑1 (MQ=255)
                                                                                                            gagagaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCa                                          <  2:1719541/150‑1 (MQ=255)
                                                                                                             agagaaaagaagaagagaaaagaagaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAg                                         <  1:136431/150‑1 (MQ=255)
                                                                                                             agagaAATGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAg                                         <  1:432721/150‑1 (MQ=255)
                                                                                                              gagaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGCCCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGa                                        >  1:994718/1‑150 (MQ=255)
                                                                                                               agaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGAc                                       <  1:864389/150‑1 (MQ=255)
                                                                                                                gaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGAct                                      >  1:710565/1‑150 (MQ=255)
                                                                                                                gaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGAct                                      <  2:603555/150‑1 (MQ=255)
                                                                                                                gaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGAct                                      <  2:238943/150‑1 (MQ=255)
                                                                                                                gaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGAct                                      >  2:982299/1‑150 (MQ=255)
                                                                                                                    agaagaagaGAAAAGAAAACGCCCCTGCCAGCAACTGCCAGGGGCGTACGTACCCGAGAGGAATAAACGCCGCCCCAGCTGGAAAATCAGCATTTCAGGTGCGCAGGCAAAAGTGACATCGATATCCAGACGCACACCGTCAGACTCTTc                                  >  1:87275/1‑150 (MQ=255)
                                                                                                                    agaagaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGTGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTc                                  <  2:370528/150‑1 (MQ=255)
                                                                                                                      aagaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAg                                >  2:1217039/1‑150 (MQ=255)
                                                                                                                       agaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGt                               <  1:971345/150‑1 (MQ=255)
                                                                                                                                aagaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACACGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCgg                      >  2:62550/1‑150 (MQ=255)
                                                                                                                                     aaCGCCCCTGCCCGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCAGGAAAATCAGCATTACAGCTACGCAGGAAAAAGTGAGATCGATATCGAGGCGCACACGGTCAGCCTCTGCAGTGAAAGTCGGGGTGa                 >  1:89758/1‑150 (MQ=255)
                                                                                                                                      aCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAt                <  2:311028/150‑1 (MQ=255)
                                                                                                                                       cGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtt               >  2:821126/1‑150 (MQ=255)
                                                                                                                                       cGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtt               >  1:448073/1‑150 (MQ=255)
                                                                                                                                        gCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAttt              >  1:1174847/1‑150 (MQ=255)
                                                                                                                                         ccccTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTACCGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGagtgt             >  1:244624/1‑146 (MQ=255)
                                                                                                                                         ccccTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGTTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtttt             >  2:819098/1‑150 (MQ=255)
                                                                                                                                         ccccTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGTATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtttt             >  2:1448377/1‑150 (MQ=255)
                                                                                                                                         ccccTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtttt             >  1:1265143/1‑150 (MQ=255)
                                                                                                                                         ccccTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCGTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTGAGTGAAAGTCGGGGTGAtttt             >  1:1634396/1‑150 (MQ=255)
                                                                                                                                         ccccTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtttt             >  1:1241775/1‑150 (MQ=255)
                                                                                                                                          cccTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTg            <  1:212169/150‑1 (MQ=255)
                                                                                                                                           ccTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGc           >  2:938089/1‑150 (MQ=255)
                                                                                                                                             tGCCAGCTACTGTCAGGGTCGTACGTTCCCGAGAGGATTTAACGGAGACCCTGCTGGATAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAATTCGATATCCAGGCGCACACCGTCTGACTCTTCAGTGAAAGTCGGGGTGATTTTGCaa         <  1:1002264/150‑1 (MQ=255)
                                                                                                                                             tGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCaa         <  1:1008711/150‑1 (MQ=255)
                                                                                                                                                cAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACCCCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGt      >  2:985722/1‑150 (MQ=255)
                                                                                                                                                cAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGt      <  2:108231/150‑1 (MQ=255)
                                                                                                                                                    aaCTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAg  <  2:241202/150‑1 (MQ=255)
                                                                                                                                                    aaCTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAg  <  1:1449969/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAG  >  NC_000913/2460471‑2460768

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