Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,381,915 G→T 15.1% intergenic (+13/‑32) ycjV → / → ompG pseudogene; putative ATP‑binding component of a transport system/outer membrane porin G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,381,9150GT15.1% 411.1 / 56.2 153intergenic (+13/‑32)ycjV/ompGpseudogene; putative ATP‑binding component of a transport system/outer membrane porin G
Reads supporting (aligned to +/- strand):  ref base G (66/63);  new base T (17/6);  total (83/70)
Fisher's exact test for biased strand distribution p-value = 6.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01

Only 100 of 160 total aligned reads displayed.

ATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGA  >  NC_000913/1381817‑1382014
                                                                                                  |                                                                                                   
aTCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAgg                                                                                                    <  1:3297171/100‑1 (MQ=255)
aTCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAg                                                                                                     >  1:1023971/1‑99 (MQ=255)
   aTGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCAt                                                                                                 >  2:2372224/1‑100 (MQ=255)
    tGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATACGCTAAATACAGGGGGAAGGCAtt                                                                                                >  2:2170673/1‑100 (MQ=255)
     gCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTc                                                                                               >  1:1909486/1‑100 (MQ=255)
      cAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTcc                                                                                              >  1:1063971/1‑100 (MQ=255)
        ggAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGCGGGAAGGCATTcccc                                                                                            >  2:464932/1‑100 (MQ=255)
         gagaAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTccccc                                                                                           <  2:2165897/100‑1 (MQ=255)
          agaAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCa                                                                                          <  2:285986/100‑1 (MQ=255)
          agaAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCa                                                                                          >  2:1632723/1‑100 (MQ=255)
           gaAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAg                                                                                         <  2:1174552/100‑1 (MQ=255)
           gaAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAg                                                                                         >  1:2644779/1‑100 (MQ=255)
           gaAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAg                                                                                         <  1:2985295/100‑1 (MQ=255)
           gaAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAg                                                                                         <  1:3873111/100‑1 (MQ=255)
           gaAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCACAg                                                                                         >  1:1516731/1‑100 (MQ=255)
            aaaaTATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCGTAAACGGAAATAGCAATTCGCTAAATACAGGGGGAATGCCTTccccctg                                                                                         <  2:3550374/99‑3 (MQ=255)
            aaaaTATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAATGCCTTccccctgt                                                                                        <  1:2046203/100‑4 (MQ=255)
             aaaTATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAATGCCTTccccctgta                                                                                       <  1:87479/100‑5 (MQ=255)
             aaaTATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGa                                                                                       >  2:3924186/1‑100 (MQ=255)
              aaTATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACTGGGGGAATGCCTTccccctgtat                                                                                      <  1:2072062/100‑6 (MQ=255)
              aaTATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGat                                                                                      <  2:3031530/100‑1 (MQ=255)
               atatCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGata                                                                                     <  1:1002098/100‑1 (MQ=255)
                tatCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGataa                                                                                    >  1:687354/1‑100 (MQ=255)
                tatCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGataa                                                                                    >  1:4051135/1‑100 (MQ=255)
                      tATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAg                                                                              >  2:879014/1‑100 (MQ=255)
                      tATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAg                                                                              <  1:2066888/100‑1 (MQ=255)
                      tATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAg                                                                              <  1:2563471/100‑1 (MQ=255)
                      tATTCATTTTGATATCACGAAATGTCATTTATTTGATGCAGAAACGGAAATAGCAATTCGCTAAATCCAGGGGGAAGGCATTCCCCCAGGATAATACAAg                                                                              <  2:3356366/100‑1 (MQ=255)
                       aTTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAgg                                                                             >  1:1734016/1‑100 (MQ=255)
                       aTTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAg                                                                              <  2:72054/99‑1 (MQ=255)
                       aTTCAATTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAAGCATTCCCCAAGGATAATACAAgg                                                                             <  2:1981485/100‑1 (MQ=255)
                           aTTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACa                                                                         >  1:1476569/1‑100 (MQ=255)
                                gATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACaataat                                                                    >  1:3210327/1‑100 (MQ=255)
                                  tatGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGa                                                                  <  1:919538/100‑1 (MQ=255)
                                   atGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGaa                                                                 >  1:763681/1‑100 (MQ=255)
                                     gACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGaaaa                                                               >  2:3289115/1‑100 (MQ=255)
                                      aCGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGaaaaa                                                              >  1:2729907/1‑100 (MQ=255)
                                         aaaTGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGtt                                                           <  2:687353/100‑1 (MQ=255)
                                         aaaTGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGtt                                                           <  1:1347818/100‑1 (MQ=255)
                                          aaTGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGtta                                                          >  2:3767453/1‑100 (MQ=255)
                                           aTGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCGGGGTAATACAAGGAACAATAATGAAAAAGttat                                                         <  2:2752414/100‑1 (MQ=255)
                                           aTGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGttat                                                         <  1:282112/100‑1 (MQ=255)
                                           aTGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGttat                                                         <  2:2528221/100‑1 (MQ=255)
                                            tGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGttatt                                                        >  1:2230214/1‑100 (MQ=255)
                                              tCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTAc                                                      <  2:4051108/100‑1 (MQ=255)
                                                aTTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTAccc                                                    >  1:1790411/1‑100 (MQ=255)
                                                aTTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTAccc                                                    >  1:2304116/1‑100 (MQ=255)
                                                aTTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTAccc                                                    <  1:2269105/100‑1 (MQ=255)
                                                     tttGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTAc                                               >  2:3357089/1‑100 (MQ=255)
                                                     tttGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTAc                                               <  2:1524933/100‑1 (MQ=255)
                                                      ttGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTAcc                                              >  1:861428/1‑100 (MQ=255)
                                                      ttGATGCAGAAACGGAAATAGCAATTCGATAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTAcc                                              <  1:689567/100‑1 (MQ=255)
                                                        gATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGc                                            >  1:1114164/1‑100 (MQ=255)
                                                        gATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGc                                            >  1:3676156/1‑100 (MQ=255)
                                                           gCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACt                                         <  1:3431685/100‑1 (MQ=255)
                                                           gCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGACATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACt                                         >  1:3423301/1‑100 (MQ=255)
                                                            cAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTg                                        >  1:2372995/1‑100 (MQ=255)
                                                            cAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTg                                        >  2:1226316/1‑100 (MQ=255)
                                                             aGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTgg                                       <  1:2904305/100‑1 (MQ=255)
                                                               aaaCGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTg                                     <  1:800118/100‑1 (MQ=255)
                                                                aaCGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGa                                    <  1:1142588/100‑1 (MQ=255)
                                                                  cGGTAATAGCAATTCGCTAAAGACAGAGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGAtg                                  <  2:1200229/100‑1 (MQ=255)
                                                                    gAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGAtgtg                                <  2:147502/100‑1 (MQ=255)
                                                                    gAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGAtgt                                 >  2:1795097/1‑99 (MQ=255)
                                                                          gCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCggg                          <  2:1489480/100‑1 (MQ=255)
                                                                          gCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCggg                          <  1:2427188/100‑1 (MQ=255)
                                                                          gCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCggg                          >  1:39245/1‑100 (MQ=255)
                                                                           cAATTCGCTAAATACAATGGGAAGGCATTTCTCCAGTATAATACAAGGAACAATCATGAAAAAGTTATTACCCAGTACCGAACTTGTGATGTGTGCGGGa                         <  2:1034195/100‑1 (MQ=255)
                                                                                cGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGc                    >  1:984687/1‑100 (MQ=255)
                                                                                cGCTAAATACAGGGGGAAGGCATCCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGc                    >  1:239176/1‑100 (MQ=255)
                                                                                 gCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGcc                   >  1:3149126/1‑100 (MQ=255)
                                                                                  cTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCt                  <  1:812417/100‑1 (MQ=255)
                                                                                  cTAAAAACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCt                  <  2:3664152/100‑1 (MQ=255)
                                                                                     tatcctgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTgcg               >  2:1247711/7‑100 (MQ=255)
                                                                                       tcctgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCATgcgc              >  2:494202/5‑99 (MQ=255)
                                                                                       tACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGca             <  1:3794599/100‑1 (MQ=255)
                                                                                       tACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGca             >  2:3169298/1‑100 (MQ=255)
                                                                                        cctgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGcac            >  2:1365274/4‑100 (MQ=255)
                                                                                        cctgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGcac            >  1:1709475/4‑100 (MQ=255)
                                                                                        cctgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGcac            >  2:1922196/4‑100 (MQ=255)
                                                                                         ctgggggAATGCTTTCCCCCAGGATAATACAAGGCACCATCATGAAATCGTTATCACTCTGGACGGCACTGGCGATGTGTGCGGGTATGCCCTgc                >  2:233325/3‑95 (MQ=255)
                                                                                         ctgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGcaca           >  2:2848206/3‑100 (MQ=255)
                                                                                         ctgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGcaca           >  1:34630/3‑100 (MQ=255)
                                                                                         ctgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGcaca           >  1:4143315/3‑100 (MQ=255)
                                                                                         ctgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGcaca           >  1:980472/3‑100 (MQ=255)
                                                                                         ctgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGcaca           >  2:4016433/3‑100 (MQ=255)
                                                                                         cAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGcaca           >  2:1713195/1‑100 (MQ=255)
                                                                                          tgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTATACCGCACTGGTGATCTGTGCGGAAATAGCCTGCGcaca           >  1:1456795/2‑99 (MQ=255)
                                                                                          aGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAg          >  2:3338872/1‑100 (MQ=255)
                                                                                           gggggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAgg         >  1:541792/1‑100 (MQ=255)
                                                                                           gggggAAGGCATTACCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAgg         >  2:2592643/1‑100 (MQ=255)
                                                                                             ggaaaGGCACTCCCCCAGGACAATACAAGGAACAATAATGAAAAAGTGATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGcc       <  2:1390665/97‑1 (MQ=255)
                                                                                               gAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGa     >  1:2072565/1‑100 (MQ=255)
                                                                                                aaGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGccatt    <  1:1277519/100‑4 (MQ=255)
                                                                                                aaGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGag    >  1:2003883/1‑100 (MQ=255)
                                                                                                aaGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGag    <  1:2564775/100‑1 (MQ=255)
                                                                                                 aGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGagg   >  2:2650507/1‑100 (MQ=255)
                                                                                                  ggCATTCCCCCAGGGTAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGccattcc  <  2:2996789/100‑6 (MQ=255)
                                                                                                  ggCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGccattcc  <  1:279837/100‑6 (MQ=255)
                                                                                                  ggCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGccattcc  <  1:2642985/100‑6 (MQ=255)
                                                                                                  |                                                                                                   
ATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGA  >  NC_000913/1381817‑1382014

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