Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,991,334 T→A 10.3% I61N (ATT→AAT) ‡ cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,991,3340TA10.3% 832.0 / 68.6 252I61N (ATT→AAT) ‡cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (113/113);  new base A (12/14);  total (125/127)
Fisher's exact test for biased strand distribution p-value = 8.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 257 total aligned reads displayed.

CTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGG  >  NC_000913/3991186‑3991476
                                                                                                                                                    |                                                                                                                                              
      gCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCt                                                                                                                                          <  1:1065654/149‑1 (MQ=255)
         aTAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTa                                                                                                                                       >  3:660891/1‑149 (MQ=255)
            aaTCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTacac                                                                                                                                    >  1:1071185/1‑149 (MQ=255)
             aTCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACg                                                                                                                                   >  4:413377/1‑149 (MQ=255)
                  ttGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtga                                                                                                                              <  1:1340520/149‑1 (MQ=255)
                   tGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtgat                                                                                                                             <  1:1146211/149‑1 (MQ=255)
                      gtgtGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGc                                                                                                                                            <  1:80000/131‑1 (MQ=255)
                      gtgtGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGc                                                                                                                                            <  1:862023/131‑1 (MQ=255)
                      gtgtGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGaa                                                                                                                          >  1:1250170/1‑149 (MQ=255)
                       tgtgGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGaaa                                                                                                                         <  1:878386/149‑1 (MQ=255)
                       cgtgGATCGCGCGCTTGCTGCTATCGGGCCTGCATTCCAACAGGTCTACACTCTACTTCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGaaa                                                                                                                         >  2:377382/2‑149 (MQ=255)
                                     ttGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                       <  3:997331/121‑2 (MQ=255)
                                       gctgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCt                                                                                                         <  3:383598/149‑1 (MQ=255)
                                         tgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGa                                                                                                       >  2:180215/1‑149 (MQ=255)
                                          gctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                                                                                                      <  3:63765/149‑1 (MQ=255)
                                          gctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                                                                                                      <  3:183826/149‑1 (MQ=255)
                                               ggggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGc                                                                                                 <  2:1071185/149‑1 (MQ=255)
                                                  gCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTg                                                                                              <  4:1466622/149‑1 (MQ=255)
                                                  gCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTg                                                                                              >  3:1342894/1‑149 (MQ=255)
                                                    cTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGaa                                                                                            <  3:413377/149‑1 (MQ=255)
                                                    cTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGaa                                                                                            <  3:1293011/149‑1 (MQ=255)
                                                     tGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTTTCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGcc                                                                                                                >  3:884025/1‑128 (MQ=255)
                                                     tGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGcc                                                                                                                <  4:884025/128‑1 (MQ=255)
                                                     tGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAAc                                                                                           >  3:151087/1‑149 (MQ=255)
                                                     tGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttgg                                                                                                                                      <  2:428863/106‑3 (MQ=255)
                                                     tGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                       <  1:1386851/105‑2 (MQ=255)
                                                     tGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                       <  1:1429493/105‑2 (MQ=255)
                                                     tGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                       >  2:1429493/1‑104 (MQ=255)
                                                         ttCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtgat                                                                                                                             >  2:648439/1‑111 (MQ=255)
                                                          tCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTAt                                                                                      <  1:377382/149‑1 (MQ=255)
                                                              aCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTg                                                                                  >  2:787980/1‑149 (MQ=255)
                                                              aCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCAACCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTg                                                                                  <  3:783781/149‑1 (MQ=255)
                                                                  gTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAAt                                                                              >  2:544663/1‑149 (MQ=255)
                                                                     tACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                       <  3:1200947/89‑2 (MQ=255)
                                                                     tACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                       <  1:819862/89‑2 (MQ=255)
                                                                     tACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                       >  2:819862/1‑88 (MQ=255)
                                                                        aGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGt                                                                        <  3:1014072/149‑1 (MQ=255)
                                                                            tACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACAACCTGAACCAGCTTGAACTGTATCGTGGAATGTCAGTACAg                                                                    >  4:1198042/1‑149 (MQ=255)
                                                                                 ccGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt                                                                                   <  4:1372466/129‑1 (MQ=255)
                                                                                 ccGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt                                                                                   >  3:1372466/1‑129 (MQ=255)
                                                                                 ccGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATcc                                                               >  4:1169248/1‑149 (MQ=255)
                                                                                  cGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATccg                                                              >  4:857016/1‑149 (MQ=255)
                                                                                    aCATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATccgcc                                                            <  2:264242/149‑1 (MQ=255)
                                                                                         gttgCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATccgc                                                             <  2:893001/143‑1 (MQ=255)
                                                                                           tgCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGt                                                     <  1:180215/149‑1 (MQ=255)
                                                                                            gCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAg                                                                                                  >  3:583125/1‑103 (MQ=255)
                                                                                            gCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGt                                                                                        >  1:917331/1‑113 (MQ=255)
                                                                                            gCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGt                                                                                        <  2:917331/113‑1 (MQ=255)
                                                                                             cACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAAt                                                                              <  4:647874/122‑1 (MQ=255)
                                                                                               cTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGc                                                 >  1:1150245/1‑149 (MQ=255)
                                                                                                 aTCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGaa                                                         >  1:562437/1‑139 (MQ=255)
                                                                                                   cACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTcc                                             >  4:1206624/1‑149 (MQ=255)
                                                                                                    aCCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGc                                                                                                 >  3:777096/1‑96 (MQ=255)
                                                                                                     ccATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTa                                                                                       <  4:525487/105‑1 (MQ=255)
                                                                                                     ccATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTc                                                                           >  2:1394506/1‑117 (MQ=255)
                                                                                                      cATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCa                                            >  1:635743/1‑147 (MQ=255)
                                                                                                        tCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTa                                        >  2:245293/1‑149 (MQ=255)
                                                                                                        tCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTa                                        <  3:857016/149‑1 (MQ=255)
                                                                                                         ccGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                                                                                                      >  3:954032/1‑86 (MQ=255)
                                                                                                         ccGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                                                                                                      <  4:954032/86‑1 (MQ=255)
                                                                                                         ccGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTAc                                       <  1:207336/149‑1 (MQ=255)
                                                                                                           gCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCa                                                                          <  1:502270/112‑1 (MQ=255)
                                                                                                           gCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                     <  2:1150245/149‑1 (MQ=255)
                                                                                                             tAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGgt                                   <  4:1342894/149‑1 (MQ=255)
                                                                                                               aTGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGaaa                                                        <  4:1057619/126‑1 (MQ=255)
                                                                                                                   cGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATcc                                                               <  3:268814/115‑1 (MQ=255)
                                                                                                                   cGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATcc                                                               >  4:268814/1‑115 (MQ=255)
                                                                                                                       ttACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCa                         <  1:1456971/149‑1 (MQ=255)
                                                                                                                         aCCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATg                       >  4:181456/1‑149 (MQ=255)
                                                                                                                          ccTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCagta                      <  2:574352/149‑4 (MQ=255)
                                                                                                                           cTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtga                                                                                                                              >  1:302150/1‑44 (MQ=255)
                                                                                                                           cTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGCTGAAACTCAAc                                                                                                                   >  4:1306838/1‑55 (MQ=255)
                                                                                                                           cTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGaa                                                                                            <  1:1016708/78‑1 (MQ=255)
                                                                                                                            ttGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAg                                                                    <  3:684264/101‑1 (MQ=255)
                                                                                                                            ttGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGt                                                     <  2:1181042/116‑1 (MQ=255)
                                                                                                                            ttGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGt                                                     >  1:1181042/1‑116 (MQ=255)
                                                                                                                            ttGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGgc                    >  4:458498/1‑149 (MQ=255)
                                                                                                                                tGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGa                                                                                             <  2:836362/72‑1 (MQ=255)
                                                                                                                                 ggTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGgcagca                >  4:1078208/1‑148 (MQ=255)
                                                                                                                                 ggTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGgcagca                <  3:1078208/148‑1 (MQ=255)
                                                                                                                                  gTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGaaa                                                                                                                         >  2:1258087/1‑42 (MQ=255)
                                                                                                                                      cGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGt          <  1:648683/149‑1 (MQ=255)
                                                                                                                                       gTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCaa                                           >  4:790828/1‑115 (MQ=255)
                                                                                                                                       gTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTc         >  1:1290651/1‑149 (MQ=255)
                                                                                                                                       gTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTc         >  4:368150/1‑149 (MQ=255)
                                                                                                                                        ttCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTGCAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTcc        >  4:702424/1‑149 (MQ=255)
                                                                                                                                        ttCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTcc        <  4:671687/149‑1 (MQ=255)
                                                                                                                                          cccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTAc                                       >  2:765558/1‑116 (MQ=255)
                                                                                                                                          cccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTAc                                       <  1:765558/116‑1 (MQ=255)
                                                                                                                                           ccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGt                                                     <  4:782779/101‑1 (MQ=255)
                                                                                                                                           ccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTc                                              <  1:396890/108‑1 (MQ=255)
                                                                                                                                           ccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTgg                                    >  3:1172892/1‑118 (MQ=255)
                                                                                                                                            cAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGaa                                                                               <  3:314629/74‑1 (MQ=255)
                                                                                                                                            cAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGaa                                                                               >  4:314629/1‑74 (MQ=255)
                                                                                                                                            cAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                     >  2:1313679/1‑116 (MQ=255)
                                                                                                                                            cAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAg    >  3:794899/1‑149 (MQ=255)
                                                                                                                                             aaaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGt                                                                            <  3:450297/76‑1 (MQ=255)
                                                                                                                                              aaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAggg  >  1:1030011/1‑149 (MQ=255)
                                                                                                                                              aaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAggg  <  2:1030011/149‑1 (MQ=255)
                                                                                                                                                    tttGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGgca                   >  3:1105392/1‑126 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGG  >  NC_000913/3991186‑3991476

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