Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,190,298 +C coding (613/2481 nt) yehB ← putative fimbrial usher protein YehB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,190,2941.C98.2% 1150.0 / 10.6 279coding (617/2481 nt)yehBputative fimbrial usher protein YehB
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base C (136/138);  total (137/142)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 309 total aligned reads displayed.

AGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAACGTACATATGTACTC  >  NC_000913/2190194‑2190389
                                                                                                      |                                                                                               
aGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCcaca                                                                                                  <  2:3327071/101‑1 (MQ=255)
aGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCcaca                                                                                                  >  2:3327331/1‑101 (MQ=255)
   aTCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACAccc                                                                                              >  2:3327337/1‑101 (MQ=255)
     cAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACAccccc                                                                                            >  1:3327339/1‑100 (MQ=255)
     cAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACAccccc                                                                                            <  2:3326970/101‑2 (MQ=255)
      aaaaaTATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACAccccct                                                                                           <  2:3326609/101‑3 (MQ=255)
      aaaaaTATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACAccccct                                                                                           >  2:3327342/1‑99 (MQ=255)
       aaaaTATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACAccccctg                                                                                          <  1:3326791/101‑4 (MQ=255)
       aaaaTATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACAccccctg                                                                                          <  1:3327003/101‑4 (MQ=255)
          atatCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGAt                                                                                        >  1:3327346/1‑101 (MQ=255)
              cGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGtt                                                                                    >  1:3327350/1‑101 (MQ=255)
               gCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTa                                                                                   >  1:3327354/1‑101 (MQ=255)
                cTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTAt                                                                                  >  2:3327356/1‑101 (MQ=255)
                  tGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTAttt                                                                                >  2:3327360/1‑101 (MQ=255)
                   gATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTg                                                                               >  1:3327362/1‑101 (MQ=255)
                     tgtgTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGtt                                                                             >  2:3327365/1‑101 (MQ=255)
                     tgtgTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGtt                                                                             <  1:3327032/101‑1 (MQ=255)
                        gtACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTa                                                                          >  1:3327370/1‑101 (MQ=255)
                          aCATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACt                                                                        >  1:3327374/1‑101 (MQ=255)
                          aCATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACt                                                                        >  2:3327373/1‑101 (MQ=255)
                           cATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTg                                                                       <  2:3327036/101‑1 (MQ=255)
                            atatCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGa                                                                      <  1:3326884/101‑1 (MQ=255)
                            atatCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGa                                                                      >  1:3327377/1‑101 (MQ=255)
                             tatCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGaa                                                                     <  2:3327085/101‑1 (MQ=255)
                             tatCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGaa                                                                     <  1:3326763/101‑1 (MQ=255)
                              atCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGaaa                                                                    >  1:3327380/1‑101 (MQ=255)
                                cACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACt                                                                  >  1:3327383/1‑101 (MQ=255)
                                cACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACt                                                                  <  2:3326991/101‑1 (MQ=255)
                                 aCCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTg                                                                 >  1:3327388/1‑101 (MQ=255)
                                    cACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCa                                                              >  2:3327390/1‑101 (MQ=255)
                                     aCGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCAt                                                             >  1:3327392/1‑101 (MQ=255)
                                     aCGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCAt                                                             <  1:3326842/101‑1 (MQ=255)
                                      cgcgAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATc                                                            <  1:3326905/101‑1 (MQ=255)
                                       gcgAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCa                                                           >  2:3327396/1‑101 (MQ=255)
                                        cgAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAg                                                          >  2:3327399/1‑101 (MQ=255)
                                         gAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGa                                                         >  1:3327403/1‑101 (MQ=255)
                                          aaGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGaa                                                        >  2:3327406/1‑101 (MQ=255)
                                          aaGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGaa                                                        <  1:3326774/101‑1 (MQ=255)
                                           aGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAAt                                                       <  2:3326957/101‑1 (MQ=255)
                                           aGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAAt                                                       >  1:3327407/1‑101 (MQ=255)
                                             cGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGc                                                     <  1:3326856/101‑1 (MQ=255)
                                              gTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGTGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACt                                                                  <  2:3327050/87‑1 (MQ=255)
                                                gCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGt                                                  >  1:3327413/1‑101 (MQ=255)
                                                gCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGt                                                  >  2:3327414/1‑101 (MQ=255)
                                                gCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGt                                                  <  1:3327100/101‑1 (MQ=255)
                                                 ccGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACGCCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGtt                                                 <  2:3326867/101‑1 (MQ=255)
                                                 ccGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGtt                                                 <  1:3326935/101‑1 (MQ=255)
                                                 ccGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGtt                                                 <  2:3326768/101‑1 (MQ=255)
                                                  cGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGa                                                         >  2:3327417/1‑92 (MQ=255)
                                                  cGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTg                                                <  2:3326917/101‑1 (MQ=255)
                                                    agaAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGcc                                              >  1:3327422/1‑101 (MQ=255)
                                                      aaGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCAc                                            >  2:3327424/1‑101 (MQ=255)
                                                      aaGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCAc                                            <  2:3327158/101‑1 (MQ=255)
                                                       aGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCAcc                                           <  1:3326976/101‑1 (MQ=255)
                                                       aGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCAc                                            <  2:3326934/100‑1 (MQ=255)
                                                        gTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCAccc                                          <  1:3326628/101‑1 (MQ=255)
                                                         ttGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCAcccc                                         >  2:3327428/1‑101 (MQ=255)
                                                           gggCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAg                                       >  1:3327431/1‑101 (MQ=255)
                                                           gggCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAg                                       <  2:3326949/101‑1 (MQ=255)
                                                           gggCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAg                                       <  1:3327023/101‑1 (MQ=255)
                                                            ggCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGt                                      <  2:3326846/101‑1 (MQ=255)
                                                              cAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTaa                                    >  2:3327433/1‑101 (MQ=255)
                                                               aaaTCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTaaa                                   >  2:3327437/1‑101 (MQ=255)
                                                                  tCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAAttt                                >  2:3327440/1‑101 (MQ=255)
                                                                   ccACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTa                               <  1:3327052/101‑1 (MQ=255)
                                                                    cACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTaa                              <  1:3326872/101‑1 (MQ=255)
                                                                    cACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTaa                              <  1:3327207/101‑1 (MQ=255)
                                                                    cACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTaa                              <  2:3326910/101‑1 (MQ=255)
                                                                    cACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTaa                              >  1:3327442/1‑101 (MQ=255)
                                                                     aCGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAAc                             <  1:3326719/101‑1 (MQ=255)
                                                                      cGTTCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAAcc                            <  2:3326804/101‑1 (MQ=255)
                                                                        ttCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCg                          <  2:3326878/101‑1 (MQ=255)
                                                                        ttCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCg                          <  1:3327039/101‑1 (MQ=255)
                                                                        ttCCAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCg                          >  1:3327446/1‑101 (MQ=255)
                                                                           cAGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTg                       <  1:3327160/101‑1 (MQ=255)
                                                                            aGATACAGGGTATTGCTTTTCCACGCCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGt                      <  2:3327127/101‑1 (MQ=255)
                                                                            aGATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGt                      <  1:3326810/101‑1 (MQ=255)
                                                                             gATACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGtt                     >  1:3327450/1‑101 (MQ=255)
                                                                               tACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTaa                   <  2:3327009/101‑1 (MQ=255)
                                                                               tACAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTaa                   >  2:3327451/1‑101 (MQ=255)
                                                                                aCAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTaaa                  <  1:3327331/101‑1 (MQ=255)
                                                                                aCAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTaaa                  <  1:3327182/101‑1 (MQ=255)
                                                                                 cAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTaaaa                 <  1:3326966/101‑1 (MQ=255)
                                                                                 cAGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTaaaa                 >  1:3327456/1‑101 (MQ=255)
                                                                                  aGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAAc                >  2:3327460/1‑101 (MQ=255)
                                                                                  aGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAAc                >  1:3327459/1‑101 (MQ=255)
                                                                                  aGGGTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAAc                <  1:3326984/101‑1 (MQ=255)
                                                                                    ggTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAACGt              >  1:3327464/1‑101 (MQ=255)
                                                                                    ggTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAACGt              <  1:3326944/101‑1 (MQ=255)
                                                                                     gTATTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAACGTa             <  1:3327086/101‑1 (MQ=255)
                                                                                       aTTGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAACGTACa           >  1:3327468/1‑101 (MQ=255)
                                                                                         tGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAACGTACata         <  2:3327187/101‑1 (MQ=255)
                                                                                         tGCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAACGTACata         >  2:3327470/1‑101 (MQ=255)
                                                                                          gCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAACGTACatat        >  1:3327473/1‑101 (MQ=255)
                                                                                          gCTTTTCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAACGTACatat        <  1:3327007/101‑1 (MQ=255)
                                                                                             tttCCACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAACGTACATATGTa     >  2:3327475/1‑101 (MQ=255)
                                                                                                ccACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAACGTACATATGTActc  <  2:3327107/101‑1 (MQ=255)
                                                                                                ccACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAACGTACATATGTActc  <  2:3327065/101‑1 (MQ=255)
                                                                                                ccACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAACGTACATATGTActc  >  1:3327481/1‑101 (MQ=255)
                                                                                                ccACACCCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCGGTTAAAACGTACATATGTActc  <  1:3326992/101‑1 (MQ=255)
                                                                                                      |                                                                                               
AGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCTGGATTGTTATTTGTTTTACTGAAACTGGCATCAGAATGCAGTTGCCACCCCAGTAAATTTAACCCGCTGTTAAAACGTACATATGTACTC  >  NC_000913/2190194‑2190389

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