Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,470,558 G→A pseudogene (42/3497 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,470,5580GA97.8% 442.4 / ‑0.8 137pseudogene (42/3497 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (54/78);  total (57/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

Only 100 of 140 total aligned reads displayed.

CAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACGGAGG  >  NC_000913/1470461‑1470655
                                                                                                 |                                                                                                 
cAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTaaa                                                                                                 <  2:2859196/100‑1 (MQ=255)
 aaaaCACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTaaaa                                                                                                <  2:2462042/100‑1 (MQ=255)
 aaaaCACCGAAAGAGAGAAGTGAAAGGGGGGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAAGTTAGGTTTTGCTTTTGTTACTGGTaaaa                                                                                                >  1:1157786/1‑100 (MQ=25)
  aaaCACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAaaaa                                                                                               <  2:1848683/100‑1 (MQ=255)
    acacCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATa                                                                                             <  1:129253/100‑1 (MQ=255)
      acCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATaca                                                                                           >  2:1832100/1‑100 (MQ=255)
       ccGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATacaa                                                                                          >  1:2673598/1‑100 (MQ=255)
       ccGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGT‑‑AAATacaaca                                                                                        >  1:320602/1‑100 (MQ=255)
          aaaGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAg                                                                                       <  2:2334962/100‑1 (MQ=255)
           aaGAGATAATTGATAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAgg                                                                                      <  1:567843/100‑1 (MQ=255)
            agagATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGt                                                                                     <  1:3408384/100‑1 (MQ=255)
             gagaTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGta                                                                                    <  1:1534989/100‑1 (MQ=255)
                 tAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATaaa                                                                                <  1:1015990/100‑1 (MQ=255)
                 tAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATaaa                                                                                >  2:614951/1‑100 (MQ=255)
                   aTTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATa                                                                              >  2:3023429/1‑100 (MQ=255)
                    ttGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAg                                                                             >  1:3083783/1‑100 (MQ=255)
                    ttGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAg                                                                             <  1:806953/100‑1 (MQ=255)
                     tGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGt                                                                            >  2:3055370/1‑100 (MQ=255)
                       aaaGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTgg                                                                          <  1:488701/100‑1 (MQ=255)
                           gggTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGt                                                                            <  1:1088197/94‑1 (MQ=255)
                            ggTGTTGCATTTACAGATTTAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACg                                                                      <  1:1394432/99‑1 (MQ=255)
                              tgtTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGAtc                                                                   <  2:522984/100‑1 (MQ=255)
                              tgtTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGAtc                                                                   <  2:3387250/100‑1 (MQ=255)
                               gtTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGAtct                                                                  >  2:206821/1‑100 (MQ=255)
                                  gCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGt                                                               <  1:678768/100‑1 (MQ=255)
                                   cATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGtt                                                              >  2:1851906/1‑100 (MQ=255)
                                    aTTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGtta                                                             <  1:202629/100‑1 (MQ=255)
                                    aTTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGtta                                                             >  2:1005203/1‑100 (MQ=255)
                                     ttGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttat                                                            >  2:2212481/1‑100 (MQ=255)
                                      tGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttatt                                                           <  1:1844046/100‑1 (MQ=255)
                                      tGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttatt                                                           <  1:3148448/100‑1 (MQ=255)
                                      tGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttatt                                                           >  1:3074026/1‑100 (MQ=255)
                                        aCAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTAc                                                         >  1:2160769/1‑100 (MQ=255)
                                         cAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACa                                                        <  2:821801/100‑1 (MQ=255)
                                           gATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGATTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACaaa                                                      >  1:1493633/1‑100 (MQ=255)
                                            aTTGAATCTACACTGTTGAAAGTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACCACAGGTATACATTGTGGCACGATCTCGTTATTACaaaa                                                     >  1:799807/1‑100 (MQ=255)
                                             ttGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAAt                                                    <  2:593940/100‑1 (MQ=255)
                                             ttGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAAt                                                    <  2:369968/100‑1 (MQ=255)
                                             ttGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAAt                                                    <  2:2842277/100‑1 (MQ=255)
                                              tGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATg                                                   >  1:1818230/1‑100 (MQ=255)
                                                 aTCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATggga                                                <  2:2348191/100‑3 (MQ=255)
                                                  tCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTa                                                <  2:2898004/99‑1 (MQ=255)
                                                    tACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTaaa                                              >  2:1901331/1‑99 (MQ=255)
                                                    tACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAg                                             >  1:2497224/1‑100 (MQ=255)
                                                    tACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAg                                             <  1:3394353/100‑1 (MQ=255)
                                                    tACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAg                                             <  2:1472022/100‑1 (MQ=255)
                                                    tACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAg                                             >  2:1349807/1‑100 (MQ=255)
                                                      cacTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGAt                                           >  1:27689/1‑100 (MQ=255)
                                                      cacTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTCAAGAt                                           >  2:1506782/1‑100 (MQ=255)
                                                      cacTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGAt                                           <  2:1708661/100‑1 (MQ=255)
                                                      cacTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGAt                                           >  2:1131088/1‑100 (MQ=255)
                                                          gTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCgg                                       <  1:542002/100‑1 (MQ=255)
                                                             gAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCAc                                    >  1:497377/1‑100 (MQ=255)
                                                              aaaTTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCAcc                                   <  1:3088243/100‑1 (MQ=255)
                                                               aaTTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCg                                  >  1:1080848/1‑100 (MQ=255)
                                                                aTTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGt                                 >  2:3223867/1‑100 (MQ=255)
                                                                   cAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtctc                              >  1:2403064/1‑100 (MQ=255)
                                                                   cAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtctc                              <  1:2508965/100‑1 (MQ=255)
                                                                    aGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTcc                             <  2:1692182/100‑1 (MQ=255)
                                                                    aGAAATTAGGTTAAGCATTTATTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtctc                              <  1:987139/99‑1 (MQ=255)
                                                                     gAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTccc                            >  2:133449/1‑100 (MQ=255)
                                                                      aaTTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCa                           <  1:36273/100‑1 (MQ=255)
                                                                       aTTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTGTAAATAGTGGCACGATCTCGGTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCAt                          >  2:1896318/1‑100 (MQ=255)
                                                                       aTTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCAt                          >  2:2802442/1‑100 (MQ=255)
                                                                       aTTTAGGTTTTACTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCAt                          <  1:2580443/100‑1 (MQ=255)
                                                                        tttAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCAtt                         <  2:1142153/100‑1 (MQ=255)
                                                                        tttAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCAtt                         <  1:13023/100‑1 (MQ=255)
                                                                         cttccgatctcTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTg                        <  2:2216703/90‑1 (MQ=255)
                                                                            ggTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGttt                     <  1:1221288/100‑1 (MQ=255)
                                                                            ggTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGttt                     <  1:2321142/100‑1 (MQ=255)
                                                                             gTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGtttt                    >  2:809914/1‑100 (MQ=255)
                                                                             gTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGtttt                    >  1:2167117/1‑100 (MQ=255)
                                                                               tttGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTAc                  >  1:1616304/1‑100 (MQ=255)
                                                                               tttGCTTTTGTTACTGGGAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATGACAAAATGGGAAAGATCCGGCACCGTCTCCCATTGTTTTAc                  >  1:640163/1‑100 (MQ=255)
                                                                                ttGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACt                 <  1:1568640/100‑1 (MQ=255)
                                                                                ttGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACt                 <  2:32060/100‑1 (MQ=255)
                                                                                ttGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACt                 >  1:3339852/1‑100 (MQ=255)
                                                                                 tGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTg                <  1:1895175/100‑1 (MQ=255)
                                                                                 tGCTTTTGTTACTGGTAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTgg               >  1:2702967/1‑100 (MQ=255)
                                                                                  gCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTgg               >  2:3205103/1‑100 (MQ=255)
                                                                                   cTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc              <  1:634465/100‑1 (MQ=255)
                                                                                   cTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc              <  1:872084/100‑1 (MQ=255)
                                                                                   cTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc              <  2:1604932/100‑1 (MQ=255)
                                                                                    ttttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGct             <  2:3359792/100‑1 (MQ=255)
                                                                                    ttttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGct             <  1:2271346/100‑1 (MQ=255)
                                                                                     tttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGcta            <  1:1940306/100‑1 (MQ=255)
                                                                                     tttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGct             <  1:761892/99‑1 (MQ=255)
                                                                                      ttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctac           <  2:1966353/100‑1 (MQ=255)
                                                                                      ttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctac           <  1:846341/100‑1 (MQ=255)
                                                                                      ttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctac           <  1:388611/100‑1 (MQ=255)
                                                                                       tGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctact          >  2:3087351/1‑100 (MQ=255)
                                                                                          tACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAAc       >  2:3039455/1‑100 (MQ=255)
                                                                                          tACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAAc       >  2:647899/1‑100 (MQ=255)
                                                                                          tACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAAc       <  2:1132545/100‑1 (MQ=255)
                                                                                          tACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAAc       <  1:824238/100‑1 (MQ=255)
                                                                                          tACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGAGCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAAc       >  2:2435010/1‑100 (MQ=255)
                                                                                              ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggag   >  1:1341178/1‑100 (MQ=255)
                                                                                              ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggag   >  2:2302221/1‑100 (MQ=255)
                                                                                              ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggag   >  1:2944218/1‑100 (MQ=255)
                                                                                              ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTACTGGCTACTAACggagg  >  2:198310/1‑100 (MQ=255)
                                                                                                 |                                                                                                 
CAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACGGAGG  >  NC_000913/1470461‑1470655

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