Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,574,259 A→G pseudogene (341/2937 nt) yjiV → pseudogene, conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,574,2590AG98.6% 274.5 / NA 73pseudogene (341/2937 nt)yjiVpseudogene, conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (40/32);  total (41/32)
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

TGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAG  >  NC_000913/4574140‑4574393
                                                                                                                       |                                                                                                                                      
tgtgAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGaa                                                                                                                     >  1:1005068/1‑139 (MQ=255)
tgtgAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGaa                                                                                                                     >  2:1013355/1‑139 (MQ=255)
 gtgAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAAt                                                                                                                    >  1:1092621/1‑139 (MQ=255)
 gtgAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAAt                                                                                                                    <  2:1092621/139‑1 (MQ=255)
 gtgAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAAt                                                                                                                    <  2:1004508/139‑1 (MQ=255)
    aaCCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGATGTTGCATACTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGTGAAAATAAAGGTGGCGATGAATtct                                                                                                                 >  1:568448/1‑139 (MQ=255)
    aaCCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATtct                                                                                                                 >  2:892693/1‑139 (MQ=255)
    aaCCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATtct                                                                                                                 >  2:765725/1‑139 (MQ=255)
      ccAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTctct                                                                                                               >  2:190003/1‑139 (MQ=255)
      ccAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTctct                                                                                                               >  1:334564/1‑139 (MQ=255)
      ccAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTctct                                                                                                               >  1:982089/1‑139 (MQ=255)
      ccAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTctct                                                                                                               >  1:891954/1‑139 (MQ=255)
      ccAGCGAAAAAATGGTGAGTGTGGAAGTCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATAAGCCTCATCTTCTGGCACGGGACAAAAACGGCAAACTAGTCCAGCGGGGAAATACAGGTGGCGATGAATTctct                                                                                                               >  1:781269/1‑139 (MQ=255)
        aGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCttt                                                                                                             >  1:918675/1‑139 (MQ=255)
        aGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCttt                                                                                                             <  1:765725/139‑1 (MQ=255)
         gCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTg                                                                                                            <  2:95663/139‑1 (MQ=255)
          cGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGc                                                                                                           >  2:601825/1‑139 (MQ=255)
              aaaaTTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGa                                                                                                       <  2:982089/139‑1 (MQ=255)
              aaaaTTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGa                                                                                                       <  1:1052258/139‑1 (MQ=255)
                       aCTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCa                                                                                                          <  1:700948/127‑1 (MQ=255)
                       aCTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCa                                                                                                          >  2:700948/1‑127 (MQ=255)
                         tgtgGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTg                                                                                            <  1:1020570/139‑1 (MQ=255)
                                     gTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAg                                                                                                                                >  1:69548/1‑91 (MQ=255)
                                     gTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAg                                                                                                                                <  2:69548/91‑1 (MQ=255)
                                     gTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACg                                                                                >  1:305022/1‑139 (MQ=255)
                                     gTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACg                                                                                <  1:727564/139‑1 (MQ=255)
                                      tataCGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCa                                                                                                          <  1:1031383/112‑1 (MQ=255)
                                      tataCGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCa                                                                                                          >  2:1031383/1‑112 (MQ=255)
                                        taCGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTa                                                                             >  2:207678/1‑139 (MQ=255)
                                        taCGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTa                                                                             <  2:333442/139‑1 (MQ=255)
                                             aGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACagaa                                                                        >  1:430822/1‑139 (MQ=255)
                                              gTTGCTTTCTGTAATGAATGTAGTGAGCCTGATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTTGGGAAATAAAGGAGGCGAAGAATTCTATTCGCAGGATGAAGTACCTATTGAACATGACACTACagaag                                                                       >  1:838152/1‑139 (MQ=255)
                                              gTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACagaag                                                                       >  1:868348/1‑139 (MQ=255)
                                                tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATtct                                                                                                                 <  2:933765/95‑1 (MQ=255)
                                                tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATtct                                                                                                                 >  1:933765/1‑95 (MQ=255)
                                                tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAgaagaa                                                                     >  2:67234/1‑139 (MQ=255)
                                                tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAgaagaa                                                                     >  1:1011272/1‑139 (MQ=255)
                                                tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTGGAACATGACGCTACAgaagaa                                                                     >  2:653347/1‑139 (MQ=255)
                                                 gCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTAc                                                                            <  1:364981/131‑1 (MQ=255)
                                                 gCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTAc                                                                            >  2:364981/1‑131 (MQ=255)
                                                  cATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGaaaa                                                                   >  1:445080/1‑139 (MQ=255)
                                                            tGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAAt                                                                                                                    <  2:496523/80‑1 (MQ=255)
                                                            tGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAAt                                                                                                                    >  1:496523/1‑80 (MQ=255)
                                                                tGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAgaga                                                     <  1:302263/139‑1 (MQ=255)
                                                                 gTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTg                                                                                            >  1:65810/1‑99 (MQ=255)
                                                                 gTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTg                                                                                            <  2:65810/99‑1 (MQ=255)
                                                                   aaTGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACa                                                  >  2:173515/1‑139 (MQ=255)
                                                                    aTGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGaa                                                                                       <  1:29639/101‑1 (MQ=255)
                                                                    aTGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGaa                                                                                       >  2:29639/1‑101 (MQ=255)
                                                                        gCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAg                                                                                                  <  1:555233/86‑1 (MQ=255)
                                                                        gCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAg                                                                                                  >  2:555233/1‑86 (MQ=255)
                                                                        gCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGtttt                                             >  2:649180/1‑139 (MQ=255)
                                                                          cTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGc                                                                                                                           >  2:927435/1‑59 (MQ=255)
                                                                          cTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGc                                                                                                                           <  1:927435/59‑1 (MQ=255)
                                                                            cATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGc                                         <  1:190003/139‑1 (MQ=255)
                                                                              tcttctGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCt                                       <  2:164900/139‑1 (MQ=255)
                                                                              tcttctGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCt                                       <  2:918675/139‑1 (MQ=255)
                                                                                 tctGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCg                                    <  2:339430/139‑1 (MQ=255)
                                                                                       aCGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGatt                              >  2:370963/1‑139 (MQ=255)
                                                                                         gggACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGattat                            >  1:477716/1‑139 (MQ=255)
                                                                                         gggACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGattat                            >  2:100003/1‑139 (MQ=255)
                                                                                          ggACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGattatt                           >  1:807181/1‑139 (MQ=255)
                                                                                             cAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCg                                    >  2:735827/1‑127 (MQ=255)
                                                                                             cAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCg                                    <  1:735827/127‑1 (MQ=255)
                                                                                                        aaCTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGAc             <  1:67234/139‑1 (MQ=255)
                                                                                                              tCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGa       >  2:778052/1‑139 (MQ=255)
                                                                                                               ctaGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTACCGCAGGAGAGACCAGCGAg      <  1:653347/137‑1 (MQ=255)
                                                                                                               ccAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAg      <  2:445080/139‑1 (MQ=255)
                                                                                                                cAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAgg     <  1:398842/139‑1 (MQ=255)
                                                                                                                cAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGGTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAgg     <  1:778052/139‑1 (MQ=255)
                                                                                                                  gTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCa   <  2:477716/139‑1 (MQ=255)
                                                                                                                   tGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACg                                                                                >  2:1076643/1‑61 (MQ=255)
                                                                                                                   tGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACg                                                                                <  1:1076643/61‑1 (MQ=255)
                                                                                                                   tGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAg  >  2:267924/1‑139 (MQ=255)
                                                                                                                       |                                                                                                                                      
TGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAG  >  NC_000913/4574140‑4574393

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: