Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 738,853 G→A *254* (TAG→TAA ybfD → H repeat‑associated putative transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913738,8530GA98.9% 284.9 / NA 87*254* (TAG→TAAybfDH repeat‑associated putative transposase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (40/46);  minor base T (0/1);  total (40/47)
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

TTTCCGCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAGAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAAC  >  NC_000913/738719‑738985
                                                                                                                                      |                                                                                                                                    
tttCCGCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAg                                                                                                                                  <  1:1106267/139‑1 (MQ=255)
tttCCGCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAg                                                                                                                                  <  2:1516109/139‑1 (MQ=255)
        gAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTa                                                                                                                          <  2:718328/139‑1 (MQ=255)
             aaTTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGtt                                                                                                                     <  2:309963/139‑1 (MQ=255)
                tAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGa                                                                                                                  <  2:801508/139‑1 (MQ=255)
                       tCCCAAGCATGAGAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTGACGTTTGATTAAAAAGGGCAGCTATGGTTAGAAACTAcc                                                                                                           >  2:1039463/1‑139 (MQ=255)
                         ccAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTg                                                                                                         >  2:929180/1‑139 (MQ=255)
                            aGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGtt                                                                                                                     <  2:537516/124‑1 (MQ=255)
                            aGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGtt                                                                                                                     >  1:537516/1‑124 (MQ=255)
                            aGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACg                                                                                                      >  2:1584679/1‑139 (MQ=255)
                                  aCAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTc                                                                                                >  1:1022475/1‑139 (MQ=255)
                                  aCAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTc                                                                                                >  1:196074/1‑139 (MQ=255)
                                   catttACGCAATCATTTAAAAGAGCTGCGGCAGACGAGCATCCCGTCTTCATATTGTTTGCGATGTCCCTGATTAGTTTATTGAATTCACGTTTGAATATAAAGGGCAACTATGGTTAGAAGCTACCTGACGTCAGTcc                                                                                               <  1:1039463/136‑1 (MQ=18)
                                        aCGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCg                                                                                          <  1:929180/139‑1 (MQ=255)
                                             aTTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAg                                                                                     <  1:150347/139‑1 (MQ=255)
                                              ttAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGc                                                                                    >  2:784092/1‑139 (MQ=255)
                                                   gAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCt                                                                               <  1:328162/139‑1 (MQ=255)
                                                            cACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAggg                                                                                                                                >  1:1568346/1‑81 (MQ=255)
                                                            cACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAggg                                                                                                                                <  2:290634/81‑1 (MQ=255)
                                                            cACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAggg                                                                                                                                <  2:1568346/81‑1 (MQ=255)
                                                            cACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATAGATTTCACGTTTGAATAAAAAggg                                                                                                                                >  1:290634/1‑81 (MQ=255)
                                                              cGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAAttt                                                                    >  2:1324398/1‑139 (MQ=255)
                                                                   gagaAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGTAGGCTTTCGTAAATTTGTCCt                                                               >  2:1553575/1‑139 (MQ=255)
                                                                    agaagaAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTg                                                              >  2:1667953/1‑139 (MQ=255)
                                                                     gaagaaACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGc                                                             >  1:1287254/1‑139 (MQ=255)
                                                                     gaagaaACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGc                                                             >  1:1599987/1‑139 (MQ=255)
                                                                         aaaCCCGTCTTTATATCGTTTGCGATGGCCCTGATGAGCTTATTGATTACACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGACCTTGCGTGGAGCAGGCTTTCGTAAATTTGTCCTGCTaca                                                         <  1:1596076/139‑1 (MQ=255)
                                                                         aaaCCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTaca                                                         >  2:1163755/1‑139 (MQ=255)
                                                                          aaCCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTacaa                                                        <  1:1123690/139‑1 (MQ=255)
                                                                              cGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAggg                                                                                                                                >  2:820878/1‑63 (MQ=255)
                                                                              cGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAggg                                                                                                                                <  1:820878/63‑1 (MQ=255)
                                                                                tCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGAt                                                  >  2:1596076/1‑139 (MQ=255)
                                                                                 cTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGAtt                                                 <  2:561017/139‑1 (MQ=255)
                                                                                  ttCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTa                                                >  2:219629/1‑139 (MQ=255)
                                                                                   tCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTaa                                               <  1:1163755/139‑1 (MQ=255)
                                                                                           tttGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACaa                                       >  1:123369/1‑139 (MQ=255)
                                                                                           tttGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACaa                                       >  1:600301/1‑139 (MQ=255)
                                                                                            ttGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTc                                                                 >  2:143995/1‑112 (MQ=255)
                                                                                            ttGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTc                                                                 <  1:143995/112‑1 (MQ=255)
                                                                                            ttGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACaaa                                      <  2:1022475/139‑1 (MQ=255)
                                                                                            ttGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACaaa                                      <  2:1599987/139‑1 (MQ=255)
                                                                                            ttGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACaaa                                      <  2:1287254/139‑1 (MQ=255)
                                                                                             tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAAtt                                                                     >  1:818664/1‑107 (MQ=255)
                                                                                             tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAAtt                                                                     <  2:818664/107‑1 (MQ=255)
                                                                                             tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTcc                                                                >  2:1671284/1‑112 (MQ=255)
                                                                                             tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTcc                                                                <  1:1671284/112‑1 (MQ=255)
                                                                                             tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACaaa                                      <  1:417840/138‑1 (MQ=255)
                                                                                             tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACaaa                                      >  2:417840/1‑138 (MQ=255)
                                                                                              gCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTacaa                                                        >  2:1855930/1‑119 (MQ=255)
                                                                                              gCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTacaa                                                        <  1:1855930/119‑1 (MQ=255)
                                                                                                  tGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCAt                                <  1:1324398/139‑1 (MQ=255)
                                                                                                     cccTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                                  <  1:308446/102‑1 (MQ=255)
                                                                                                     cccTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                                  >  2:308446/1‑102 (MQ=255)
                                                                                                      ccTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACGGATTAACTTCACAAATATCATTtctc                           <  2:196074/139‑1 (MQ=255)
                                                                                                       cTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGTAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtctc                           <  1:1553575/139‑1 (MQ=255)
                                                                                                       cTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtctc                           <  1:784092/139‑1 (MQ=255)
                                                                                                       cTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtctc                           <  2:600301/139‑1 (MQ=255)
                                                                                                       cTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtctc                           <  1:615423/139‑1 (MQ=255)
                                                                                                             aGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTACCTTCACAAATATCATTTCTCAACGTc                     <  2:575688/139‑1 (MQ=255)
                                                                                                              gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTTAATTTGt                                                                  <  1:1230926/93‑1 (MQ=255)
                                                                                                              gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                                  >  2:1230926/1‑93 (MQ=255)
                                                                                                              gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTaca                                                         >  2:24415/1‑102 (MQ=255)
                                                                                                              gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTaca                                                         <  1:24415/102‑1 (MQ=255)
                                                                                                              gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTcaca                                        <  1:1377642/119‑1 (MQ=255)
                                                                                                              gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTcaca                                        >  2:1377642/1‑119 (MQ=255)
                                                                                                              gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtctc                           <  2:1846898/132‑1 (MQ=255)
                                                                                                              gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtctc                           >  1:677606/1‑132 (MQ=255)
                                                                                                              gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtctc                           <  2:677606/132‑1 (MQ=255)
                                                                                                              gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtctc                           >  1:1846898/1‑132 (MQ=255)
                                                                                                               cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTAcc                                                                                                           >  1:514896/1‑51 (MQ=255)
                                                                                                               cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTAcc                                                                                                           <  2:514896/51‑1 (MQ=255)
                                                                                                               cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCt                                                                                                          <  2:1156989/52‑1 (MQ=255)
                                                                                                               cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCt                                                                                                          >  1:1156989/1‑52 (MQ=255)
                                                                                                               cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTa                   >  1:262102/1‑139 (MQ=255)
                                                                                                               cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTa                   >  1:1237434/1‑139 (MQ=255)
                                                                                                                  aTTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTc                                                                 <  2:387658/90‑1 (MQ=255)
                                                                                                                  aTTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTc                                                                 >  1:387658/1‑90 (MQ=255)
                                                                                                                     gATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTa             <  2:1347275/139‑1 (MQ=255)
                                                                                                                         tCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTaca                                                         <  2:85596/91‑1 (MQ=255)
                                                                                                                         tCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTaca                                                         >  1:85596/1‑91 (MQ=255)
                                                                                                                         tCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTcaca                                        >  2:1288733/1‑108 (MQ=255)
                                                                                                                         tCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTcaca                                        <  1:1288733/108‑1 (MQ=255)
                                                                                                                               ttGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTaaa   <  1:1667953/139‑1 (MQ=255)
                                                                                                                               ttGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTaaa   <  1:1873480/139‑1 (MQ=255)
                                                                                                                                tGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTacaa                                                        >  2:125684/1‑85 (MQ=255)
                                                                                                                                tGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTacaa                                                        <  1:125684/85‑1 (MQ=255)
                                                                                                                                tGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAAc  >  2:477647/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                    
TTTCCGCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAGAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAAC  >  NC_000913/738719‑738985

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: