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,5580GA99.1% 404.0 / ‑4.5 114pseudogene (42/3497 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (47/65);  total (48/66)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 116 total aligned reads displayed.

CAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACGGAG  >  NC_000913/1470461‑1470654
                                                                                                 |                                                                                                
cAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTaaa                                                                                                <  2:231928/100‑1 (MQ=255)
 aaaaCACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTaaaa                                                                                               >  2:752088/1‑100 (MQ=255)
  aaaCACCGAAAGAGATAATTGAACGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAaaaa                                                                                              >  1:397204/1‑100 (MQ=255)
  aaaCACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAaaaa                                                                                              >  1:194639/1‑100 (MQ=255)
    acacCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATa                                                                                            <  1:1907407/100‑1 (MQ=255)
    acacCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGGAAAAATa                                                                                            <  2:1027925/100‑1 (MQ=255)
     cacCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATac                                                                                           >  2:800098/1‑100 (MQ=255)
     cacCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATac                                                                                           >  2:797635/1‑100 (MQ=255)
      acCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATaca                                                                                          <  2:1229192/100‑1 (MQ=255)
      acCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATaca                                                                                          <  1:793434/100‑1 (MQ=255)
        cGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATacaac                                                                                        >  1:1059249/1‑100 (MQ=255)
         gcaaGAGATAATTAAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATacaaca                                                                                       >  2:1309154/3‑100 (MQ=255)
           aaGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAgg                                                                                     <  2:1409057/100‑1 (MQ=255)
            agagATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGt                                                                                    >  2:599866/1‑100 (MQ=255)
             gagaTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGta                                                                                   >  1:767942/1‑100 (MQ=255)
               gaTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGtata                                                                                 >  2:1457001/1‑100 (MQ=255)
               gaTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGtata                                                                                 >  1:70641/1‑100 (MQ=255)
                 tAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATaaa                                                                               >  2:543703/1‑100 (MQ=255)
                  aaTTGAAAGGGGTGTTGCATTGACAGATTGAATCTACATTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAAt                                                                              <  1:370169/100‑1 (MQ=255)
                  aaTTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAAt                                                                              <  1:22750/100‑1 (MQ=255)
                   aTTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATa                                                                             >  2:623511/1‑100 (MQ=255)
                     tGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGt                                                                           <  1:1252121/100‑1 (MQ=255)
                      gAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTg                                                                          <  1:918383/100‑1 (MQ=255)
                      gAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTg                                                                          >  1:836468/1‑100 (MQ=255)
                        aaGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGc                                                                        >  2:952683/1‑100 (MQ=255)
                        aaGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGc                                                                        <  2:1100991/100‑1 (MQ=255)
                        aaGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGc                                                                        <  1:1853579/100‑1 (MQ=255)
                           gggTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACg                                                                     <  1:1585734/100‑1 (MQ=255)
                             gtgtTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCAc                                                                      <  1:1881419/97‑1 (MQ=255)
                                 tGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCg                                                               >  1:838589/1‑100 (MQ=255)
                                 tGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCg                                                               >  1:1871960/1‑100 (MQ=255)
                                 tGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTCACGTGTAAACATACAACAGGTATAAATAGTGGCACGATCTCg                                                               >  1:1363243/1‑100 (MQ=255)
                                  gCATTGACAGATTGAGTCTACAATGTTGAACTTCGGAATTTAGGTTTTGCTTTTGTTCCTGGGACACATACAACGAGTATACATAGTgtccgcctct                                                                 >  1:1545647/1‑88 (MQ=255)
                                    aTTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGtta                                                            >  1:1643473/1‑100 (MQ=255)
                                    aTTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGtta                                                            <  1:1069483/100‑1 (MQ=255)
                                     ttGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttat                                                           <  1:364467/100‑1 (MQ=255)
                                     ttGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttat                                                           <  2:353673/100‑1 (MQ=255)
                                     ttGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttat                                                           <  2:888299/100‑1 (MQ=255)
                                     ttGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttat                                                           >  1:1127861/1‑100 (MQ=255)
                                     ttGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttat                                                           <  1:797633/100‑1 (MQ=255)
                                     ttGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGATACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttat                                                           <  1:800096/100‑1 (MQ=255)
                                       gACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttatta                                                         >  1:1403071/1‑100 (MQ=255)
                                       gACAGATTGAATCAACACTGTTGAAATTCAGAATTTAGCTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttatta                                                         >  1:1199530/1‑100 (MQ=255)
                                        aCAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTAc                                                        <  1:599864/100‑1 (MQ=255)
                                        aCAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTAc                                                        <  1:870689/100‑1 (MQ=255)
                                        aCAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTAc                                                        <  2:739698/100‑1 (MQ=255)
                                           gATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACaaa                                                     >  1:1050473/1‑100 (MQ=255)
                                            aTTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACaaaa                                                    >  2:724241/1‑100 (MQ=255)
                                            aTTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACaaa                                                     >  2:911100/1‑99 (MQ=255)
                                              tGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAAt                                                   <  1:957846/99‑1 (MQ=255)
                                              tGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATg                                                  >  1:44374/1‑100 (MQ=255)
                                                aaTCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGt                                                >  2:613024/1‑100 (MQ=255)
                                                 aTCTACTCTGCTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTa                                               <  2:1816646/100‑1 (MQ=255)
                                                   cTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTaaa                                             <  1:998183/100‑1 (MQ=255)
                                                   cTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTaaa                                             <  1:752087/100‑1 (MQ=255)
                                                   cTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTaaa                                             <  1:682636/100‑1 (MQ=255)
                                                   cTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTaaa                                             <  1:295559/100‑1 (MQ=255)
                                                    tACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAg                                            <  1:1556161/100‑1 (MQ=255)
                                                    tACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAg                                            <  1:968389/100‑1 (MQ=255)
                                                    tACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAg                                            <  2:1189515/100‑1 (MQ=255)
                                                    tACACTGTTGAAATTCAGAATTTAGGTTTTGCGTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAg                                            >  2:1842727/1‑100 (MQ=255)
                                                     acacTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGa                                           <  2:194640/100‑1 (MQ=255)
                                                      cacTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGAt                                          <  1:1673277/100‑1 (MQ=255)
                                                         tGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCg                                       <  2:105822/100‑1 (MQ=255)
                                                          gTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCgg                                      <  2:834703/100‑1 (MQ=255)
                                                            tGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCa                                    >  2:1158228/1‑100 (MQ=255)
                                                            tGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCa                                    >  1:1096675/1‑100 (MQ=255)
                                                             gAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCAc                                   <  2:896843/100‑1 (MQ=255)
                                                             gAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCAc                                   >  1:1650896/1‑100 (MQ=255)
                                                              aaaTTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCAcc                                  >  1:1942266/1‑100 (MQ=255)
                                                                 ttCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtc                               <  1:1515624/100‑1 (MQ=255)
                                                                  tCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtct                              <  1:308204/100‑1 (MQ=255)
                                                                  tCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtct                              >  1:1698510/1‑100 (MQ=255)
                                                                    agagTTCAGGTTTTGCTTTTGTTACTGGTCAAAATACAACAGGTATAAATAGTGGCACGATCTCGTCATTACAAAATGGTAAAGACCCGGCACCGTCTcc                            <  1:1648639/100‑1 (MQ=255)
                                                                    aGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAACGATCCGGCACCGTCTcc                            >  1:1528131/1‑100 (MQ=255)
                                                                     taaTTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTGCAAAATGGTAAAGATCCGGCACCGTCTccn                           <  1:1722924/99‑2 (MQ=255)
                                                                     gAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTccc                           <  1:1415231/100‑1 (MQ=255)
                                                                     gAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTccc                           <  2:650027/100‑1 (MQ=255)
                                                                     gAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTccc                           <  1:489243/100‑1 (MQ=255)
                                                                      aaTTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTccaa                          >  1:1842039/1‑98 (MQ=255)
                                                                      aaTTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCa                          <  2:842340/100‑1 (MQ=255)
                                                                       aTTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCAt                         <  2:1942279/100‑1 (MQ=255)
                                                                        tttAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCAtt                        <  1:173169/100‑1 (MQ=255)
                                                                        tttAGGTTTTGCTTTTGTTACTCGTAAAAATCCAACATGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCAtt                        <  1:1309147/100‑1 (MQ=255)
                                                                                 tGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTg               <  1:226418/100‑1 (MQ=255)
                                                                                  gCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTgg              <  1:656814/100‑1 (MQ=255)
                                                                                   cTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc             <  2:310020/100‑1 (MQ=255)
                                                                                   cTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc             >  2:1559556/1‑100 (MQ=255)
                                                                                   cTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc             <  1:1494667/100‑1 (MQ=255)
                                                                                    ttttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGct            <  2:1130636/100‑1 (MQ=255)
                                                                                        gTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctacta        >  2:1108254/1‑100 (MQ=255)
                                                                                        gTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctacta        <  1:1523954/100‑1 (MQ=255)
                                                                                         ttACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTaa       >  2:968395/1‑100 (MQ=255)
                                                                                          tACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAAc      <  1:1017653/100‑1 (MQ=255)
                                                                                           aCTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACg     <  1:911099/100‑1 (MQ=255)
                                                                                           aCTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACg     <  1:681535/100‑1 (MQ=255)
                                                                                             tGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACgga   <  1:1523798/100‑1 (MQ=255)
                                                                                             tGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACgga   <  1:1102918/100‑1 (MQ=255)
                                                                                              ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggag  >  2:525765/1‑100 (MQ=255)
                                                                                              ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggag  >  1:133358/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACGGAG  >  NC_000913/1470461‑1470654

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