Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,995,556 C→A 6.5% pseudogene (205/447 nt) ygeN ← pseudogene, orgB family, part of T3SS PAI ETT2 remnant

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,995,5560CA 6.5% 597.8 / 21.6 193pseudogene (205/447 nt)ygeNpseudogene, orgB family, part of T3SS PAI ETT2 remnant
Reads supporting (aligned to +/- strand):  ref base C (86/94);  new base A (7/6);  total (93/100)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

Only 100 of 195 total aligned reads displayed.

CAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAATTTTGATAAAAACATCAACTTCT  >  NC_000913/2995408‑2995704
                                                                                                                                                    |                                                                                                                                                    
cAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGAc                                                                                                                                                      <  2:1494588/149‑1 (MQ=255)
  gATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACAt                                                                                                                                                    >  1:466449/1‑149 (MQ=255)
   aTTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACCAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACAtt                                                                                                                                                   <  1:161538/149‑1 (MQ=255)
        gTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCata                                                                                                                                              <  2:1309465/149‑1 (MQ=255)
         ttAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCatat                                                                                                                                             <  1:905657/149‑1 (MQ=255)
         ttAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCatat                                                                                                                                             >  2:261837/1‑149 (MQ=255)
          tAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCTGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATAtt                                                                                                                                            <  2:1219441/149‑1 (MQ=255)
          tAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATAtt                                                                                                                                            <  2:476268/149‑1 (MQ=255)
            aCCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGAATGTCCATAttac                                                                                                                                          <  2:379335/149‑3 (MQ=255)
            aCCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGAATGTCCATAttac                                                                                                                                          <  2:298480/149‑3 (MQ=255)
             ccTTGTCAGGGGAGAAACAATGATTATTAATTAAATCTTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGAt                                                                                                                                         <  2:850474/149‑1 (MQ=255)
             ccTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGAt                                                                                                                                         >  1:1543442/1‑149 (MQ=255)
             ccTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCAATATTGAt                                                                                                                                         >  1:1091048/1‑149 (MQ=255)
             ccTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTAACAATTATCCACAAATTCCTGCGGACTAAATTAAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGAt                                                                                                                                         <  1:1031909/149‑1 (MQ=255)
             ccTTGTCAGGGGAGAAACAATGATTATAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGAATGTCCATAttacat                                                                                                                                        <  2:1458474/149‑5 (MQ=255)
                tGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGAATGTCCATAttacattt                                                                                                                                      <  1:885087/149‑7 (MQ=255)
                   cAGGGGTGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTAAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGAAATACAAAATTGCTTTGTCg                                                                                                                                   >  2:120698/1‑149 (MQ=255)
                         agaAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACaaaa                                                                                                                             <  1:112535/149‑1 (MQ=255)
                          gaAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAat                                                                                                                            <  2:819838/149‑1 (MQ=255)
                          gaAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAat                                                                                                                            <  1:751182/149‑1 (MQ=255)
                           aaaCAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAata                                                                                                                           >  2:1058872/1‑149 (MQ=255)
                              cAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGa                                                                                                                        <  2:960752/149‑1 (MQ=255)
                              cAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGa                                                                                                                        <  2:396404/149‑1 (MQ=255)
                              cAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGa                                                                                                                        <  2:804747/149‑1 (MQ=255)
                                    ttattaATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACAt                                                                                                                                                    >  1:155491/1‑115 (MQ=255)
                                    ttattaATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACAt                                                                                                                                                    <  2:155491/115‑1 (MQ=255)
                                    ttattaATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTAc                                                                                                                  >  1:1252528/1‑149 (MQ=255)
                                      attaattaAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACtt                                                                                                                <  1:596493/149‑1 (MQ=255)
                                      attaattaAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTAGACAAAATATGACTCTACtt                                                                                                                <  1:86250/149‑1 (MQ=255)
                                       ttaattaaATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTc                                                                                                               >  2:667037/1‑149 (MQ=255)
                                         aattaaATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCaa                                                                                                             >  2:581193/1‑149 (MQ=255)
                                         aattaaATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCaa                                                                                                             <  1:849042/149‑1 (MQ=255)
                                         aattaaATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCaa                                                                                                             <  1:1171968/149‑1 (MQ=255)
                                         aattaaATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCaa                                                                                                             <  2:1316258/149‑1 (MQ=255)
                                          attaaATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAag                                                                                                            <  1:986011/149‑1 (MQ=255)
                                           ttaaATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCGATATTGATTTGTCGACAAAATATGACTCTACTTCAaga                                                                                                           <  1:1470092/149‑1 (MQ=255)
                                              aaTATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGAc                                                                                                        >  1:1414352/1‑149 (MQ=255)
                                                     cTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCtt                                                                                                 <  2:659546/149‑1 (MQ=255)
                                                     cTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCtt                                                                                                 >  1:1453456/1‑149 (MQ=255)
                                                      tCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCttt                                                                                                <  2:463159/149‑1 (MQ=255)
                                                        acaATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGa                                                                                              <  2:141432/149‑1 (MQ=255)
                                                        acaATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGa                                                                                              <  1:1427473/149‑1 (MQ=255)
                                                        acaATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGa                                                                                              <  2:1119657/149‑1 (MQ=255)
                                                         caATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGaa                                                                                             <  1:1239926/149‑1 (MQ=255)
                                                          aaTTATCCACAAATTCCTGCGGACTAAATTAAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGaaa                                                                                            >  2:885279/1‑149 (MQ=255)
                                                           aTTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAAc                                                                                           >  2:1535586/1‑149 (MQ=255)
                                                                ccACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTggta                                                                                      <  2:895030/149‑1 (MQ=255)
                                                                ccACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTggta                                                                                      <  2:925674/149‑1 (MQ=255)
                                                                  acaAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTggtagg                                                                                    >  2:1480949/1‑149 (MQ=255)
                                                                   caAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTggtaggt                                                                                   >  1:1044381/1‑149 (MQ=255)
                                                                   aaaaTTCCTTCTGTCTAAATTCAGCAATAAATTGATCGGTAAAAAAGACAAAACGCTTATCATCATGAAATGTAATATGGACTTTCCATATTGATTTCTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTggtaggt                                                                                   <  2:1436605/148‑1 (MQ=255)
                                                                   aaaaTTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTggtaggt                                                                                   <  1:1113892/148‑1 (MQ=255)
                                                                          cTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGga                                                                            >  2:1278933/1‑149 (MQ=255)
                                                                               gACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGAttt                                                                       >  2:621798/1‑149 (MQ=255)
                                                                                  tAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTTGTAGGTATATGGAGATTTACa                                                                    >  1:1541206/1‑149 (MQ=255)
                                                                                  tAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACa                                                                    >  1:148825/1‑149 (MQ=255)
                                                                                   aaaTTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAg                                                                   <  1:660943/149‑1 (MQ=255)
                                                                                   aaaTTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAg                                                                   <  2:592138/149‑1 (MQ=255)
                                                                                    aaTTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGTTAGGTATATGGAGATTTACatt                                                                  >  1:1138422/1‑147 (MQ=255)
                                                                                    aaTTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAgg                                                                  >  1:709330/1‑149 (MQ=255)
                                                                                    aaTTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAgg                                                                  >  1:509961/1‑149 (MQ=255)
                                                                                    aaTTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAgg                                                                  >  1:607788/1‑149 (MQ=255)
                                                                                       tCAGAAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGAcc                                                               <  1:302480/149‑1 (MQ=255)
                                                                                          gCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGa                                                            <  2:1167678/149‑1 (MQ=255)
                                                                                           cAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGat                                                           >  1:307300/1‑149 (MQ=255)
                                                                                            aaTAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGata                                                          <  1:1369584/149‑1 (MQ=255)
                                                                                            aaTAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGata                                                          >  2:569142/1‑149 (MQ=255)
                                                                                              tAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATAtt                                                        >  1:1058866/1‑149 (MQ=255)
                                                                                              tAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATAtt                                                        >  1:376697/1‑149 (MQ=255)
                                                                                                aaTTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTCGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGa                                                      <  1:346919/149‑1 (MQ=255)
                                                                                                  ttGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATg                                                    <  2:277454/149‑1 (MQ=255)
                                                                                                  ttGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATCCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATg                                                    <  2:193322/149‑1 (MQ=255)
                                                                                                        ggTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAAtt                                              <  1:1208746/149‑1 (MQ=255)
                                                                                                        ggTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAAtt                                              <  1:1160711/149‑1 (MQ=255)
                                                                                                        ggTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACCCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAAtt                                              >  2:962483/1‑149 (MQ=255)
                                                                                                          tAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAAtttt                                            <  2:1041266/149‑1 (MQ=255)
                                                                                                          tAAAAAACACAAAACGCTTATCATCATGAAAAGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAAATCGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAAtttt                                            >  1:19291/1‑149 (MQ=255)
                                                                                                              aaaCACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTa                                        <  2:338510/149‑1 (MQ=255)
                                                                                                                acacAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTAcc                                      <  2:1134648/149‑1 (MQ=255)
                                                                                                                   caAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACACaa                                   >  2:746164/1‑149 (MQ=255)
                                                                                                                    aaaaCGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAAc                                  >  1:1196642/1‑149 (MQ=255)
                                                                                                                    aaaaCGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACATCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTAAATGACCTGATATTGATGGTAATTTTGTTACCCAAc                                  >  1:95442/1‑149 (MQ=255)
                                                                                                                                atcatGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAAttt                      <  1:1278933/149‑1 (MQ=255)
                                                                                                                                 tcatGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAAtttt                     >  1:429341/1‑149 (MQ=255)
                                                                                                                                  catGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAAttttt                    >  2:1599215/1‑148 (MQ=255)
                                                                                                                                   atGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAATTTTGa                   >  1:1287643/1‑149 (MQ=255)
                                                                                                                                    tGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAATTTTGAt                  <  1:621798/149‑1 (MQ=255)
                                                                                                                                    tGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAATTTTGAt                  >  2:115715/1‑149 (MQ=255)
                                                                                                                                    tGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTAACCAACTTTCAAAAATTTTGAt                  >  2:1401305/1‑149 (MQ=255)
                                                                                                                                     gAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTTATTTTGTTACCCAACTTTCAAAAATTTTGATa                 >  2:248286/1‑149 (MQ=255)
                                                                                                                                      aaaTGTAATATTGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAATTTTGATaa                <  1:524353/149‑1 (MQ=255)
                                                                                                                                       aaatcaaTATGGAATGTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAATTTTGATaaa               >  2:1081738/6‑149 (MQ=255)
                                                                                                                                       aaTGTTATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATTTGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAATTTTGATaaa               <  2:1570411/149‑1 (MQ=255)
                                                                                                                                       aaTGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAATTTTGATaaa               >  1:174271/1‑149 (MQ=255)
                                                                                                                                        aTGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAAATGGTAGGTATATGGAGATTTACACGACCTGATATTGATGGTAATTTTGTTACACAACTATCAAAAATTTTTATaaaa              >  2:98821/1‑149 (MQ=255)
                                                                                                                                          tcaaTATGGAATGTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAATTTTGATAAAAAc            >  1:695778/3‑149 (MQ=255)
                                                                                                                                           tAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAATTTTGATAAAAACa           >  1:531463/1‑149 (MQ=255)
                                                                                                                                             atatGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAATTTTGATAAAAACATc         <  2:314050/149‑1 (MQ=255)
                                                                                                                                                tGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAATTTTGATAAAAACATCAAc      <  2:346357/149‑1 (MQ=255)
                                                                                                                                                    cATTCCATATTGATTTGTCGACAAAACATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAATTTTGATAAAAACATCAActtct  >  2:1330788/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCACAATTATCCACAAATTCCTGCGGACTAAATTCAGCAATAAATTGATCGGTAAAAAACACAAAACGCTTATCATCATGAAATGTAATATGGACATTCCATATTGATTTGTCGACAAAATATGACTCTACTTCAAGAGACTTATCTTTGAAACTGGTAGGTATATGGAGATTTACAGGACCTGATATTGATGGTAATTTTGTTACCCAACTTTCAAAAATTTTGATAAAAACATCAACTTCT  >  NC_000913/2995408‑2995704

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