Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,786,746 Δ1 bp pseudogene (350/2253 nt) ygaQ → uncharacterized protein, putative enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,786,7460C.98.8% 336.1 / ‑3.1 82pseudogene (350/2253 nt)ygaQuncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (39/42);  minor base T (1/0);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACATTCCTGACTATGAG  >  NC_000913/2786615‑2786879
                                                                                                                                   |                                                                                                                                     
gTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTCGAATAAAGCAGCACAGCGTAGCATTACAGTGCCAGACAGTATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTg                                                                                                                               >  2:66308/1‑139 (MQ=255)
 ttAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGa                                                                                                                              >  2:107636/1‑139 (MQ=255)
   aaaaaaCGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTATCAACGAGATAAATGATTGa                                                                                                                              >  1:1022089/1‑136 (MQ=255)
   aaaaaaCGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTATCAACGAGATAAATGATTGa                                                                                                                              <  2:1022089/136‑1 (MQ=255)
          gATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGa                                                                                                                     <  1:327152/139‑1 (MQ=255)
           aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGAt                                                                                                                    >  1:730132/1‑139 (MQ=255)
           aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGAt                                                                                                                    >  2:1001134/1‑139 (MQ=255)
                  aaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCg                                                                                                             <  1:1001134/139‑1 (MQ=255)
                  aaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTAAAAGAGGTGATAACAgag                                                                                                             <  2:730132/139‑3 (MQ=255)
                        ccTTTCTCGAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                                       <  2:802062/139‑1 (MQ=255)
                        ccTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                                       <  1:193306/139‑1 (MQ=255)
                                       aaCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTg                                                                                        >  2:307372/1‑139 (MQ=255)
                                        aCTTGCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGc                                                                                       >  2:764000/1‑139 (MQ=255)
                                        aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGc                                                                                       >  2:239130/1‑139 (MQ=255)
                                        aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGc                                                                                       >  2:315522/1‑139 (MQ=255)
                                          ttACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTg                                                                                     <  2:26614/139‑1 (MQ=255)
                                            aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGaa                                                                                   >  2:569618/1‑139 (MQ=255)
                                             ccGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAAc                                                                                  >  2:711715/1‑139 (MQ=255)
                                             ccGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGACTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAAc                                                                                  >  2:827039/1‑139 (MQ=255)
                                              cGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACa                                                                                 >  2:1004460/1‑139 (MQ=255)
                                                 gagTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGaaa                                                                                                                            <  2:75279/93‑1 (MQ=255)
                                                 gagTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGaaa                                                                                                                            >  1:75279/1‑93 (MQ=255)
                                                 gagTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACAttt                                                                              <  1:1194162/139‑1 (MQ=255)
                                                 gagTTGAATAAAGCAGCAAAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACAttt                                                                              >  2:729620/1‑139 (MQ=255)
                                                    ttGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                                                                           <  2:761443/139‑1 (MQ=255)
                                                     tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                          <  1:569618/139‑1 (MQ=255)
                                                         tAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                                      <  1:711715/139‑1 (MQ=255)
                                                         tAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGCATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGACTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                                      <  1:827039/139‑1 (MQ=255)
                                                            agcagcACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATgg                                                                                             >  1:610298/1‑113 (MQ=255)
                                                            agcagcACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATgg                                                                                             <  2:610298/113‑1 (MQ=255)
                                                              cagcaCAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtata                                                                 <  1:942771/139‑1 (MQ=255)
                                                                 cacaGCGTAGCATTACAGCGCCAGACAGGCTAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAAAGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTTCTGAACATTTCACCCACATTATATCt                                                              >  2:188812/1‑139 (MQ=255)
                                                                      cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCgg                                                                                                            >  1:392101/1‑88 (MQ=255)
                                                                      cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCgg                                                                                                            <  2:392101/88‑1 (MQ=255)
                                                                      cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaa                                                         >  2:1248686/1‑139 (MQ=255)
                                                                                cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                          >  1:984535/1‑112 (MQ=255)
                                                                                cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                          <  2:984535/112‑1 (MQ=255)
                                                                                cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCa                                               <  2:130373/139‑1 (MQ=255)
                                                                                 aGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaa                                              <  1:764000/139‑1 (MQ=255)
                                                                                  gcgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaaa                                             >  1:1296571/1‑139 (MQ=255)
                                                                                   cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                               >  1:862989/1‑72 (MQ=255)
                                                                                   cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                               <  2:862989/72‑1 (MQ=255)
                                                                                   cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCt                                                     >  1:574950/1‑130 (MQ=255)
                                                                                   cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCt                                                     <  2:574950/130‑1 (MQ=255)
                                                                                    gcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGAc                                                           >  2:1209729/1‑123 (MQ=255)
                                                                                    gcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGAc                                                           <  1:1209729/123‑1 (MQ=255)
                                                                                     ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTa                                                                                                >  1:1022843/1‑85 (MQ=255)
                                                                                     ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTa                                                                                                <  2:1022843/85‑1 (MQ=255)
                                                                                     ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTg                                                                                        <  2:1224301/93‑1 (MQ=255)
                                                                                     ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTg                                                                                        >  1:1224301/1‑93 (MQ=255)
                                                                                        gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAccc                                                                         <  2:1177649/105‑1 (MQ=255)
                                                                                        gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAccc                                                                         >  1:1177649/1‑105 (MQ=255)
                                                                                        gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATaa                                       >  1:1208313/1‑139 (MQ=255)
                                                                                        gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATaa                                       <  1:1004460/139‑1 (MQ=255)
                                                                                            ggATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAatat                                   <  2:1063931/139‑1 (MQ=255)
                                                                                             gATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATg                                  <  1:188812/139‑1 (MQ=255)
                                                                                                  gagaAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                                                                           <  2:1151866/93‑1 (MQ=255)
                                                                                                  gagaAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                                                                           >  1:1151866/1‑93 (MQ=255)
                                                                                                   agaaagaaAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCat                            <  1:417048/139‑1 (MQ=255)
                                                                                                    gaaagaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata                           <  1:107636/139‑1 (MQ=255)
                                                                                                     aaagaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCatat                          <  1:160761/139‑1 (MQ=255)
                                                                                                       agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                               >  1:56237/1‑52 (MQ=255)
                                                                                                       agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                               <  2:56237/52‑1 (MQ=255)
                                                                                                       agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgt                                                                                                           >  1:532562/1‑56 (MQ=255)
                                                                                                       agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgt                                                                                                           <  2:532562/56‑1 (MQ=255)
                                                                                                       agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                          <  1:835607/89‑1 (MQ=255)
                                                                                                       agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                          >  2:835607/1‑89 (MQ=255)
                                                                                                       agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata                           <  1:489407/136‑1 (MQ=255)
                                                                                                       agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata                           >  2:489407/1‑136 (MQ=255)
                                                                                                       agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata                           >  1:178169/1‑136 (MQ=255)
                                                                                                       agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata                           <  2:178169/136‑1 (MQ=255)
                                                                                                       agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc                        >  2:1066028/1‑139 (MQ=255)
                                                                                                       agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc                        >  1:585179/1‑139 (MQ=255)
                                                                                                       agaaaAACATAAGTTTATCAACGAGATTAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAGCATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATGTGTTGCATATGc                        >  2:783749/1‑139 (MQ=255)
                                                                                                        gaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCt                       <  1:307372/139‑1 (MQ=255)
                                                                                                              cATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAAc                 <  1:1066028/139‑1 (MQ=255)
                                                                                                                taagTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAAc                 <  2:1034217/134‑1 (MQ=255)
                                                                                                                taagTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAAc                 >  1:1034217/4‑137 (MQ=255)
                                                                                                                   gTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAg                                           <  2:1138250/108‑1 (MQ=255)
                                                                                                                   gTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAg                                           >  1:1138250/1‑108 (MQ=255)
                                                                                                                        tCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACATTCCTGACt       <  1:451163/139‑1 (MQ=255)
                                                                                                                             gagaTAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACATTCCTGACTATGAg  >  1:156269/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                     
GTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACATTCCTGACTATGAG  >  NC_000913/2786615‑2786879

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