Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,406,098 C→T 5.8% D40D (GAC→GAT)  smrA → DNA endonuclease SmrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,406,0980CT 5.8% 272.0 / 4.3 86D40D (GAC→GAT) smrADNA endonuclease SmrA
Reads supporting (aligned to +/- strand):  ref base C (47/34);  new base T (3/2);  total (50/36)
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.60e-01

ATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCT  >  NC_000913/1405949‑1406247
                                                                                                                                                     |                                                                                                                                                     
atCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGac                                                                                                                                                       >  2:870218/1‑150 (MQ=255)
atCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGac                                                                                                                                                       >  2:125881/1‑150 (MQ=255)
      ctctTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTg                                                                                                                                                 >  1:1007363/1‑150 (MQ=255)
       tctTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGc                                                                                                                                                <  1:1342511/150‑1 (MQ=255)
       tctTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGc                                                                                                                                                >  2:740994/1‑150 (MQ=255)
        ctTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCa                                                                                                                                               >  2:337704/1‑150 (MQ=255)
         ttAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                                                                              <  1:687408/150‑3 (MQ=255)
          tAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGc                                                                                                                                             <  1:1542552/150‑1 (MQ=255)
            aTTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCtt                                                                                                                                           >  1:110172/1‑150 (MQ=255)
             ttGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTg                                                                                                                                          <  2:767602/150‑1 (MQ=255)
             ttGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTg                                                                                                                                          <  2:743584/150‑1 (MQ=255)
             ttGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTg                                                                                                                                          <  1:1486062/150‑1 (MQ=255)
                  ggaggaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAAt                                                                                                                                     >  1:646201/1‑150 (MQ=255)
                         gTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCAGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCACAGTGCGCCGAAGCGTATCGACACGCTGCAGCTTGATAATTTCCTca                                                                                                                              >  1:1605235/1‑150 (MQ=255)
                           aaCATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTcacc                                                                                                                            <  1:1033643/150‑1 (MQ=255)
                           aaCATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTcacc                                                                                                                            >  2:776197/1‑150 (MQ=255)
                               tGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCg                                                                                                                        >  2:364457/1‑150 (MQ=255)
                                     ttGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTtc                                                                                                                  >  2:5875/1‑150 (MQ=255)
                                          gacAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAc                                                                                                             <  2:29248/150‑1 (MQ=255)
                                          gacAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAc                                                                                                             <  2:239261/150‑1 (MQ=255)
                                          gacAAATCGCTGTTTCTTGACGCCAAGGAAGATGACCAGCCGCTGAAACGTGCTACCGATGTCAACTGGCATACAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCATGATAATTACCTAACCACCGGATTTCTCGAc                                                                                                             >  2:1072522/1‑150 (MQ=255)
                                           acAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAca                                                                                                            >  1:1069463/1‑150 (MQ=255)
                                                 cGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATcc                                                                                                      >  1:216891/1‑150 (MQ=255)
                                                   cTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATcccc                                                                                                    >  2:408584/1‑149 (MQ=255)
                                                    tGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                                                                                   <  1:48217/150‑1 (MQ=255)
                                                       ttCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTaa                                                                                                >  1:685373/1‑150 (MQ=255)
                                                           tGACGCCATTGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCAAAAGCCTAACCAACTTGCGAAGCAGCGTATCGACACGCTCCAGCTTGATAATTTCCTCAACAACGGATTTCTAGACATCATCCCACTAAGTCa                                                                                            >  2:676642/1‑150 (MQ=255)
                                                           tGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCa                                                                                            >  1:1606478/1‑150 (MQ=255)
                                                              cGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGcc                                                                                         >  1:741562/1‑150 (MQ=255)
                                                              cGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGcc                                                                                         >  1:1380530/1‑150 (MQ=255)
                                                              cGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGcc                                                                                         <  2:486741/150‑1 (MQ=255)
                                                                ccATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCg                                                                                        >  1:1488670/1‑149 (MQ=255)
                                                                ccATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCg                                                                                        <  2:1488670/149‑1 (MQ=255)
                                                                 cATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGTTGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCt                                                                                      >  1:80728/1‑150 (MQ=255)
                                                                  aTGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTg                                                                                     <  1:776197/150‑1 (MQ=255)
                                                                  aTGGAAGATGTCCAGACGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTg                                                                                     <  1:803417/150‑1 (MQ=255)
                                                                   tGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTgg                                                                                    >  2:336480/1‑150 (MQ=255)
                                                                    ggAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGa                                                                                   >  1:159180/1‑150 (MQ=255)
                                                                     gAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAg                                                                                  <  1:43958/150‑1 (MQ=255)
                                                                      aaGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCCCTAAGTCAGCCGCTGGAGt                                                                                 >  1:403840/1‑150 (MQ=255)
                                                                           gTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcgg                                                                            >  1:1758/1‑150 (MQ=255)
                                                                            tCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcggc                                                                           >  1:712251/1‑150 (MQ=255)
                                                                             ccAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcggcg                                                                          <  1:1329457/150‑1 (MQ=255)
                                                                              cAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcggcgg                                                                         <  2:1209495/150‑1 (MQ=255)
                                                                              cAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcggcgg                                                                         <  1:52188/150‑1 (MQ=255)
                                                                               aGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCggg                                                                        >  2:638252/1‑150 (MQ=255)
                                                                                gCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGa                                                                       >  2:801074/1‑150 (MQ=255)
                                                                                   gCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAGATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAgg                                                                    <  1:300055/150‑1 (MQ=255)
                                                                                      gAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGtt                                                                 <  1:1588868/150‑1 (MQ=255)
                                                                                        aaCGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGc                                                               <  1:1011081/150‑1 (MQ=255)
                                                                                          cGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCaa                                                             >  2:1100005/1‑150 (MQ=255)
                                                                                                 ccGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgggg                                                      >  1:1416505/1‑150 (MQ=255)
                                                                                                 ccGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgggg                                                      <  1:125881/150‑1 (MQ=255)
                                                                                                  cGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGt                                                     <  2:1007363/150‑1 (MQ=255)
                                                                                                      gTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTg                                                 <  1:1463430/150‑1 (MQ=255)
                                                                                                      gTCCACTGGCATCCAACGCGTAACCAACGGGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTg                                                 <  1:1100005/150‑1 (MQ=255)
                                                                                                       tCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTgg                                                >  1:201921/1‑150 (MQ=255)
                                                                                                         cACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTTGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGAt                                              <  1:676642/150‑1 (MQ=255)
                                                                                                             ggCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGc                                          >  1:320861/1‑150 (MQ=255)
                                                                                                              gCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCt                                         >  1:202957/1‑150 (MQ=255)
                                                                                                              gCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCt                                         <  2:646201/150‑1 (MQ=255)
                                                                                                               cATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTg                                        >  2:1070661/1‑150 (MQ=255)
                                                                                                                aTCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgc                                       <  1:801074/150‑1 (MQ=255)
                                                                                                                 tCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcg                                      <  2:110172/150‑1 (MQ=255)
                                                                                                                  ccAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcgc                                     <  1:740994/150‑1 (MQ=255)
                                                                                                                    aaCGCGAAACCAAAGTGCGCCGCAGCGAAACGACACGCTGCCGCTTGATAATTTCCTCACCACCGCATTTCTCGACATCATCCCACTAATACAGCAGCTCGACTTTAGGCCGGAACTGATGAAACATTTGCTGCAGCATAAGCTGCGCat                                   >  2:557867/1‑149 (MQ=255)
                                                                                                                        cgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGt                               <  2:1380530/150‑1 (MQ=255)
                                                                                                                        cgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGt                               >  1:1247580/1‑150 (MQ=255)
                                                                                                                            aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAat                           >  1:791693/1‑150 (MQ=255)
                                                                                                                             ccAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGTAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAata                          >  2:324323/1‑150 (MQ=255)
                                                                                                                             ccAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAata                          >  1:1006441/1‑150 (MQ=255)
                                                                                                                             ccAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAata                          <  2:216891/150‑1 (MQ=255)
                                                                                                                              cAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTCGAGTTTCGTCGCGAAGGGTTGCAACATCTGGTGCTGGATAAGCTGCGCAGTGGTAAataa                         >  2:185619/1‑149 (MQ=255)
                                                                                                                               aaCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGT‑GTAAATATcc                       >  1:950861/1‑150 (MQ=255)
                                                                                                                                aCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATcc                       <  2:1069463/150‑1 (MQ=255)
                                                                                                                                   tGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCa                    <  1:364457/150‑1 (MQ=255)
                                                                                                                                    gcgcCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTCCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCaa                   >  2:633949/1‑150 (MQ=255)
                                                                                                                                    aagctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCaa                   >  1:185597/6‑150 (MQ=255)
                                                                                                                                      gctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCCCCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACa                 >  1:1109720/4‑150 (MQ=255)
                                                                                                                                      gccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCAGCGGATTTCTCGACATCATGCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCGGGATAAGCTGCGCAGTGGTAAATATCCGCAAcg                 >  2:1590073/1‑149 (MQ=255)
                                                                                                                                      gccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACa                 >  1:559399/1‑150 (MQ=255)
                                                                                                                                       ctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCGGTGGTAAATATCCGCAACAg                >  1:721946/3‑150 (MQ=255)
                                                                                                                                            gCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAg           <  2:1606478/150‑1 (MQ=255)
                                                                                                                                             cGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGc          <  2:438150/150‑1 (MQ=255)
                                                                                                                                               tATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCt        >  1:81715/1‑150 (MQ=255)
                                                                                                                                                aTCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTg       >  1:359672/1‑150 (MQ=255)
                                                                                                                                                aTCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTg       >  1:1601389/1‑150 (MQ=255)
                                                                                                                                                aTCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTg       >  2:1099420/1‑150 (MQ=255)
                                                                                                                                                     cacGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTCCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCTCGCCTgttgcg  <  1:870218/150‑6 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCT  >  NC_000913/1405949‑1406247

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