Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,993,145 C→A 100% R665S (CGT→AGT)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,993,1450CA98.8% 282.7 / ‑3.4 80R665S (CGT→AGT) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (44/35);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAA  >  NC_000913/3992997‑3993290
                                                                                                                                                    |                                                                                                                                                 
gatgatCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGc                                                                                                                                                   <  2:583439/149‑1 (MQ=255)
      cGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCttt                                                                                                                                             >  1:464387/1‑149 (MQ=255)
         aGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCa                                                                                                                                          <  3:61780/149‑1 (MQ=255)
         aGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCa                                                                                                                                          <  1:441702/149‑1 (MQ=255)
              tGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCAcc                                                                                                                                     <  2:329375/149‑1 (MQ=255)
              tGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCAcc                                                                                                                                     >  1:234214/1‑149 (MQ=255)
                 aaaCTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGc                                                                                                                                  <  4:331368/149‑1 (MQ=255)
                   aCTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGTACCCGCCa                                                                                                                                <  3:629209/149‑1 (MQ=255)
                     tATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAgg                                                                                                                              >  2:11550/1‑149 (MQ=255)
                     tATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAgg                                                                                                                              >  4:53534/1‑149 (MQ=255)
                     tATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAgg                                                                                                                              >  4:226342/1‑149 (MQ=255)
                          tcGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAAcc                                                                                                                         >  4:535482/1‑149 (MQ=255)
                             gCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCggg                                                                                                                      >  1:65566/1‑149 (MQ=255)
                              cAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGc                                                                                                                     >  4:245403/1‑149 (MQ=255)
                                     cATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCaa                                                                                                              <  4:505411/149‑1 (MQ=255)
                                     cATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCaa                                                                                                              <  4:539576/149‑1 (MQ=255)
                                     cATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCaa                                                                                                              <  3:440687/149‑1 (MQ=255)
                                      aTCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAg                                                                                                             <  1:194201/149‑1 (MQ=255)
                                         aGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGgcg                                                                                                          >  3:288418/1‑149 (MQ=255)
                                          ggACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGgcgc                                                                                                         >  3:37835/1‑149 (MQ=255)
                                                 gcACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGTCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGt                                                                                                  <  2:65566/149‑1 (MQ=255)
                                                      gCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTg                                                                                             >  2:83892/1‑149 (MQ=255)
                                                      gCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACAGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTg                                                                                             <  3:465136/149‑1 (MQ=255)
                                                       ccAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTgg                                                                                            <  2:584439/149‑1 (MQ=255)
                                                          gATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATGGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCa                                                                                         >  2:49956/1‑149 (MQ=255)
                                                            tAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCaaa                                                                                       <  4:672590/149‑1 (MQ=255)
                                                             aGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAAGCGCTGCGCGTTTCTGGTCAAAc                                                                                      >  4:433100/1‑149 (MQ=255)
                                                              gCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAAcc                                                                                     >  2:470532/1‑149 (MQ=255)
                                                                gTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTg                                                                                   <  3:272460/149‑1 (MQ=255)
                                                                gTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTg                                                                                   <  1:470532/149‑1 (MQ=255)
                                                                gTGGAAGTCTTCTGTTAAAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTg                                                                                   >  1:102289/1‑149 (MQ=255)
                                                                   gAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTgggg                                                                                <  4:185531/149‑1 (MQ=255)
                                                                   gAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTgggg                                                                                >  2:38767/1‑149 (MQ=255)
                                                                   gAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTgggg                                                                                >  2:351493/1‑149 (MQ=255)
                                                                   gAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTgggg                                                                                >  1:104849/1‑149 (MQ=255)
                                                                    aaGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTggggg                                                                               >  3:318720/1‑149 (MQ=255)
                                                                    aaGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTggggg                                                                               <  3:296947/149‑1 (MQ=255)
                                                                    aaGTCTTCTGTTATAGCCAGCATCCGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTggggg                                                                               >  1:595692/1‑149 (MQ=255)
                                                                     aGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgt                                                                              <  3:22174/149‑1 (MQ=255)
                                                                         ttctGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgttgtt                                                                          <  3:509/149‑1 (MQ=255)
                                                                         ttctGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGGTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgttgtt                                                                          >  1:545350/1‑149 (MQ=255)
                                                                            tGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGttctt                                                                       <  1:26089/149‑1 (MQ=255)
                                                                             gTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGttcttc                                                                      <  3:552784/149‑1 (MQ=255)
                                                                              ttATAGCCAGTATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCg                                                                     <  2:240077/149‑1 (MQ=255)
                                                                              ttATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCg                                                                     >  4:298387/1‑149 (MQ=255)
                                                                                 taGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCttt                                                                                                                                             >  4:250746/1‑74 (MQ=255)
                                                                                 taGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCttt                                                                                                                                             <  3:250746/74‑1 (MQ=255)
                                                                                 taGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAAc                                                                  <  2:234214/149‑1 (MQ=255)
                                                                                   gCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGc                                                                >  4:665982/1‑149 (MQ=255)
                                                                                    ccAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGcc                                                               >  2:417603/1‑149 (MQ=255)
                                                                                     cAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCt                                                              <  3:330168/149‑1 (MQ=255)
                                                                                         aTCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAAt                                                          <  1:203458/149‑1 (MQ=255)
                                                                                         aTCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAAt                                                          >  2:486416/1‑149 (MQ=255)
                                                                                             gcgcGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTAt                                                      <  3:114539/149‑1 (MQ=255)
                                                                                               gcgGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCg                                                    >  1:266047/1‑149 (MQ=255)
                                                                                                    ttAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACa                                               >  2:268556/1‑149 (MQ=255)
                                                                                                       aTTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGaa                                            >  3:558012/1‑149 (MQ=255)
                                                                                                        ttCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGaaa                                           >  2:266820/1‑149 (MQ=255)
                                                                                                        ttCGTACTCGCGTGCAGCAACTGGTGGCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGTGCGTTTCTGGTCAAACCTGGGGATTGTTCTTCGAACGCCTGAATGTATCGGTACAGaaa                                           >  4:231161/1‑149 (MQ=255)
                                                                                                            tACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTgg                                       <  2:104849/149‑1 (MQ=255)
                                                                                                               tCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTa                                                 >  1:414315/1‑136 (MQ=255)
                                                                                                               tCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTa                                                 <  2:414315/136‑1 (MQ=255)
                                                                                                                   gTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGc                                >  4:577125/1‑149 (MQ=255)
                                                                                                                        gcaACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCg                           >  1:474350/1‑149 (MQ=255)
                                                                                                                           aCTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGc                                                                                                                     >  1:665426/1‑56 (MQ=255)
                                                                                                                           aCTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGc                                                                                                                     <  2:665426/56‑1 (MQ=255)
                                                                                                                           aCTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGt                        <  4:37835/149‑1 (MQ=255)
                                                                                                                            cTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGtt                       >  4:455417/1‑149 (MQ=255)
                                                                                                                              ggTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGttta                     >  1:620574/1‑148 (MQ=255)
                                                                                                                                tttCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTAt                   <  4:496688/149‑1 (MQ=255)
                                                                                                                                 ttCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATg                  <  1:486416/149‑1 (MQ=255)
                                                                                                                                 ttCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATg                  <  4:558012/149‑1 (MQ=255)
                                                                                                                                     gAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGgcgc              >  1:140109/1‑149 (MQ=255)
                                                                                                                                          tATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGAttt         >  1:538914/1‑149 (MQ=255)
                                                                                                                                            ttGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCg       >  2:661427/1‑149 (MQ=255)
                                                                                                                                            ttGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCg       >  1:48561/1‑149 (MQ=255)
                                                                                                                                            ttGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCg       >  1:235131/1‑149 (MQ=255)
                                                                                                                                            ttGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCg       >  1:123499/1‑149 (MQ=255)
                                                                                                                                               aaTTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCAt    >  3:401636/1‑149 (MQ=255)
                                                                                                                                                 ttGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATaa  <  1:96635/149‑1 (MQ=255)
                                                                                                                                                 ttGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATaa  >  1:497058/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAA  >  NC_000913/3992997‑3993290

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