Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,014,287 C→T *257* (TAG→TAA yqeC ← putative selenium‑dependent hydroxylase accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,014,2870CT98.9% 287.1 / NA 89*257* (TAG→TAAyqeCputative selenium‑dependent hydroxylase accessory protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (44/44);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACG  >  NC_000913/3014153‑3014418
                                                                                                                                      |                                                                                                                                   
aGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAAtt                                                                                                                                >  2:782364/1‑139 (MQ=255)
     aGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACaa                                                                                                                           >  2:394160/1‑139 (MQ=255)
              tATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                  <  2:383376/139‑1 (MQ=255)
              tATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                  <  1:1696040/139‑1 (MQ=255)
                   acaaTTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGAtt                                                                                                             <  1:394160/139‑1 (MQ=255)
                       tttGCGCCATAGCACACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGCGAAAATATTCATACCAACCTCAAAGGATATCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAg                                                                                                         >  2:72867/1‑139 (MQ=255)
                       tttGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAg                                                                                                         >  1:149663/1‑139 (MQ=255)
                        ttGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                                        <  2:149663/139‑1 (MQ=255)
                        ataagagacaGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTCAAAATATTCATACCAACCTCAATTGATATCTTTCAGTAATCCGGAATACCTGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                                        <  2:1196142/130‑1 (MQ=255)
                         tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                       <  1:997196/139‑1 (MQ=255)
                         tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                       <  2:154307/139‑1 (MQ=255)
                           cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAGTCTGCGCGCGAtt                                                                                                             >  1:1507774/1‑131 (MQ=255)
                           cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAGTCTGCGCGCGAtt                                                                                                             <  2:1507774/131‑1 (MQ=255)
                           cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                     <  2:1529088/139‑1 (MQ=255)
                                 aGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                                                    <  2:786901/118‑1 (MQ=255)
                                 aGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                                                    >  1:786901/1‑118 (MQ=255)
                                    aaaCGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGtat                                                                                            <  1:782364/139‑1 (MQ=255)
                                      aCGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTatat                                                                                          >  2:362337/1‑139 (MQ=255)
                                        ggACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCg                                                                                        >  2:920195/1‑139 (MQ=255)
                                        ggACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCg                                                                                        <  1:1837266/139‑1 (MQ=255)
                                          aCAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATCCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGccc                                                                                     >  2:374066/1‑139 (MQ=255)
                                              ttttGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCGGCGCGCGATTTCAGGATGTTCTTGtactgtctcttat                                                                                  >  1:1196142/1‑128 (MQ=255)
                                                   ctctTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAt                                                                             <  1:362337/139‑1 (MQ=255)
                                                        cGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGtat                                                                                            >  2:291533/1‑119 (MQ=255)
                                                        cGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGtat                                                                                            <  1:291533/119‑1 (MQ=255)
                                                         gAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATCCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCctct                                                                      <  1:374066/139‑1 (MQ=255)
                                                           gTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTa                                                                     <  1:466330/139‑1 (MQ=255)
                                                               tGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtc                                                                 >  2:1401063/1‑139 (MQ=255)
                                                                                  ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                     >  2:1706086/1‑84 (MQ=255)
                                                                                  ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                     <  1:1706086/84‑1 (MQ=255)
                                                                                  ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGa                                              <  1:920195/139‑1 (MQ=255)
                                                                                    cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTa                                                                     <  1:826871/114‑1 (MQ=255)
                                                                                    cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTa                                                                     >  2:826871/1‑114 (MQ=255)
                                                                                    cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGttgttg                                                        >  2:202098/1‑127 (MQ=255)
                                                                                    cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGttgttg                                                        <  1:202098/127‑1 (MQ=255)
                                                                                    cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAg                                            >  2:1642402/1‑139 (MQ=255)
                                                                                    cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAg                                            >  2:656164/1‑139 (MQ=255)
                                                                                    cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGACTGCAGGa                                                                                                       <  2:406808/80‑1 (MQ=255)
                                                                                    cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGACTGCAGGa                                                                                                       >  1:406808/1‑80 (MQ=255)
                                                                                           aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCa                                                                                                          >  1:76915/1‑70 (MQ=255)
                                                                                           aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCa                                                                                                          <  2:76915/70‑1 (MQ=255)
                                                                                           aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                                        >  1:144509/1‑72 (MQ=255)
                                                                                           aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                                        <  2:144509/72‑1 (MQ=255)
                                                                                           aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtcgt                                                               <  2:454382/113‑1 (MQ=255)
                                                                                           aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtcgt                                                               >  1:454382/1‑113 (MQ=255)
                                                                                           aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                                                      >  2:1816368/1‑122 (MQ=255)
                                                                                           aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                                                      <  1:1816368/122‑1 (MQ=255)
                                                                                           aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg                                     >  2:313291/1‑139 (MQ=255)
                                                                                            ccTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                                        >  2:923682/1‑71 (MQ=255)
                                                                                            ccTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                                        <  1:923682/71‑1 (MQ=255)
                                                                                                 aTGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTtgt                               >  1:392441/1‑139 (MQ=255)
                                                                                                       tatcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAAt                          <  1:313291/136‑1 (MQ=255)
                                                                                                        atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATc                         >  2:1181833/3‑139 (MQ=255)
                                                                                                         tcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg                        <  2:392441/138‑1 (MQ=255)
                                                                                                             tCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAAt                          <  2:1222202/133‑1 (MQ=255)
                                                                                                             tCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAAt                          >  1:1222202/1‑133 (MQ=255)
                                                                                                             tCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCAtt                    <  1:1261090/139‑1 (MQ=255)
                                                                                                                   aCCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTGCGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTcaca              >  2:1000607/1‑139 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                  <  2:1564532/38‑1 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                  >  1:1564532/1‑38 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                  <  2:533397/38‑1 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                  >  1:533397/1‑38 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                       >  1:510322/1‑49 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                       <  2:510322/49‑1 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCtt                                                                                                <  1:416665/56‑1 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCtt                                                                                                >  2:416665/1‑56 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAt                                                                             >  2:1780731/1‑75 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAt                                                                             <  2:765797/75‑1 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAt                                                                             <  1:1780731/75‑1 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAt                                                                             >  1:765797/1‑75 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGcc                                                                         >  1:1318703/1‑79 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGcc                                                                         <  2:1318703/79‑1 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg                        >  2:1354298/1‑128 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg                        <  1:1354298/128‑1 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACAt             >  1:289455/1‑139 (MQ=255)
                                                                                                                    cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACAt             >  1:228847/1‑139 (MQ=255)
                                                                                                                         aaTACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCtt                                  <  2:1044403/113‑1 (MQ=255)
                                                                                                                         aaTACCCGGGCCGTTAATACACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCtt                                  >  1:1044403/1‑113 (MQ=255)
                                                                                                                           tACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                                        >  1:587651/1‑41 (MQ=255)
                                                                                                                           tACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                                        <  2:587651/41‑1 (MQ=255)
                                                                                                                              ccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtcg                                                                >  2:495818/1‑78 (MQ=255)
                                                                                                                              ccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtcg                                                                <  1:495818/78‑1 (MQ=255)
                                                                                                                               cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACg  >  1:1568360/1‑139 (MQ=255)
                                                                                                                                  gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                     >  1:1774439/1‑37 (MQ=255)
                                                                                                                                  gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                     <  2:1774439/37‑1 (MQ=255)
                                                                                                                                  gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGcc                                                                         <  1:9330/65‑1 (MQ=255)
                                                                                                                                  gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGcc                                                                         >  2:9330/1‑65 (MQ=255)
                                                                                                                                  gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTtgtg                              <  1:1765313/108‑1 (MQ=255)
                                                                                                                                  gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTtgtg                              >  2:1765313/1‑108 (MQ=255)
                                                                                                                                      |                                                                                                                                   
AGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACG  >  NC_000913/3014153‑3014418

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