Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,467,477 A→C pseudogene (2086/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,467,4770AC96.5% 686.5 / 5.9 203pseudogene (2086/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base C (96/100);  total (100/103)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

TGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTG  >  NC_000913/1467383‑1467573
                                                                                              |                                                                                                
tGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAAc                                                                                             <  2:2467612/100‑1 (MQ=255)
tGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAAc                                                                                             <  1:4098929/100‑1 (MQ=255)
tGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAAc                                                                                             >  2:3282304/1‑100 (MQ=255)
   gAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGtt                                                                                          <  1:4915701/100‑1 (MQ=255)
      tACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTAc                                                                                       <  2:2146915/100‑1 (MQ=255)
      tACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTAc                                                                                       <  1:2399515/100‑1 (MQ=255)
      tACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTAc                                                                                       >  1:2361844/1‑100 (MQ=255)
      tACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTAc                                                                                       <  2:2942452/100‑1 (MQ=255)
       aCAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACg                                                                                      >  2:2009278/1‑100 (MQ=255)
         aGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGCCAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGc                                                                                    <  2:2515002/100‑1 (MQ=255)
         aGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGc                                                                                    <  1:4577050/100‑1 (MQ=255)
           tGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAg                                                                                  <  2:234512/100‑1 (MQ=255)
            gAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGt                                                                                 <  2:1519784/100‑1 (MQ=255)
            gAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGt                                                                                 >  2:3532084/1‑100 (MQ=255)
              atatCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGa                                                                               <  1:4686487/100‑1 (MQ=255)
                atCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGata                                                                             >  1:4237722/1‑100 (MQ=255)
                atCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGata                                                                             >  1:5573288/1‑100 (MQ=255)
                 tCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGataa                                                                            >  1:3131346/1‑100 (MQ=255)
                    cTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGAtaataa                                                                         <  1:708959/100‑1 (MQ=255)
                    cTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGAtaataa                                                                         <  1:3535362/100‑1 (MQ=255)
                    cTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGAtaataa                                                                         <  2:38481/100‑1 (MQ=255)
                     tGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAAc                                                                        <  1:3416989/100‑1 (MQ=255)
                     tGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAAc                                                                        <  2:209481/100‑1 (MQ=255)
                     tGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAAc                                                                        >  1:900827/1‑100 (MQ=255)
                      ggTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACg                                                                       <  1:4764897/100‑1 (MQ=255)
                      ggTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACg                                                                       <  1:3916973/100‑1 (MQ=255)
                       gTAACGTTCTGGTTGATACGGCTATAACCGCAGACAATGAGGCGGTATATTTCTTCGATACATCCGTGGGTCTAAAAGTTTCCGGCAGTGATCCTAACGt                                                                      >  2:1356617/1‑100 (MQ=255)
                       gTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGt                                                                      <  2:2827071/100‑1 (MQ=255)
                       gTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGt                                                                      >  2:4851920/1‑100 (MQ=255)
                       gTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGt                                                                      <  1:1824819/100‑1 (MQ=255)
                       gTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGt                                                                      >  2:1495462/1‑100 (MQ=255)
                        tAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGATCTCAACGTTTACGGCAGTGATAATAACGTg                                                                     <  1:4460725/100‑1 (MQ=255)
                         aaCGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGa                                                                    <  1:3993255/100‑1 (MQ=255)
                         aaCGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGa                                                                    <  1:1514618/100‑1 (MQ=255)
                          aCGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGAc                                                                   >  1:4842564/1‑100 (MQ=255)
                           cGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACg                                                                  >  2:2892276/1‑100 (MQ=255)
                            gTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGc                                                                 >  1:2460167/1‑100 (MQ=255)
                             ttCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCt                                                                >  1:538855/1‑100 (MQ=255)
                             ttCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCt                                                                <  1:2919808/100‑1 (MQ=255)
                             tgcTGGTGGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCt                                                                <  2:4088901/98‑1 (MQ=255)
                              tCTGGTTGACACGGATCCAACCGCAGACCCTGTGGCGGAATACTTATAAGATCCAACCGTGGGCATCAACGTTTCCGGCAGTGATAATAACGTGACGCtg                                                               <  2:5492415/100‑1 (MQ=255)
                               cTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCtgg                                                              <  2:2062571/100‑1 (MQ=255)
                                 ggTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCtggat                                                            >  2:1267553/1‑100 (MQ=255)
                                 ggTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCtggat                                                            >  2:3124998/1‑100 (MQ=255)
                                 ggTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCtggat                                                            >  1:1063129/1‑100 (MQ=255)
                                 ggTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCtggat                                                            >  2:2228124/1‑100 (MQ=255)
                                 ggTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCtggat                                                            >  1:3019140/1‑100 (MQ=255)
                                 ggTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCtggat                                                            >  1:3628438/1‑100 (MQ=255)
                                 ggTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTAAGAATAACGTGACGCtggat                                                            >  1:5219832/1‑100 (MQ=255)
                                  gTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCtggatg                                                           <  1:1011004/100‑1 (MQ=255)
                                  gTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCtggatg                                                           <  1:438375/100‑1 (MQ=255)
                                   ttGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCtggatgg                                                          <  1:2279613/100‑1 (MQ=255)
                                    tGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCtggatgga                                                         >  1:2087908/1‑100 (MQ=255)
                                       tAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAg                                                      >  2:4869431/1‑100 (MQ=255)
                                        aaGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGt                                                     <  1:4319350/100‑1 (MQ=255)
                                          ggATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATGATAACGTGACGCTGGATGGAAAGTTa                                                   <  2:4900102/100‑1 (MQ=255)
                                            aTAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAAc                                                 >  2:1652783/1‑100 (MQ=255)
                                            aTAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAAc                                                 >  1:87956/1‑100 (MQ=255)
                                               aaaCCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACtgt                                              >  2:4264012/1‑100 (MQ=255)
                                               aaaCCGCAGACAATGCGGCGGAATATTTCTTCGAACCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACTTGACGCTGGATGGAAAGTTAACtgt                                              >  1:1989548/1‑100 (MQ=255)
                                                 aCCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACtgttg                                            >  1:1304344/1‑100 (MQ=255)
                                                 aCCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACtgttg                                            <  2:4683112/100‑1 (MQ=255)
                                                 aCCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACtgttg                                            <  2:3250160/100‑1 (MQ=255)
                                                 aCCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACtgttg                                            <  1:782715/100‑1 (MQ=255)
                                                  ccGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACtgttgt                                           >  1:1595697/1‑100 (MQ=255)
                                                  ccGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACtgttgt                                           >  2:710785/1‑100 (MQ=255)
                                                   cGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTa                                          >  2:2797946/1‑100 (MQ=255)
                                                    gCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAACGTTAACTGTTGTAt                                         >  2:3936873/1‑100 (MQ=255)
                                                    gCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTAt                                         >  1:2308273/1‑100 (MQ=255)
                                                     cAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATc                                        <  2:397496/100‑1 (MQ=255)
                                                      aGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCa                                       <  1:5504757/100‑1 (MQ=255)
                                                       gACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCan                                      <  1:5081208/100‑2 (MQ=255)
                                                       gACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAg                                      >  2:2948399/1‑100 (MQ=255)
                                                       gACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAg                                      <  1:1388482/100‑1 (MQ=255)
                                                         cAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCGGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGAc                                    <  1:4128112/100‑1 (MQ=255)
                                                         cAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGAc                                    >  2:5219878/1‑100 (MQ=255)
                                                           aTGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTTATCGTTTGCGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAg                                  <  1:3666332/100‑1 (MQ=255)
                                                           aTGCGGCGGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAg                                  <  1:4177835/100‑1 (MQ=255)
                                                           aTGCGGCGGAACATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAg                                  <  1:1443536/100‑1 (MQ=255)
                                                             gcggcgGAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTg                                <  2:4697943/100‑1 (MQ=255)
                                                              cggcggAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGa                               <  2:3465030/100‑1 (MQ=255)
                                                              cggcggAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGa                               <  1:3313061/100‑1 (MQ=255)
                                                                gcggAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAgg                             <  1:5003024/100‑1 (MQ=255)
                                                                 cggAATATTTCTTCGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGt                            >  2:4358494/1‑100 (MQ=255)
                                                                    aaTGTTTCTTCGAGCCAACCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTAc                         <  2:2652847/100‑1 (MQ=255)
                                                                          tcttcGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGCCAGTGAGGTTACTTCTCg                   >  2:5257093/1‑100 (MQ=255)
                                                                          tcttcGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCg                   <  1:2757939/100‑1 (MQ=255)
                                                                           cttcGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGt                  >  2:4548924/1‑100 (MQ=255)
                                                                           cttcGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGt                  >  2:5102661/1‑100 (MQ=255)
                                                                           cttcGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGt                  <  1:1727014/100‑1 (MQ=255)
                                                                              cGATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCag               <  2:347307/100‑1 (MQ=255)
                                                                               gATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCaga              >  2:3805354/1‑100 (MQ=255)
                                                                               gATCCATCCGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCaga              >  2:5346931/1‑100 (MQ=255)
                                                                                   catccGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTaa          <  1:5173678/100‑1 (MQ=255)
                                                                                    atccGTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTaac         >  2:4425406/1‑99 (MQ=255)
                                                                                       cgctggTCTCAACGTTTACGGCAGTGATAATAACGTGTCGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTTattta      <  2:1393525/96‑1 (MQ=255)
                                                                                        gTGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAatttat     <  2:1640377/100‑1 (MQ=255)
                                                                                         tGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAatttatt    <  1:1374673/100‑1 (MQ=255)
                                                                                         tGGGTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAatttatt    >  2:3430092/1‑100 (MQ=255)
                                                                                           ggTCTCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTg  <  2:1574925/100‑1 (MQ=255)
                                                                                              |                                                                                                
TGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTG  >  NC_000913/1467383‑1467573

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