Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 985,853 C→T 7.7% intergenic (‑144/+41) aspC ← / ← ompF aspartate aminotransferase, PLP‑dependent/outer membrane porin 1a (Ia;b;F)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913985,8530CT 7.7% 274.2 / 8.1 90intergenic (‑144/+41)aspC/ompFaspartate aminotransferase, PLP‑dependent/outer membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base C (45/38);  new base T (3/4);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

AAACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGT  >  NC_000913/985704‑986000
                                                                                                                                                      |                                                                                                                                                   
aaaCATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGAcc                                                                                                                                                     >  2:1322848/1‑150 (MQ=255)
 aaCATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCa                                                                                                                                                    >  1:404942/1‑150 (MQ=255)
   cATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCaaa                                                                                                                                                  >  2:198076/1‑150 (MQ=255)
    aTGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAg                                                                                                                                                 <  1:691876/150‑1 (MQ=255)
     tGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGt                                                                                                                                                >  2:1528948/1‑150 (MQ=255)
     tGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGt                                                                                                                                                >  1:1685205/1‑150 (MQ=255)
      gACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAAACAGGACCAAAGt                                                                                                                                                >  1:1565947/1‑150 (MQ=255)
       aCGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTcc                                                                                                                                              <  1:962661/150‑1 (MQ=255)
         gatgTTCCATTATGGTTACAGAAGCGAAGTCCGCTATCAGGGTAACTGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTTTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTg                                                                                                                                            <  1:121577/147‑1 (MQ=255)
         gAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTg                                                                                                                                            <  1:185329/150‑1 (MQ=255)
            gTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACTTTGGTCCTGttt                                                                                                                                         <  1:1307257/150‑1 (MQ=255)
               ccATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACTTTGGTCCTGtttttt                                                                                                                                      <  1:81579/150‑1 (MQ=255)
                cATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACTTTGGTCCTGttttttt                                                                                                                                     <  2:375244/150‑2 (MQ=255)
                 aTTATGGTTACAGAAGGGAAGTCCGCTATCATGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCg                                                                                                                                    <  1:232859/150‑1 (MQ=255)
                 aTTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCg                                                                                                                                    <  1:658632/150‑1 (MQ=255)
                  ttATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCgg                                                                                                                                   >  2:1417004/1‑150 (MQ=255)
                     tGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGTGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGTCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCAt                                                                                                                                <  1:1452016/150‑1 (MQ=255)
                      ggTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCAtt                                                                                                                               >  1:1214616/1‑150 (MQ=255)
                       gTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGt                                                                                                                                           <  1:887905/137‑1 (MQ=255)
                       gTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGt                                                                                                                                           >  2:887905/1‑137 (MQ=255)
                           cAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACa                                                                                                                          >  1:363045/1‑150 (MQ=255)
                             gAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACaaa                                                                                                                        >  2:946133/1‑150 (MQ=255)
                             gAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACaaa                                                                                                                        >  1:117909/1‑150 (MQ=255)
                                gggAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAagag                                                                                                                     <  1:520463/150‑1 (MQ=255)
                                  gAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgt                                                                                                                   >  1:1682708/1‑150 (MQ=255)
                                   aaGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCAAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtg                                                                                                                  <  1:1522425/150‑1 (MQ=255)
                                     gTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTTCTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTATCAAAGAGGtgtg                                                                                                                <  1:490448/150‑1 (MQ=255)
                                       ccGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCt                                                                                                              >  1:955221/1‑150 (MQ=255)
                                       ccGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCt                                                                                                              >  2:1186958/1‑150 (MQ=255)
                                       ccGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCt                                                                                                              >  1:1665153/1‑150 (MQ=255)
                                         gCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGAAGGATGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCGTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTAt                                                                                                            >  1:851321/1‑150 (MQ=255)
                                           tATCAGTGTAACGGGAGATTTACAAAAATCCAACTATAACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTa                                                                                                          <  2:209792/150‑1 (MQ=255)
                                           tATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCGAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTa                                                                                                          <  2:542100/150‑1 (MQ=255)
                                           tATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTa                                                                                                          >  2:1572239/1‑150 (MQ=255)
                                            aTCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGTAGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAg                                                                                                         <  1:151521/150‑1 (MQ=255)
                                               aGGGTTACTGTAGATTTACAAAATTCCAACTATTACTGATGAAAACGCTGGCTTTTTTTGCAAGACGTGAGTTTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAAc                                                                                                      <  2:340753/150‑1 (MQ=255)
                                               aGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAAc                                                                                                      <  1:198076/150‑1 (MQ=255)
                                                 ggTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTg                                                                                                    >  2:1129238/1‑150 (MQ=255)
                                                  gTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTATCAAAGAGGTGTGCTATTAGAACTgg                                                                                                   >  1:1551034/1‑150 (MQ=255)
                                                  gTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTgg                                                                                                   <  2:637433/150‑1 (MQ=255)
                                                  gTAAAGGGAGATTTTCAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTTGAAGGGAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTgg                                                                                                   <  2:302338/150‑1 (MQ=255)
                                                   ttaCGGGAGTTTTTCAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAAGAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGt                                                                                                  <  1:1239840/148‑1 (MQ=255)
                                                        ggAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtg                                                                                                                  <  2:1465033/129‑1 (MQ=255)
                                                        ggAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtg                                                                                                                  >  1:1465033/1‑129 (MQ=255)
                                                        ggAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCGGCTTTTGCAAGCCGGGAGATTGCTCTGGAAGGTAT‑AAAAAAAACACGCCCAAAGTCCTGTTTTTTCGGCATTTAACAAAGCGGCGTGCTATTAGAACTGGTAAACg                                                                                             <  2:1169644/150‑1 (MQ=255)
                                                            aTTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATAc                                                                                         <  1:1417004/150‑1 (MQ=255)
                                                               tACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCca                                                                                      <  2:1230494/150‑1 (MQ=255)
                                                                 cAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCcaca                                                                                    >  2:1697129/1‑150 (MQ=255)
                                                                  aaaaTTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAg                                                                                   >  2:1534912/1‑150 (MQ=255)
                                                                   aaaTTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGc                                                                                  >  1:635002/1‑150 (MQ=255)
                                                                     aTTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCaa                                                                                <  1:946133/150‑1 (MQ=255)
                                                                           aCTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGgtgt                                                                          >  1:935888/1‑150 (MQ=255)
                                                                           aCTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGgtgt                                                                          <  1:1528948/150‑1 (MQ=255)
                                                                            cTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTgtc                                                                         >  2:490716/1‑150 (MQ=255)
                                                                             tATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTgtcg                                                                        <  2:955221/150‑1 (MQ=255)
                                                                                    gatgaAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAAc                                                                 >  1:1066458/1‑150 (MQ=255)
                                                                                          aaCGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGaa                                                                  <  1:1213512/143‑1 (MQ=255)
                                                                                          aaCGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTAGTCTGaa                                                                  >  2:1213512/1‑143 (MQ=255)
                                                                                            cGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGccc                                                         >  2:773086/1‑150 (MQ=255)
                                                                                              cAGGCTGTTTTTGCAAGCCGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATAAGAACTGCTAAACGATACCCACAGAAAAGGTGTCGTCTCAACCTACGGCCAg                                                       >  1:1611113/1‑150 (MQ=255)
                                                                                                  ctctttttGCATGACGTGAGATTGCTCTGGAAGGAAT‑AAAAAAAACAGGACCACAGTCCGGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTttg                                                   <  1:1091109/147‑1 (MQ=255)
                                                                                                  cTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTttg                                                   <  1:566665/150‑1 (MQ=255)
                                                                                                       tttGCATGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTc                                              <  2:1123930/150‑1 (MQ=255)
                                                                                                        ttGCAAGTCGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCa                                             <  1:1350407/150‑1 (MQ=255)
                                                                                                         tgtaaGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAg                                            <  2:1570889/147‑1 (MQ=255)
                                                                                                          gCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTAGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGACCCTACGCCCAGTTTGTTGTCAGa                                           >  2:1001354/1‑150 (MQ=255)
                                                                                                            aaGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTAAGCCCAGTTTGTTGTCAGAAt                                         >  1:1002161/1‑150 (MQ=255)
                                                                                                             aGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATc                                        <  2:802193/150‑1 (MQ=255)
                                                                                                                     gaTTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGt                                >  1:590974/1‑150 (MQ=255)
                                                                                                                          ctctGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatg                           >  2:1353705/1‑150 (MQ=255)
                                                                                                                          ctctGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatg                           >  2:746964/1‑150 (MQ=255)
                                                                                                                          ctctGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatg                           >  2:265323/1‑150 (MQ=255)
                                                                                                                          ctctGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatg                           <  2:635002/150‑1 (MQ=255)
                                                                                                                          ctctGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGATATTAGAACTGGTAAAAGATACCAACAGCAACGGTGTCGTCTGAACCTACGCCAAGTTTGTTGTCAGAATCGATCTGGTtgatt                           >  2:1598679/1‑149 (MQ=255)
                                                                                                                           tctGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatga                          <  2:117909/150‑1 (MQ=255)
                                                                                                                            ctGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTgatgat                         >  2:242133/1‑150 (MQ=255)
                                                                                                                               gAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTa                      >  1:1712227/1‑150 (MQ=255)
                                                                                                                                   gTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGTATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCa                  >  1:887291/1‑150 (MQ=255)
                                                                                                                                   gTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCa                  >  1:1543227/1‑150 (MQ=255)
                                                                                                                                    tat‑aAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCaa                 >  1:1247470/1‑150 (MQ=255)
                                                                                                                                     at‑aAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAAc                <  1:488018/150‑1 (MQ=255)
                                                                                                                                     at‑aAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAAc                <  2:404942/150‑1 (MQ=255)
                                                                                                                                      gccgaaaaaaCAGGACTTTGGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGtt                               <  1:149550/131‑1 (MQ=255)
                                                                                                                                      gccgaaaaaaCAGGACTTTGGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGtt                               >  2:149550/5‑135 (MQ=255)
                                                                                                                                       ccgaaaaaaCAGGACTGTGGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACa               >  1:1305977/4‑150 (MQ=255)
                                                                                                                                        aaaaaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATAAACACAGCAACAGTGACGTCTGAACCTACGCCAACTATGTAGTCAGAATCGATCTGCTTAATTATGAAGTCAACAt              >  2:1257451/1‑150 (MQ=255)
                                                                                                                                         aaaaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATa             >  2:1710536/1‑150 (MQ=255)
                                                                                                                                            aaaaCAGGACTTTGGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATTATGTAGTCAACATAGGt          >  1:1703146/1‑150 (MQ=255)
                                                                                                                                              aaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTgg        <  1:490716/150‑1 (MQ=255)
                                                                                                                                                  ggACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACAt    >  1:396296/1‑150 (MQ=255)
                                                                                                                                                    aCCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGt  >  2:753067/1‑150 (MQ=255)
                                                                                                                                                     ccAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGatgatg                        <  2:337759/127‑1 (MQ=255)
                                                                                                                                                     ccAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGatgatg                        >  1:337759/1‑127 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
AAACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGT  >  NC_000913/985704‑986000

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