Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 293 A→G intergenic (+38/‑44) thrL → / → thrA thr operon leader peptide/aspartokinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132930AG99.7% 1000.1 / ‑6.0 317intergenic (+38/‑44)thrL/thrAthr operon leader peptide/aspartokinase I
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (125/191);  total (126/191)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 343 total aligned reads displayed.

GCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAC  >  NC_000913/203‑382
                                                                                          |                                                                                          
gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGt                                                                                          <  1:6379821/93‑1 (MQ=255)
acaccacCATTACCCCCCCCATCACCATTACAACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGt                                                                                          <  2:7367443/92‑1 (MQ=255)
   ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCg                                                                                       <  2:4511212/93‑1 (MQ=255)
   ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCg                                                                                       <  2:4767897/93‑1 (MQ=255)
   ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCg                                                                                       >  1:8037358/1‑93 (MQ=255)
   ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCg                                                                                       >  2:7753603/1‑93 (MQ=255)
   ccaccaTAACCCCCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCg                                                                                       <  2:819256/93‑1 (MQ=255)
     cccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCggg                                                                                     <  1:3593610/92‑1 (MQ=255)
     accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGGGCggg                                                                                     >  1:9419874/1‑93 (MQ=255)
       caTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCt                                                                                  >  2:7935243/1‑93 (MQ=255)
        aTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtt                                                                                 >  1:7401727/1‑93 (MQ=255)
        aTTACCAACACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTAAAGGAAACACAAAAAAAAGCCCGCACCTGACGGTGCGGGCtt                                                                                 <  1:8982622/93‑1 (MQ=255)
         ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttt                                                                                <  2:3177732/93‑1 (MQ=255)
         ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttt                                                                                <  1:6291165/93‑1 (MQ=255)
           accaccaccATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttttt                                                                              <  2:2284744/93‑1 (MQ=255)
             caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttttttt                                                                            >  1:3241633/1‑93 (MQ=255)
              accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtttttttt                                                                           <  2:1520261/93‑1 (MQ=255)
              accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtttttttt                                                                           >  2:6717544/1‑93 (MQ=255)
                caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTc                                                                         <  2:6474435/93‑1 (MQ=255)
                 accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGGCGGTGCGGTCTTTTTTTTTCg                                                                        >  1:8611131/1‑93 (MQ=255)
                 accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCg                                                                        <  1:10005434/93‑1 (MQ=255)
                  ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtttttttttttg                                                                       >  2:4839869/1‑90 (MQ=255)
                    atcaccatTACCACAGGTAACGGTGCGGGTTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACGGTGCGGGCTTTTTTTTTCGAcc                                                                     <  2:3071328/93‑1 (MQ=255)
                     tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCa                                                                    >  2:1755294/1‑93 (MQ=255)
                       cccatTACCACAGGTAACGGTGGGGGCTGTCGCGTACAGGAAACGCAGAAAAAAGCCCCCCCCTGACGGTGCGGGCTTTTTTTTTCGACCaaa                                                                  <  2:6009703/92‑1 (MQ=255)
                       accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCaaa                                                                  >  1:5184612/1‑93 (MQ=255)
                         catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAgg                                                                <  1:7071814/93‑1 (MQ=255)
                           ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCAACTGACGGTGCGGGGGTTTTTTTTTGTCCaaaagga                                                              >  2:3865106/1‑89 (MQ=255)
                            tACCACAGGTAACGGTGCGGTCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTaa                                                             <  2:7011657/93‑1 (MQ=255)
                            tACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTaa                                                             <  1:1375387/93‑1 (MQ=255)
                            tACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTaa                                                             <  1:2271891/93‑1 (MQ=255)
                             aCCACAGGTACCGGTGCGGGCTCACGCGTACAGGAAACACAGAAAAACGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAAc                                                            <  2:7487773/93‑1 (MQ=255)
                             aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAGGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAAc                                                            <  2:4853906/93‑1 (MQ=255)
                             aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCTGGCTTTTTTTTTCGACCAAAGGGAAc                                                            >  2:4462342/1‑93 (MQ=255)
                              ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTCGACCAAAGGTAACGa                                                          <  2:4324693/93‑1 (MQ=255)
                               cacaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGa                                                          <  1:6672412/93‑1 (MQ=255)
                                acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTTGACCAAAAGTGACGAg                                                         >  2:4107004/1‑93 (MQ=255)
                                 agggTAACGGGGCGGGCTGACGCGTACAGAAAACAAAGAAAAAACCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg                                                        <  2:8826692/91‑1 (MQ=255)
                                       aCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                  <  2:2825559/93‑1 (MQ=255)
                                       aCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                  <  1:46601/93‑1 (MQ=255)
                                       aCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                  >  2:2555940/1‑93 (MQ=255)
                                         ggTGGGGGCTGAAGGGAACGGAAACCACAGAAAAAACCCCCCCCCGGCCGGTGCGGGTTTTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                 <  2:9372555/93‑1 (MQ=255)
                                         ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAAcc                                                <  1:9221972/93‑1 (MQ=255)
                                          gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                               <  2:9620056/93‑1 (MQ=255)
                                             cGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                            <  1:3442589/93‑1 (MQ=255)
                                               gccTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                          <  1:3810013/91‑1 (MQ=255)
                                                gCTGCCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                         <  1:407381/93‑1 (MQ=255)
                                                 cTGACGCGTACAGAAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                                        <  2:5760147/93‑1 (MQ=255)
                                                  tGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtg                                       >  1:4967807/1‑93 (MQ=255)
                                                      gcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                   <  1:12174/93‑1 (MQ=255)
                                                      gcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                   <  1:3307990/93‑1 (MQ=255)
                                                      gcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                   <  1:2816474/93‑1 (MQ=255)
                                                        gTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                 >  1:9514440/1‑93 (MQ=255)
                                                         tACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGt                                >  2:8314971/1‑93 (MQ=255)
                                                          accggAAAAACAGAAAAAAGCCCGCCCCTGACGGTGCGGGCTTTTTTTGTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGtt                               <  1:7616599/90‑1 (MQ=255)
                                                          aCAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGtt                               <  1:5508758/93‑1 (MQ=255)
                                                          aCAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGtt                               <  2:730510/93‑1 (MQ=255)
                                                          aCAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGtt                               <  1:1660358/93‑1 (MQ=255)
                                                           cAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                              <  1:483353/93‑1 (MQ=255)
                                                            aGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                             <  1:5913818/93‑1 (MQ=255)
                                                            aGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                             <  2:3303612/93‑1 (MQ=255)
                                                            aGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTca                             >  2:6807566/1‑92 (MQ=255)
                                                             ggAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                            <  2:1285701/93‑1 (MQ=255)
                                                             ggAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                            <  1:7940796/93‑1 (MQ=255)
                                                             gaaaaCACAGAAAAAAGCCCGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                            <  2:93143/91‑1 (MQ=255)
                                                              gAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                           <  1:2503199/93‑1 (MQ=255)
                                                               aaaGACAGAAAATAGCCCTCACCTGGCGGTGCGGTCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcg                          <  1:4151982/93‑1 (MQ=255)
                                                                aaCACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcgg                         <  2:5897239/93‑1 (MQ=255)
                                                                  cacaGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                       <  2:8755930/93‑1 (MQ=255)
                                                                  cacaGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                       <  1:8562804/93‑1 (MQ=255)
                                                                   acaGAAAAAAGCCGGCACCTGACGGTGCGGGCTTTTTTTTGCGACCAAAGGTAACGAGGTAACACCCATGCGAGTGTTGAAGTTCGGCGGTAc                      <  1:5209806/93‑1 (MQ=255)
                                                                   acaGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                      >  2:9778590/1‑93 (MQ=255)
                                                                   acaGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACa                     <  1:2371549/93‑1 (MQ=255)
                                                                    caGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACa                     <  2:3644611/93‑1 (MQ=255)
                                                                    caGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACa                     <  1:4744462/93‑1 (MQ=255)
                                                                      gAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                   <  2:4645871/93‑1 (MQ=255)
                                                                      gAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                   >  1:737074/1‑93 (MQ=255)
                                                                      gAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                   <  1:5341418/93‑1 (MQ=255)
                                                                       aaaaagCCCCGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCa                  <  2:1943672/93‑1 (MQ=255)
                                                                       aaaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCa                  <  2:2607201/93‑1 (MQ=255)
                                                                        aaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAgg                >  2:6311440/1‑92 (MQ=255)
                                                                         aaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTTGACCAAAAGTAACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGGGACATCAgg                >  2:5703072/1‑92 (MQ=255)
                                                                           aaccccGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg              <  2:7437161/90‑1 (MQ=255)
                                                                            aGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc             >  1:1759520/1‑93 (MQ=255)
                                                                              cccGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGGAACGAGGTAACAACCATGCGAGGGTTGAAGGTCGGCGGGACATCAGTGGCaa           >  2:5910952/1‑93 (MQ=255)
                                                                               ccGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaaa          <  1:3097432/93‑1 (MQ=255)
                                                                               ccGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAAATCAGTGGCaaa          >  1:9783375/1‑93 (MQ=255)
                                                                                 gCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg        <  2:6574555/93‑1 (MQ=255)
                                                                                  cACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc       <  1:5712554/93‑1 (MQ=255)
                                                                                  cACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc       <  1:6503162/93‑1 (MQ=255)
                                                                                  cACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc       <  2:6142151/93‑1 (MQ=255)
                                                                                  cACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc       >  2:2860183/1‑93 (MQ=255)
                                                                                   aCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa      >  1:2746049/1‑93 (MQ=255)
                                                                                    ccTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg     <  1:4468550/93‑1 (MQ=255)
                                                                                    ccTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg     >  1:8135146/1‑93 (MQ=255)
                                                                                     cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTGACGAGGTAACAACCATGCGAGTGTTGAAGTTCAGCGGTACATCAGTGGCAAATGCAGa    >  2:1845999/1‑93 (MQ=255)
                                                                                     cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa    <  2:8482038/93‑1 (MQ=255)
                                                                                     cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa    <  1:8118300/93‑1 (MQ=255)
                                                                                     cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa    >  1:8358981/1‑93 (MQ=255)
                                                                                       gACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc  >  1:4332700/1‑93 (MQ=255)
                                                                                          |                                                                                          
GCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAC  >  NC_000913/203‑382

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: