Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,322,398 C→G 7.1% V183V (GTG→GTC)  trpE ← component I of anthranilate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,322,3980CG 7.1% 295.8 / 11.1 85V183V (GTG→GTC) trpEcomponent I of anthranilate synthase
Reads supporting (aligned to +/- strand):  ref base C (39/40);  new base G (4/2);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

GCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTC  >  NC_000913/1322250‑1322545
                                                                                                                                                    |                                                                                                                                                   
gCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATc                                                                                                                                                     <  1:640436/149‑1 (MQ=255)
 cGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCa                                                                                                                                                    <  2:989731/149‑1 (MQ=255)
   gCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCAcc                                                                                                                                                  <  4:1155441/149‑1 (MQ=255)
   gCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCAcc                                                                                                                                                  >  2:17483/1‑149 (MQ=255)
            ccACTGGCAGCGGCGGCGCGGCTTCGGTCGGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAg                                                                                                                                             <  3:696438/145‑1 (MQ=255)
            ccACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAg                                                                                                                                             >  4:696438/1‑145 (MQ=255)
                          ggcgcggcTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTc                                                                                                                                       <  4:1293006/137‑1 (MQ=255)
                          ggcgcggcTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTc                                                                                                                                       >  3:1293006/1‑137 (MQ=255)
                           gcgcggcTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATaaaaa                                                                                                                          <  2:1248884/149‑1 (MQ=255)
                           gcgcggcTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATaaaaa                                                                                                                          <  1:652684/149‑1 (MQ=255)
                                  ttCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGa                                                                                                                      >  2:1295667/1‑146 (MQ=255)
                                  ttCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGa                                                                                                                      <  1:1295667/146‑1 (MQ=255)
                                     ggTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAg                                                                                                                <  3:1419375/149‑1 (MQ=255)
                                     ggTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAg                                                                                                                <  2:334136/149‑1 (MQ=255)
                                               tGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTAtt                                                                                                      <  1:17483/149‑1 (MQ=255)
                                                      gttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAGCAGAAATCAGGGCAGTTATTTTCCGCt                                                                                               >  3:770753/1‑149 (MQ=255)
                                                      gttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGtt                                                                                                         >  2:590830/1‑139 (MQ=255)
                                                      gttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGtt                                                                                                         <  1:590830/139‑1 (MQ=255)
                                                              aGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTg                                                                                       >  3:1096284/1‑149 (MQ=255)
                                                              aGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTg                                                                                       >  3:190333/1‑149 (MQ=255)
                                                                     cAGTGAGACGTAGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTaaa                                                                                >  3:860055/1‑149 (MQ=255)
                                                                       gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATc                                                                              <  4:190333/149‑1 (MQ=255)
                                                                          agaCGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAGCAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTc                                                                           <  4:770753/149‑1 (MQ=255)
                                                                              gttgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCgaga                                                                                                                                >  1:965980/1‑92 (MQ=255)
                                                                              gttgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCgaga                                                                                                                                <  2:965980/92‑1 (MQ=255)
                                                                              gttgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGc                                                                                      >  1:800859/1‑134 (MQ=255)
                                                                              gttgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGc                                                                                      <  2:800859/134‑1 (MQ=255)
                                                                               ttgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACa                                                                                           >  1:787064/1‑128 (MQ=255)
                                                                               ttgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACa                                                                                           <  2:787064/128‑1 (MQ=255)
                                                                               ttgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATc                                                                      >  1:1436690/1‑149 (MQ=255)
                                                                               ttgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATc                                                                      >  4:107872/1‑149 (MQ=255)
                                                                                 gttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCga                                                                                                                                  <  4:817038/87‑1 (MQ=255)
                                                                                 gttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCga                                                                                                                                  >  3:817038/1‑87 (MQ=255)
                                                                                   ttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGc                                                                  <  3:107872/149‑1 (MQ=255)
                                                                                         tcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAg                                                            <  2:1436690/149‑1 (MQ=255)
                                                                                         tcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAg                                                            <  4:860055/149‑1 (MQ=255)
                                                                                          cttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAgg                                                           >  3:920533/1‑149 (MQ=255)
                                                                                           ttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGc                                                                                                <  2:1222531/111‑1 (MQ=255)
                                                                                           ttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGc                                                                                                >  1:1222531/1‑111 (MQ=255)
                                                                                           ttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGt                                                          >  4:1264528/1‑149 (MQ=255)
                                                                                           ttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGt                                                          >  3:123197/1‑149 (MQ=255)
                                                                                            tcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCt                                                                             <  3:802720/129‑1 (MQ=255)
                                                                                            tcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCt                                                                             >  4:802720/1‑129 (MQ=255)
                                                                                                 cGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGc                                                                  <  2:553392/135‑1 (MQ=255)
                                                                                                 cGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGc                                                                  >  1:553392/1‑135 (MQ=255)
                                                                                                    aGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCTGTAAATCTTCa                                                                          <  3:836103/124‑1 (MQ=255)
                                                                                                    aGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCa                                                                          >  4:836103/1‑124 (MQ=255)
                                                                                                        aaCAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAg                                             <  4:1096284/149‑1 (MQ=255)
                                                                                                             gCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAgg                                                                                                               <  2:147609/78‑1 (MQ=255)
                                                                                                             gCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAgg                                                                                                               >  1:147609/1‑78 (MQ=255)
                                                                                                                 gCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCgaa                                    <  1:1288192/149‑1 (MQ=255)
                                                                                                                 gCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCgaa                                    >  2:1288192/1‑149 (MQ=255)
                                                                                                                   cTGAATACGGGTGCTTTTTTTCTGATGGTCAATCAGCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGtg                                                                                        >  2:780577/1‑94 (MQ=255)
                                                                                                                   cTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGtt                                                                                        <  1:780577/95‑1 (MQ=255)
                                                                                                                   cTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCtt                                                                            >  1:57129/1‑107 (MQ=255)
                                                                                                                   cTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCtt                                                                            <  2:57129/107‑1 (MQ=255)
                                                                                                                    tGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTAtttt                                                                                                    >  1:736428/1‑82 (MQ=255)
                                                                                                                    tGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTAtttt                                                                                                    <  2:736428/82‑1 (MQ=255)
                                                                                                                        tACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATg                             <  3:1264528/149‑1 (MQ=255)
                                                                                                                                 tttttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGt                    >  4:26972/1‑149 (MQ=255)
                                                                                                                                 tttttttCTGATGGTCAATCACCATAAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGt                    >  3:150/1‑149 (MQ=255)
                                                                                                                                         tGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTtctc                      >  1:520407/1‑139 (MQ=255)
                                                                                                                                         tGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTtctc                      <  2:520407/139‑1 (MQ=255)
                                                                                                                                         tGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttcctt            >  4:866410/1‑149 (MQ=255)
                                                                                                                                         tGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttcctt            >  4:1328958/1‑149 (MQ=255)
                                                                                                                                         tGATGGTCAATCACCATCAGCGTTTCAGCGAAATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttcctt            <  4:123197/149‑1 (MQ=255)
                                                                                                                                          gATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAAc                                               >  2:288998/1‑113 (MQ=255)
                                                                                                                                          gATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAAc                                               <  1:288998/113‑1 (MQ=255)
                                                                                                                                            tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac                                                               <  2:489455/95‑1 (MQ=255)
                                                                                                                                            tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac                                                               >  1:489455/1‑95 (MQ=255)
                                                                                                                                            tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAgg                                                           <  4:667396/99‑1 (MQ=255)
                                                                                                                                            tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAgg                                                           >  3:667396/1‑99 (MQ=255)
                                                                                                                                             ggTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGt        >  4:104045/1‑149 (MQ=255)
                                                                                                                                               tCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTaaa                                                                                >  4:1192302/1‑75 (MQ=255)
                                                                                                                                               tCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTaaa                                                                                <  3:1192302/75‑1 (MQ=255)
                                                                                                                                                 aaTCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGc                                                                                      >  2:466181/1‑67 (MQ=255)
                                                                                                                                                 aaTCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGc                                                                                      <  1:466181/67‑1 (MQ=255)
                                                                                                                                                 aaTCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAaga                                                   >  2:232904/1‑102 (MQ=255)
                                                                                                                                                 aaTCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAaga                                                   <  1:232904/102‑1 (MQ=255)
                                                                                                                                                  aTCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttccttc           >  3:789321/1‑141 (MQ=255)
                                                                                                                                                  aTCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCAGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCtt                         <  3:450650/127‑1 (MQ=255)
                                                                                                                                                  aTCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCAGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCtt                         >  4:450650/1‑127 (MQ=255)
                                                                                                                                                  aTCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCACTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttccttc           <  4:789321/141‑1 (MQ=255)
                                                                                                                                                   tCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACAt                              <  3:625221/121‑1 (MQ=255)
                                                                                                                                                   tCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACAt                              >  4:625221/1‑121 (MQ=255)
                                                                                                                                                   tCACCATCAGCGTTTCAGAGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTc  >  4:1406616/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTC  >  NC_000913/1322250‑1322545

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