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

Only 100 of 105 total aligned reads displayed.

TCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACGGAG  >  NC_000913/1470460‑1470654
                                                                                                  |                                                                                                
tCAAAACACCGAAAGAGCTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTTAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGa                                                                                                 >  2:1763455/1‑100 (MQ=255)
 cAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTaaa                                                                                                >  1:1142046/1‑100 (MQ=255)
  naaaCACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTaaaa                                                                                               >  2:1395843/2‑100 (MQ=255)
  aaaaCACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTaaaa                                                                                               >  2:222331/1‑100 (MQ=255)
   aaaCACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAaaaa                                                                                              >  2:40690/1‑100 (MQ=255)
   aaaCACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAaaaa                                                                                              >  2:347808/1‑100 (MQ=255)
    aaCACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAAt                                                                                             <  1:2272369/100‑1 (MQ=255)
       acCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATaca                                                                                          <  2:2263585/100‑1 (MQ=255)
       acCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATaca                                                                                          <  1:2186300/100‑1 (MQ=255)
        ccGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATacaa                                                                                         >  2:1650434/1‑100 (MQ=255)
        ccGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATacaa                                                                                         >  2:1673167/1‑100 (MQ=255)
         cGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTTAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATacaa                                                                                         <  1:794130/99‑1 (MQ=255)
          gAAAGATATAATTTATAGTGGTGTTGCATTGACAGATTTAATCTACACTGTTGAAATTCATAATTTATGTTTTGCTTTTGTTACTGGTAAAAATacaa                                                                                         <  2:2032240/98‑1 (MQ=25)
          gAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATacaaca                                                                                       >  1:2771776/1‑100 (MQ=255)
          gAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATacaac                                                                                        <  2:2115997/99‑1 (MQ=255)
           aaaGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAg                                                                                      <  1:2488122/100‑1 (MQ=255)
            aaGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAgg                                                                                     <  1:1074599/100‑1 (MQ=255)
                gaTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGtata                                                                                 >  1:1174044/1‑100 (MQ=255)
                 aTAATTGAAAGGGGTGTTGTATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATaa                                                                                <  1:1778548/100‑1 (MQ=255)
                   aaTTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAAt                                                                              >  2:2509995/1‑100 (MQ=255)
                   aaTTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATATAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAAt                                                                              >  2:2084699/1‑100 (MQ=255)
                    aTTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATa                                                                             >  2:1393807/1‑100 (MQ=255)
                    aTTGAAAGGAGGGTTGCATTGACAGAATGAATCTACACTGTTGAAATTCATCAGTTAGGGATTGCTTTTGTTGCTGGTAAAAAAAAAATAGGAAAAAATa                                                                             >  2:1296302/1‑100 (MQ=16)
                     ttGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAg                                                                            <  2:150010/100‑1 (MQ=255)
                       gAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTg                                                                          <  2:419594/100‑1 (MQ=255)
                       gAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTg                                                                          >  1:650658/1‑100 (MQ=255)
                       gAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTg                                                                          >  1:2503836/1‑100 (MQ=255)
                        aaaGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTgg                                                                         <  1:347811/100‑1 (MQ=255)
                         aaGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGc                                                                        >  1:1816427/1‑100 (MQ=255)
                         aaGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGc                                                                        <  1:1118097/100‑1 (MQ=255)
                           ggggTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCAc                                                                      >  2:885824/1‑100 (MQ=255)
                            gggTGTTGCATTGACAGATTGAATCTACACTTTTGAAATGCGGAATTGAGTTTTTGCTTTTGTTACTGGTTAAAATATAACAAGTATAAATAGTGGCACg                                                                     <  2:2587172/100‑1 (MQ=37)
                            gggTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACg                                                                     >  1:207049/1‑100 (MQ=255)
                              gtgtTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGAt                                                                   >  2:221425/1‑100 (MQ=255)
                                gtTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGAt                                                                   >  1:741629/1‑98 (MQ=255)
                                   gCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGt                                                              >  1:2962453/1‑100 (MQ=255)
                                   gCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGt                                                              >  2:1207767/1‑100 (MQ=255)
                                   gCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGt                                                              <  1:545031/100‑1 (MQ=255)
                                    cATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACTATCTCGtt                                                             <  2:1802571/100‑1 (MQ=255)
                                    cATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGtt                                                             <  1:1341184/100‑1 (MQ=255)
                                    cATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGtt                                                             <  1:3009614/100‑1 (MQ=255)
                                     aTTGACATATTGAATCTACACTTTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGt                                                              <  1:1397154/98‑1 (MQ=255)
                                      ttGACATATTTAATCTACACTGTTTAAATTCATAATTTATGTTTTTCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttat                                                           <  1:878603/100‑1 (MQ=255)
                                      ttGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttat                                                           >  2:799122/1‑100 (MQ=255)
                                       tGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACTACAGGTATAAATAGTGGCACGATCTCGttatt                                                          <  1:2889951/100‑1 (MQ=255)
                                        gACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttatta                                                         >  1:2391214/1‑100 (MQ=255)
                                        gACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttatta                                                         <  2:2825294/100‑1 (MQ=255)
                                         aCAGATTGAATCTACATTTTTTTAATTCAGAATTTAGGTTTTGCTTTTGTTACTAGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttatta                                                         <  1:790326/99‑1 (MQ=255)
                                          cAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACa                                                       >  1:1406565/1‑100 (MQ=255)
                                            gATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACaaa                                                     >  2:3039618/1‑100 (MQ=255)
                                            gATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACaaa                                                     >  2:1410817/1‑100 (MQ=255)
                                                gAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATg                                                  >  2:2865303/1‑99 (MQ=255)
                                                gAATATGCACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATgg                                                 <  2:1822922/100‑1 (MQ=255)
                                                 aaTCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGt                                                <  1:1794006/100‑1 (MQ=255)
                                                  aTCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTa                                               <  2:1208741/100‑1 (MQ=255)
                                                     tACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAATAGCACGAACTCGTTATTACAAAATGGTaaa                                             >  1:755371/1‑99 (MQ=255)
                                                      acacTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGa                                           <  2:1308592/100‑1 (MQ=255)
                                                       cacTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGTTTAAg                                            >  2:1730536/1‑98 (MQ=255)
                                                       cacTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGAt                                          <  2:2658511/100‑1 (MQ=255)
                                                        acTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGAt                                          <  2:1669469/99‑1 (MQ=255)
                                                         cTGTTGAAATTTAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATcc                                        <  1:2477781/100‑1 (MQ=255)
                                                         cTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATcc                                        >  1:1739372/1‑100 (MQ=255)
                                                           gTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCgg                                      >  2:34166/1‑100 (MQ=255)
                                                            ttGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACTGGTATAAATAGTGTCACGATCTCGTTATTACAAAATTGTAAAGATc                                         >  2:1652653/1‑96 (MQ=255)
                                                             tGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCa                                    >  2:2626728/1‑100 (MQ=255)
                                                              gAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCAc                                   >  1:356248/1‑100 (MQ=255)
                                                              gAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCAc                                   >  1:2176731/1‑100 (MQ=255)
                                                                  ttCAGTATTTATGTTTTTGTTTTGTTACTGGTAAAAATACAACAGGTAGAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtc                               <  1:414600/100‑1 (MQ=255)
                                                                  ttCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtc                               >  1:875188/1‑100 (MQ=255)
                                                                     aGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTcc                            <  2:2339677/100‑1 (MQ=255)
                                                                     aGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTcc                            >  1:1303282/1‑100 (MQ=255)
                                                                        atttatttTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCAt                         <  1:1952766/100‑1 (MQ=255)
                                                                         tttAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGGGGCATGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCAt                         <  1:1040434/99‑1 (MQ=255)
                                                                          ttAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTg                       <  2:1953616/100‑1 (MQ=255)
                                                                           tAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGt                      <  1:1475516/100‑1 (MQ=255)
                                                                           tAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGt                      <  1:3024014/100‑1 (MQ=255)
                                                                             ggTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGttt                    >  2:2019164/1‑100 (MQ=255)
                                                                             ggTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGttt                    >  1:2356018/1‑100 (MQ=255)
                                                                              gTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGtttt                   >  2:1090457/1‑100 (MQ=255)
                                                                                tttGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTAc                 >  1:1885132/1‑100 (MQ=255)
                                                                                  tGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTg               >  1:133936/1‑100 (MQ=255)
                                                                                  tGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAATGGCACGATATCGTTATTACAAAATGGTAAAGAGCCGGCACCGTCTCCCATTGTTTTACTg               >  2:290474/1‑100 (MQ=255)
                                                                                    cTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc             >  2:1705963/1‑100 (MQ=255)
                                                                                    cTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc             <  1:2992221/100‑1 (MQ=255)
                                                                                     ttttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGct            <  2:90095/100‑1 (MQ=255)
                                                                                       ttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGct            >  2:1253451/1‑98 (MQ=255)
                                                                                        tGTTATTTGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctact         <  1:569540/100‑1 (MQ=255)
                                                                                         gTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAAATCCGGCACCGTCTCCCATTGTTTTACTGGctacta        >  1:1515865/1‑100 (MQ=255)
                                                                                          ttACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTaa       >  1:949869/1‑100 (MQ=255)
                                                                                           tACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTTCAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTATACTGGCTACTAAc      >  1:745695/1‑100 (MQ=255)
                                                                                            aCTGGTAAAAATACAACAGGTATAAATAGTGGCATGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACg     <  1:2457045/100‑1 (MQ=255)
                                                                                             ntGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCACATTGTTTTACTGGCTACTaa       >  2:1250420/2‑97 (MQ=255)
                                                                                             cTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACgg    >  2:30232/1‑100 (MQ=255)
                                                                                              tGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACgga   <  2:1441771/100‑1 (MQ=255)
                                                                                              tGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACgg    >  1:1743916/1‑99 (MQ=255)
                                                                                               ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggag  >  1:947802/1‑100 (MQ=255)
                                                                                               ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggag  >  2:2978252/1‑100 (MQ=255)
                                                                                               ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggag  >  2:1517949/1‑100 (MQ=255)
                                                                                               ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTAGCTACTAACggag  >  2:1641714/1‑100 (MQ=255)
                                                                                               ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGAACCGGCAACGTCTCCCATTGTTTTACTGGCTAATAACgga   >  2:1593681/1‑99 (MQ=255)
                                                                                                  |                                                                                                
TCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACGGAG  >  NC_000913/1470460‑1470654

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