Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 985,853 C→T 6.3% 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 6.3% 251.0 / 5.2 79intergenic (‑144/+41)aspC/ompFaspartate aminotransferase, PLP‑dependent/outer membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base C (34/40);  new base T (2/3);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

ACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGT  >  NC_000913/985706‑986000
                                                                                                                                                    |                                                                                                                                                   
aCATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCa                                                                                                                                                    >  4:145340/1‑149 (MQ=255)
  aTGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCaaa                                                                                                                                                  <  3:619356/149‑1 (MQ=255)
       gAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCt                                                                                                                                             >  1:685221/1‑149 (MQ=255)
         ggTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGATCAAAGTCCTGt                                                                                                                                           <  4:316634/149‑1 (MQ=255)
         ggTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGt                                                                                                                                           <  1:260110/149‑1 (MQ=255)
          gTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGtt                                                                                                                                          <  4:454534/149‑1 (MQ=255)
             ccATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACTTTGGTCCTGttttt                                                                                                                                       <  1:482895/149‑1 (MQ=255)
              cATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGtttttt                                                                                                                                      >  3:672873/1‑149 (MQ=255)
               aTTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTc                                                                                                                                     <  2:625738/149‑1 (MQ=255)
               aTTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTc                                                                                                                                     <  4:381642/149‑1 (MQ=255)
                ttATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCg                                                                                                                                    >  3:310720/1‑149 (MQ=255)
                 tATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACTTTGGTCCTGtttttttta                                                                                                                                   <  3:378057/149‑4 (MQ=255)
                  aTGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGc                                                                                                                                  >  3:413160/1‑149 (MQ=255)
                      ttACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCAttt                                                                                                                              >  2:650941/1‑149 (MQ=255)
                          aGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACa                                                                                                                          >  3:132624/1‑149 (MQ=255)
                           gAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACaa                                                                                                                         >  4:454/1‑149 (MQ=255)
                            aaGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACaaa                                                                                                                        >  3:565265/1‑149 (MQ=255)
                            aaGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACaaa                                                                                                                        >  1:400397/1‑149 (MQ=255)
                            aaGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACaaa                                                                                                                        <  4:565265/149‑1 (MQ=255)
                              gggAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAaga                                                                                                                      <  3:309533/149‑1 (MQ=255)
                                 aaGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgt                                                                                                                   <  3:124387/149‑1 (MQ=255)
                                    tCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGGCGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGtgtg                                                                                                                >  2:142749/1‑149 (MQ=255)
                                     ccGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAACTCCTGTTTTTTCTGCATTTAACAAAGAGGTGTGc                                                                                                               <  4:244059/149‑1 (MQ=255)
                                          aTCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCaaa                                                                                                                                                  <  2:171164/109‑1 (MQ=255)
                                          aTCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCaaa                                                                                                                                                  >  1:171164/1‑109 (MQ=255)
                                          aTCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTa                                                                                                          <  2:272214/149‑1 (MQ=255)
                                          aTCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCAGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTa                                                                                                          <  1:385947/149‑1 (MQ=255)
                                           tCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAg                                                                                                         >  1:40780/1‑149 (MQ=255)
                                           tCAGGGAAAAGGGAGATTTACAAAATTCCAACTATTACTGAAGAAAACGCAGGCTGTTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTATCGGCATTTAACAAAGAGGTGTGCTATTa                                                                                                          <  1:498917/149‑1 (MQ=255)
                                             aGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGaa                                                                                                       <  2:37172/149‑1 (MQ=255)
                                                     gggAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaaa                                                                                               >  2:182779/1‑149 (MQ=255)
                                                     gggAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaaa                                                                                               <  2:437333/149‑1 (MQ=255)
                                                          aTTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATa                                                                                          <  3:614460/149‑1 (MQ=255)
                                                                aaaaTTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCcaca                                                                                    <  4:310720/149‑1 (MQ=255)
                                                                aaaaTTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCcaca                                                                                    <  4:378147/149‑1 (MQ=255)
                                                                  aaTTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACTTTGGTCCTGttt                                                                                                                                         <  4:497968/94‑1 (MQ=39)
                                                                       cAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACgg                                                                             >  3:113316/1‑149 (MQ=255)
                                                                        aaCTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACTTTGGTCCTGttttt                                                                                                                                       >  4:547697/1‑90 (MQ=38)
                                                                        aaCTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGgt                                                                            >  3:243218/1‑149 (MQ=255)
                                                                                 tgatgaAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtgt                                                                                                                 >  3:603869/1‑103 (MQ=255)
                                                                                 tgatgaAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtgt                                                                                                                 <  4:603869/103‑1 (MQ=255)
                                                                                 tgatgaAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGa                                                                   >  2:433283/1‑149 (MQ=255)
                                                                                  gatgaAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGaa                                                                  <  4:672873/149‑1 (MQ=255)
                                                                                   atgaAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAAc                                                                 <  1:329112/149‑1 (MQ=255)
                                                                                     gaAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgt                                                                                                                   >  2:63369/1‑97 (MQ=255)
                                                                                     gaAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgt                                                                                                                   <  1:63369/97‑1 (MQ=255)
                                                                                        aaCGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACg                                                            >  4:155826/1‑149 (MQ=255)
                                                                                         aCGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGc                                                           >  4:354778/1‑149 (MQ=255)
                                                                                             aGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAg                                                       >  1:109868/1‑149 (MQ=255)
                                                                                             aGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAg                                                       >  2:601489/1‑149 (MQ=255)
                                                                                                 tGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTttg                                                   >  3:340769/1‑149 (MQ=255)
                                                                                                  gTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTttgt                                                  >  1:438411/1‑149 (MQ=255)
                                                                                                        tGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCa                                             <  2:297772/149‑1 (MQ=255)
                                                                                                             gACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATc                                        <  1:433283/149‑1 (MQ=255)
                                                                                                             gACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGAACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATc                                        >  4:645900/1‑149 (MQ=255)
                                                                                                                gTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACACAGAGATGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGAt                                     <  3:494286/149‑1 (MQ=255)
                                                                                                                 tGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATc                                    <  2:528018/149‑1 (MQ=255)
                                                                                                                     aTTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGt                                <  1:51652/149‑1 (MQ=255)
                                                                                                                       tGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtg                              <  3:645900/149‑1 (MQ=255)
                                                                                                                         ctctGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtgt                                                                                                                 <  1:202414/64‑1 (MQ=255)
                                                                                                                         ctctGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtgt                                                                                                                 >  2:202414/1‑64 (MQ=255)
                                                                                                                          tctGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatg                           >  1:430494/1‑149 (MQ=255)
                                                                                                                            tGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTgatgat                         >  4:65247/1‑149 (MQ=255)
                                                                                                                             ggAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGaa                                          <  2:306899/131‑1 (MQ=255)
                                                                                                                             ggAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGaa                                          >  1:306899/1‑131 (MQ=255)
                                                                                                                              gAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGt                       <  1:128907/149‑1 (MQ=255)
                                                                                                                                   tataAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCa                  >  4:11183/1‑149 (MQ=255)
                                                                                                                                    ataAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCaa                 <  2:109868/149‑1 (MQ=255)
                                                                                                                                    ataAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCaa                 <  3:65247/149‑1 (MQ=255)
                                                                                                                                    ataAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTATCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCaa                 <  1:601489/149‑1 (MQ=255)
                                                                                                                                      cgaaaaaaCAGGACTTTGGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACa               >  1:515262/3‑149 (MQ=255)
                                                                                                                                      aaaaaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACa               <  2:369223/149‑1 (MQ=255)
                                                                                                                                       aaaaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatg                           <  1:528209/136‑1 (MQ=255)
                                                                                                                                       aaaaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatg                           >  2:528209/1‑136 (MQ=255)
                                                                                                                                       aaaaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACAt              <  4:132624/149‑1 (MQ=255)
                                                                                                                                       aaaaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACAt              >  2:536368/1‑149 (MQ=255)
                                                                                                                                          aaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAgg           >  2:480950/1‑149 (MQ=255)
                                                                                                                                            aaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTg         <  3:106763/149‑1 (MQ=255)
                                                                                                                                               aGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGAc      <  3:499063/149‑1 (MQ=255)
                                                                                                                                                   ccAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGt  <  1:142749/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGT  >  NC_000913/985706‑986000

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: