Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,519,341 A→G noncoding (69/76 nt) lysV → tRNA‑Lys

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,519,3410AG99.0% 238.0 / ‑5.1 96noncoding (69/76 nt)lysVtRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (49/46);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGC  >  NC_000913/2519251‑2519429
                                                                                          |                                                                                        
tCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGcc                                                                                        >  1:3780090/1‑93 (MQ=37)
tCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGcc                                                                                        >  2:1425081/1‑93 (MQ=37)
  ccGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCcca                                                                                      >  2:57790/1‑93 (MQ=37)
  ccGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCcca                                                                                      >  1:2923802/1‑93 (MQ=37)
  ccGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCcca                                                                                      >  2:3406027/1‑93 (MQ=37)
    gTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCccacc                                                                                    >  1:1947700/1‑93 (MQ=37)
    gTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTCCACGGCccacc                                                                                    >  1:2314333/1‑93 (MQ=255)
      catcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCaa                                                                                  >  1:4120872/1‑93 (MQ=255)
       atcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAAt                                                                                 <  2:2992175/93‑1 (MQ=255)
          acccaccACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTa                                                                              >  1:1277940/1‑93 (MQ=255)
           cccaccACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaa                                                                             <  1:4237787/93‑1 (MQ=255)
            ccaccaCTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaaa                                                                            >  2:1507675/1‑93 (MQ=255)
            ccaccaCTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaaa                                                                            <  2:3075594/93‑1 (MQ=255)
            ccaccaCTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaaa                                                                            >  2:4351680/1‑93 (MQ=255)
             caccacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaaaa                                                                           <  2:3568446/93‑1 (MQ=255)
             caccacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaaaa                                                                           >  2:4048061/1‑93 (MQ=255)
             caccacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaaaa                                                                           >  2:4048062/1‑93 (MQ=255)
             caccacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaaaa                                                                           >  1:1322208/1‑93 (MQ=255)
              accacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCAGGGCCCACCAATGTaaaaa                                                                          >  2:4152437/1‑93 (MQ=255)
               ccacTTCGGTTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaaaaaa                                                                         <  2:3415702/93‑1 (MQ=255)
               ccacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaaaaaa                                                                         <  2:2715811/93‑1 (MQ=255)
               ccacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaaaaaa                                                                         >  2:31566/1‑93 (MQ=255)
                cacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAg                                                                        <  1:1104320/93‑1 (MQ=37)
                 ccTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAgc                                                                       <  2:1426387/92‑1 (MQ=19)
                 ccTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAgc                                                                       <  2:588046/92‑1 (MQ=19)
                  cTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAgcg                                                                      >  1:491274/1‑93 (MQ=34)
                   ttCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAgcgc                                                                     >  2:1656687/1‑93 (MQ=18)
                   ttCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAgcgc                                                                     >  1:4101748/1‑93 (MQ=18)
                    tCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGcc                                                                    <  2:493325/93‑1 (MQ=17)
                    tCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGcc                                                                    <  2:4275791/93‑1 (MQ=17)
                                     ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTa                                                   <  1:1782386/93‑1 (MQ=14)
                                       ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACt                                                 <  2:1497471/93‑1 (MQ=18)
                                         tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTAt                                               >  2:2040480/1‑93 (MQ=18)
                                         tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTAt                                               >  2:786057/1‑93 (MQ=18)
                                          agagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATc                                              <  2:4425008/93‑1 (MQ=18)
                                           gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCt                                             >  1:3156294/1‑93 (MQ=18)
                                           gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCt                                             <  1:3967935/93‑1 (MQ=21)
                                             gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGc                                           >  2:3995017/1‑93 (MQ=18)
                                              cAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCg                                          >  1:1420174/1‑93 (MQ=18)
                                              cAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCg                                          >  2:593949/1‑93 (MQ=18)
                                               aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGa                                         >  1:3526060/1‑93 (MQ=18)
                                               aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGa                                         >  1:1388489/1‑93 (MQ=18)
                                               aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGa                                         >  1:908713/1‑93 (MQ=18)
                                               aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGa                                         <  2:313679/93‑1 (MQ=21)
                                               aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGa                                         >  1:279655/1‑93 (MQ=18)
                                                gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGAt                                        >  1:3199167/1‑93 (MQ=18)
                                                 ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATa                                       >  1:4125561/1‑93 (MQ=18)
                                                 ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATa                                       >  1:1603847/1‑93 (MQ=18)
                                                 ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATa                                       <  1:821956/93‑1 (MQ=18)
                                                  tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAc                                      >  1:3292572/1‑93 (MQ=18)
                                                  tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAc                                      >  2:697406/1‑93 (MQ=18)
                                                  tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAc                                      <  1:2956663/93‑1 (MQ=18)
                                                   gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACt                                     <  1:3288361/93‑1 (MQ=18)
                                                   gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACt                                     >  2:2565014/1‑93 (MQ=18)
                                                   gACTTTTAAACAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACt                                     >  1:2600601/1‑93 (MQ=255)
                                                    aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTc                                    >  2:191407/1‑93 (MQ=18)
                                                      ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaa                                  >  2:1671532/1‑93 (MQ=18)
                                                      ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaa                                  <  2:712314/93‑1 (MQ=21)
                                                       tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaaa                                 <  2:558339/93‑1 (MQ=21)
                                                       tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaaa                                 <  2:2926767/93‑1 (MQ=21)
                                                        ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAg                                >  2:3126814/1‑93 (MQ=18)
                                                         tAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGa                               <  2:2327630/93‑1 (MQ=18)
                                                         tAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGa                               <  1:4737899/93‑1 (MQ=21)
                                                         tAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGa                               <  1:610726/93‑1 (MQ=18)
                                                         tAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGa                               <  1:2775043/93‑1 (MQ=21)
                                                          aaTCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAATAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAt                              >  2:4245178/1‑93 (MQ=16)
                                                           aTCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAtt                             <  2:3790581/93‑1 (MQ=21)
                                                              aaTTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGa                          <  2:1202963/93‑1 (MQ=21)
                                                              aaTTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGa                          <  1:1795425/93‑1 (MQ=18)
                                                                ttGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAAc                        >  1:4052628/1‑93 (MQ=18)
                                                                ttGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAAc                        >  1:4052627/1‑93 (MQ=18)
                                                                 tGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAAcc                       >  1:4176655/1‑93 (MQ=18)
                                                                   gTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTg                     <  1:202826/93‑1 (MQ=18)
                                                                   gTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTg                     >  2:220018/1‑93 (MQ=18)
                                                                    tCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgc                    <  2:4106339/93‑1 (MQ=18)
                                                                      gCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcag                  <  1:3987779/93‑1 (MQ=18)
                                                                      gCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcag                  >  1:956780/1‑93 (MQ=18)
                                                                       cAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcagc                 >  2:1406315/1‑93 (MQ=18)
                                                                        aGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcagca                >  2:2496801/1‑93 (MQ=18)
                                                                         ggTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGcagcag               <  1:3492428/93‑1 (MQ=21)
                                                                         ggTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGcagcag               <  2:2572407/93‑1 (MQ=21)
                                                                         ggTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGcagcag               <  2:2138491/93‑1 (MQ=18)
                                                                           ttCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGt             >  1:3283840/1‑93 (MQ=18)
                                                                           ttCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGt             <  2:2113805/93‑1 (MQ=18)
                                                                           ttCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGt             <  2:2006057/93‑1 (MQ=18)
                                                                            tCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGtt            <  1:3816987/93‑1 (MQ=18)
                                                                            tCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGtt            <  2:1290292/93‑1 (MQ=18)
                                                                             cGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt           <  2:1724005/93‑1 (MQ=18)
                                                                               aaTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga         >  1:1294385/1‑93 (MQ=18)
                                                                                aTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag        >  2:1205429/1‑93 (MQ=18)
                                                                                 tCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt       >  2:1203786/1‑93 (MQ=18)
                                                                                   cTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttg     <  2:806163/93‑1 (MQ=18)
                                                                                   cTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttg     >  2:428465/1‑93 (MQ=18)
                                                                                    tGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga    >  2:4576859/1‑93 (MQ=18)
                                                                                    tGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga    <  2:3984166/93‑1 (MQ=18)
                                                                                    tGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga    <  2:3984165/93‑1 (MQ=21)
                                                                                    tGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga    <  1:1075953/93‑1 (MQ=18)
                                                                                    tGCACGCCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga    <  2:4536716/93‑1 (MQ=18)
                                                                                     gCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag   <  1:1892417/93‑1 (MQ=18)
                                                                                      cACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgc  <  1:3570170/93‑1 (MQ=21)
                                                                                          |                                                                                        
TCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGC  >  NC_000913/2519251‑2519429

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