Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,311,831 C→T 19.7% intergenic (+23/+17) tonB → / ← yciA membrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,311,8310CT19.7% 169.3 / 30.7 80intergenic (+23/+17)tonB/yciAmembrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base C (23/41);  new base T (12/4);  total (35/45)
Fisher's exact test for biased strand distribution p-value = 9.63e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

ATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTT  >  NC_000913/1311732‑1311930
                                                                                                   |                                                                                                   
aTGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAc                                                                                                     <  1:550084/100‑1 (MQ=255)
    gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGa                                                                                                 <  1:2318049/100‑1 (MQ=255)
     ccGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAg                                                                                                <  2:1021428/100‑1 (MQ=255)
        ggTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCt                                                                                             >  1:1865182/1‑100 (MQ=255)
        ggTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCt                                                                                             <  1:1330721/100‑1 (MQ=255)
          tAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCttt                                                                                           <  1:1208214/100‑1 (MQ=255)
           aaGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                                                                          <  1:929959/100‑1 (MQ=255)
           aaGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                          <  2:1834955/100‑1 (MQ=255)
           aaGCAAGGAAGTAGGATTGTGGTGAATATCATGTTTAAAATTACCGACACCACCGAGATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                          <  2:479842/100‑1 (MQ=255)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                         <  2:2466591/100‑1 (MQ=255)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                         <  2:305991/100‑1 (MQ=255)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                         <  2:1945192/100‑1 (MQ=255)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                         <  2:210348/100‑1 (MQ=255)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                         <  2:691786/100‑1 (MQ=255)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                         <  2:2356785/100‑1 (MQ=255)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                         <  1:988828/100‑1 (MQ=255)
              ccAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGAc                                                                                       >  1:903735/1‑100 (MQ=255)
              ccAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGAc                                                                                       <  1:1101989/100‑1 (MQ=255)
               cAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                         <  2:1224505/97‑1 (MQ=255)
                  gCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTAtt                                                                                   <  1:724965/100‑1 (MQ=255)
                      tGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTc                                                                               <  2:684840/100‑1 (MQ=255)
                      tGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTc                                                                               <  1:685175/100‑1 (MQ=255)
                         gATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAAc                                                                            <  2:2511712/100‑1 (MQ=255)
                             gtggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAAGCTTTTGACTATTACTCAACAGGt                                                                        >  2:143815/1‑100 (MQ=255)
                              tggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTa                                                                       <  1:1200690/100‑1 (MQ=255)
                              tcttgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTa                                                                       <  1:1694400/97‑1 (MQ=255)
                               ggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGCCTCATACTCAACAGGTaa                                                                      >  2:33357/1‑100 (MQ=255)
                               ggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTaa                                                                      >  1:131558/1‑100 (MQ=255)
                               ggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCACAGGATTTACACTATTACTCAACAAGTaa                                                                      >  2:666858/1‑100 (MQ=255)
                               ctcttccgatcTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTaa                                                                      <  1:1224129/90‑1 (MQ=255)
                                gtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAg                                                                     <  1:1695316/100‑1 (MQ=255)
                                 tgAATATCCTGTTTAAAATTAACGGCACCACCGCAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAga                                                                    >  2:1694388/1‑99 (MQ=255)
                                   aaTATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGgcg                                                                  <  1:1670360/100‑1 (MQ=255)
                                        ccTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAgg                                                             <  2:1970710/100‑1 (MQ=255)
                                           gTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                          >  1:364800/1‑100 (MQ=255)
                                            tttAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGtttt                                                         >  1:1165123/1‑100 (MQ=255)
                                             ttAAAGATAACGGGACCACCTAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTc                                                        <  2:1518103/100‑1 (MQ=255)
                                              tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTcc                                                       >  2:1866741/1‑100 (MQ=255)
                                              tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTcc                                                       <  1:1299906/100‑1 (MQ=255)
                                              tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTcc                                                       <  2:2526153/100‑1 (MQ=255)
                                              tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTcc                                                       <  2:131556/100‑1 (MQ=255)
                                               aaaaTTAACTGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCt                                                      >  2:2254697/1‑100 (MQ=255)
                                               aaaaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCt                                                      <  2:909411/100‑1 (MQ=255)
                                                aaaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCtt                                                     >  1:459871/1‑100 (MQ=255)
                                                 aaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTc                                                    <  2:1772703/100‑1 (MQ=255)
                                                   ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAg                                                  <  2:994648/100‑1 (MQ=255)
                                                     aaCGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAACGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAAGCGCGAGGTTTTCCTTCAg                                                  >  2:2190550/1‑98 (MQ=255)
                                                        ggCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAg                                                  <  2:2508902/95‑1 (MQ=255)
                                                         gCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCaa                                            >  2:1947624/1‑100 (MQ=255)
                                                           accaccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAAcc                                          >  1:109788/1‑100 (MQ=255)
                                                             caccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACcgc                                        <  2:1792999/100‑1 (MQ=255)
                                                             caccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACcgc                                        >  1:1334377/1‑100 (MQ=255)
                                                               ccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACcgcg                                       >  1:785425/1‑99 (MQ=255)
                                                                cGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGAc                                     <  2:740127/100‑1 (MQ=255)
                                                                cGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGAc                                     <  2:737973/100‑1 (MQ=255)
                                                                cGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGAc                                     <  2:1168722/100‑1 (MQ=255)
                                                                 gAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACa                                    >  1:1830349/1‑100 (MQ=255)
                                                                    aTTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATAc                                 >  1:1157364/1‑100 (MQ=255)
                                                                     ttCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACt                                <  1:2295556/100‑1 (MQ=255)
                                                                      tCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACt                                <  1:2294592/99‑1 (MQ=255)
                                                                        aGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGCCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTa                              <  2:1342155/99‑1 (MQ=255)
                                                                          tAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaat                           <  2:120149/100‑1 (MQ=255)
                                                                           naGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaata                          >  2:1996623/2‑100 (MQ=255)
                                                                           aaGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaata                          <  1:2169906/100‑1 (MQ=255)
                                                                            aGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaataa                         >  1:2057073/1‑100 (MQ=255)
                                                                             agtaataGTCAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaataat                        >  2:120801/7‑100 (MQ=37)
                                                                                 ataGTCAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCtt                    >  1:412738/3‑100 (MQ=255)
                                                                                  aaGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTc                   <  1:593958/100‑1 (MQ=255)
                                                                                  aaGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTc                   >  1:2070816/1‑100 (MQ=255)
                                                                                  aaGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTc                   <  2:67967/100‑1 (MQ=255)
                                                                                    gTCAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTg                 >  1:780323/1‑100 (MQ=255)
                                                                                     tCAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGt                >  1:891580/1‑100 (MQ=255)
                                                                                      cAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGCTCAACCGCGACATACTTAAATAATGCTTATGTc               >  2:323003/1‑100 (MQ=38)
                                                                                      cAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc               >  1:906898/1‑100 (MQ=255)
                                                                                      cAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc               >  1:1528470/1‑100 (MQ=255)
                                                                                       aaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCg              >  2:104340/1‑100 (MQ=255)
                                                                                       aaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCg              >  2:434530/1‑100 (MQ=255)
                                                                                        aaaGCCTCCGTTCGGAACCTTTCGTCTATTACTCAACAGTTAAGGCGCGAGGCTTTACTTCAGGTTCAACCGCGACATACTTAAATAATGCTTCTGTCGc             >  2:738479/1‑100 (MQ=255)
                                                                                        aaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGCCATACTTAAATAATGCTTCTGTCGc             >  2:293312/1‑100 (MQ=255)
                                                                                        aaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGc             >  1:482221/1‑100 (MQ=255)
                                                                                                cGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGc     >  2:1933598/1‑100 (MQ=255)
                                                                                                   ccGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGtt  >  1:1794325/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTT  >  NC_000913/1311732‑1311930

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