Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 783,000 C→A 5.3% Q306K (CAG→AAG)  nadA → quinolinate synthase, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913783,0000CA 5.3% 670.1 / 17.1 187Q306K (CAG→AAG) nadAquinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base C (92/85);  new base A (3/7);  total (95/92)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 191 total aligned reads displayed.

GAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTAAC  >  NC_000913/782852‑783146
                                                                                                                                                    |                                                                                                                                                  
gagGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTc                                                                                                                                                    >  3:944832/1‑149 (MQ=255)
gagGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTc                                                                                                                                                    >  1:1153813/1‑149 (MQ=255)
  ggCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAg                                                                                                                                                  <  1:558910/149‑1 (MQ=255)
   gCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAgg                                                                                                                                                 >  1:1345823/1‑149 (MQ=255)
     ttATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGcc                                                                                                                                               >  2:1154236/1‑149 (MQ=255)
          gtgGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTGAAGGCCAttca                                                                                                                                          <  1:945983/149‑4 (MQ=255)
           tgGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCa                                                                                                                                         >  1:1157864/1‑149 (MQ=255)
           tgGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTGAAGGCCAttcat                                                                                                                                         <  2:169754/149‑5 (MQ=255)
            ggCAACCGATCGGTGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTGCGTGGATGGCCATGAATGGCCTGAAGGCCagtcatg                                                                                                                                        <  2:1083749/149‑7 (MQ=255)
            ggCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTGAAGGCCAttcatg                                                                                                                                        <  1:1255070/149‑6 (MQ=255)
              cAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCagag                                                                                                                                      <  3:266329/149‑1 (MQ=255)
               aaCCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGGATGGCCTGAAGGCCAttcatg                                                                                                                                        >  4:1114807/1‑141 (MQ=255)
               aaCCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGGATGGCCTGAAGGCCAttcatg                                                                                                                                        <  3:1114807/146‑6 (MQ=255)
               aaCCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAgg                                                                                                                                     >  2:934125/1‑149 (MQ=255)
                    aTCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTa                                                                                                                                <  4:1035213/149‑1 (MQ=255)
                       ggggTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGaa                                                                                                                             <  2:521831/149‑1 (MQ=255)
                        gggTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAAc                                                                                                                            >  2:632886/1‑149 (MQ=255)
                           tATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACagg                                                                                                                         >  1:1100912/1‑149 (MQ=255)
                            aTTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACagga                                                                                                                        <  2:57519/149‑1 (MQ=255)
                             ttttCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACaggaa                                                                                                                       >  4:140716/1‑149 (MQ=255)
                             ttttCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACaggaa                                                                                                                       >  3:500959/1‑149 (MQ=255)
                             ttttCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACaggaa                                                                                                                       <  3:572299/149‑1 (MQ=255)
                                  tACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGaa                                                                                                                             <  4:1220631/138‑1 (MQ=255)
                                  tACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGaa                                                                                                                             >  3:1220631/1‑138 (MQ=255)
                                  tACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGgaaggaag                                                                                                                  >  3:1006494/1‑149 (MQ=255)
                                     aaaaTGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCaa                                                                                                               >  4:127868/1‑149 (MQ=255)
                                       aaTGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATc                                                                                                             >  3:635554/1‑149 (MQ=255)
                                        aTGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCa                                                                                                            >  4:367541/1‑149 (MQ=255)
                                           cagcagGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGa                                                                                                         <  3:1323557/149‑1 (MQ=255)
                                           cagcagGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGa                                                                                                         <  3:148588/149‑1 (MQ=255)
                                            agcagGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAg                                                                                                        <  2:4661/149‑1 (MQ=255)
                                            agcagGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAg                                                                                                        >  2:1417297/1‑149 (MQ=255)
                                             gcagGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTGAAGGCCAttca                                                                                                                                          >  2:1113097/1‑111 (MQ=255)
                                              cagGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGt                                                                                                      <  4:633809/149‑1 (MQ=255)
                                                  cGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCAt                                                                                                  <  1:632886/149‑1 (MQ=255)
                                                       ccAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTtga                                                                                             <  1:353305/149‑1 (MQ=255)
                                                       ccAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTtga                                                                                             <  1:896687/149‑1 (MQ=255)
                                                       ccAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTtga                                                                                             <  2:511347/149‑1 (MQ=255)
                                                       ccAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAACGAAGCAATCACGAGGTTCATGTtga                                                                                             <  4:1176035/149‑1 (MQ=255)
                                                        cAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTtgat                                                                                            <  1:1136068/149‑1 (MQ=255)
                                                          gATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTtgatga                                                                                          >  1:57010/1‑149 (MQ=255)
                                                            tAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAg                                                                                                                                                  <  4:698876/91‑1 (MQ=255)
                                                            tAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGaaa                                                                                        >  4:638719/1‑149 (MQ=255)
                                                             aaaGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGTCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAg                                                                                       <  2:1389717/149‑1 (MQ=255)
                                                              aaGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAgg                                                                                      <  1:970644/149‑1 (MQ=255)
                                                                  gTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGc                                                                                  >  3:789803/1‑149 (MQ=255)
                                                                   ttACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAggaaggaa                                                                                                                   >  4:1204439/1‑115 (MQ=255)
                                                                   ttACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCg                                                                                 >  4:554986/1‑149 (MQ=255)
                                                                      cTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCgaga                                                                              <  4:490843/149‑1 (MQ=255)
                                                                      cTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCgaga                                                                              <  3:1264170/149‑1 (MQ=255)
                                                                        ggAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCgagaga                                                                            <  1:1408120/149‑1 (MQ=255)
                                                                         gAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGagagag                                                                           >  3:1203628/1‑149 (MQ=255)
                                                                         gAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGagagag                                                                           >  1:1108815/1‑149 (MQ=255)
                                                                         gAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGagagag                                                                           >  2:896327/1‑149 (MQ=255)
                                                                            gCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGgc                                                                        >  1:204954/1‑149 (MQ=255)
                                                                               ccAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAggaaggaa                                                                                                                   <  3:720603/103‑1 (MQ=255)
                                                                                   ccGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTg                                                                 <  1:305797/149‑1 (MQ=255)
                                                                                         gTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTc                                                           >  2:27703/1‑149 (MQ=255)
                                                                                             gggTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATc                                                       <  2:994052/149‑1 (MQ=255)
                                                                                                tGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAg                                                                                                        <  1:216272/97‑1 (MQ=255)
                                                                                                tGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAg                                                                                                        >  2:216272/1‑97 (MQ=255)
                                                                                                tGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTa                                                    <  1:268236/149‑1 (MQ=255)
                                                                                                tGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTa                                                    <  2:874214/149‑1 (MQ=255)
                                                                                                 gCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTAt                                                   <  2:340404/149‑1 (MQ=255)
                                                                                                 gCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTAt                                                   >  3:522322/1‑149 (MQ=255)
                                                                                                  cAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATg                                                  <  1:723749/149‑1 (MQ=255)
                                                                                                   aaCCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGc                                                 >  3:674643/1‑149 (MQ=255)
                                                                                                   aaCCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGc                                                 >  2:1176436/1‑149 (MQ=255)
                                                                                                    aCCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCt                                                >  4:75943/1‑149 (MQ=255)
                                                                                                      cTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTgg                                              >  3:565041/1‑149 (MQ=255)
                                                                                                       tGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGa                                             <  4:223448/149‑1 (MQ=255)
                                                                                                       tGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGa                                             <  4:674643/149‑1 (MQ=255)
                                                                                                           gCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCCGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGAtttt                                         >  3:1181193/1‑149 (MQ=255)
                                                                                                               cTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCgg                                     <  4:722212/149‑1 (MQ=255)
                                                                                                                   gcgcATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTac                                 >  4:1364733/1‑149 (MQ=255)
                                                                                                                   gcgcATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTac                                 >  1:40302/1‑149 (MQ=255)
                                                                                                                     gcATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTacac                               <  2:1141024/149‑1 (MQ=255)
                                                                                                                        ttGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTAc                            >  1:1330320/1‑149 (MQ=255)
                                                                                                                           tccatcgATGGCCATGAATGGCCTTCATGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTg                         <  4:941076/143‑1 (MQ=255)
                                                                                                                            ccGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTgg                        >  3:1049295/1‑149 (MQ=255)
                                                                                                                              gTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGAt                      >  1:1331688/1‑149 (MQ=255)
                                                                                                                                tgaaggCCATTCATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATaa                    >  2:869808/5‑149 (MQ=255)
                                                                                                                                tgaaggCCATTCATGGCCATCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTg                         <  3:958638/140‑1 (MQ=255)
                                                                                                                                    ggCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGaa                <  4:565041/149‑1 (MQ=255)
                                                                                                                                       cATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGaataa             >  4:749169/1‑149 (MQ=255)
                                                                                                                                       cATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGaataa             >  2:425039/1‑149 (MQ=255)
                                                                                                                                       cATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGaataa             >  1:1368345/1‑149 (MQ=255)
                                                                                                                                        aTGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGaataat            <  1:1176436/149‑1 (MQ=255)
                                                                                                                                           gaTGGCCTGAAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAg         >  4:1174575/2‑149 (MQ=255)
                                                                                                                                           aaTGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAg         <  2:1426378/149‑1 (MQ=255)
                                                                                                                                           aaTGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAg         <  3:749169/149‑1 (MQ=255)
                                                                                                                                             tGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGc       <  1:1298698/149‑1 (MQ=255)
                                                                                                                                             tGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGc       <  1:449228/149‑1 (MQ=255)
                                                                                                                                               gCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGt     >  1:361029/1‑149 (MQ=255)
                                                                                                                                                ccTTCAGGCCATCGCGGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAAc                   <  1:794650/134‑1 (MQ=255)
                                                                                                                                                 cTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTaa   <  4:789803/149‑1 (MQ=255)
                                                                                                                                                  ttCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGc                                       >  4:1183621/1‑112 (MQ=255)
                                                                                                                                                  ttCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTAAc  >  1:1218956/1‑149 (MQ=255)
                                                                                                                                                  ttCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTAAc  >  3:1406407/1‑149 (MQ=255)
                                                                                                                                                  ttCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTAAc  >  4:543455/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTAAC  >  NC_000913/782852‑783146

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: